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Abstract 

It is thought that a region of pseudo- vacuum, where the chiral order parameter is 
misaligned from its vacuum orientation in isospin space, might occasionally form in 
high energy hadronic or nuclear collisions. The possible detection of such disoriented 
chiral condensate (DCC) would provide useful information about the chiral structure 
of the QCD vacuum and/or the chiral phase transition of strong interactions at high 
temperature. We review the theoretical developments concerning the possible DCC 
formation in high-energy collisions as well as the various experimental searches that 
have been performed so far. We discuss future prospects for upcoming DCC searches, 
e.g. in high-energy heavy-ion collision experiments at RHIC and LHC. 
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1 Introduction 



In very high energy hadronic and/or nuclear collisions highly excited states 
are produced and subsequently decay toward vacuum via incoherent multi- 

* Corresponding Author. Tel: -^33 1 69 15 70 36, Fax: -^33 1 69 15 82 87 

Email addresses: bmohantySveccal . ernet . in (B. Mohanty), 
Julien.Serreau@th.u-psud.fr (J. Serreau). 

^ APC is unite mixte de recherche UMR 7164 (CNRS, Universite Paris 7, CEA, 
Observatoire de Paris). LPT is UMR 8627 (CNRS, Universtite Paris-Sud 11). 



Preprint submitted to Physics Reports 



2 February 2008 



particle emission. Due to the approximate SUr{2) xSUl{2) chiral symmetry of 
strong interactions, there exists a contimmm of nearly degenerate low-energy 
(pseudo- vacuum) states. These correspond to collective excitations where the 
chiral quark condensate is rotated from its vacuum orientation in chiral space 
a ~ (qq) and can be seen as semi-classical configurations of the pion field 
TT ~ {qT^^q) , where q = {u, d) denotes the light-quark doublet and r are the 
usual Pauli matrices. An interesting possibility is that the decay of highly ex- 
cited states produced in high-energy collisions may proceed via one of these 
collective states characterized by a disoriented chiral condensate (DCC). The 
latter would subsequently decay toward ordinary vacuum through coherent 
emission of low-momentum pions. Due to the semi-classical nature of the 
corresponding emission process, this may lead to specific signatures, such as 
anomalously large event-by-event fiuctuations of the charged-to-neutral ratio 
of produced pions. If the space-time region where this happens is large enough 
the phenomenon might be experimentally observable, thereby providing an in- 
teresting opportunity to study the chiral structure of QCD. 

Although very speculative, the idea that a DCC may form in high-energy 
collisions is quite appealing. Since it has been proposed in the early 1990's 
[1,2] (see also [3,4,5]) it has attracted a lot of interest and has generated an 
intense theoretical and experimental activity. One of the main motivations - 
beyond its appealing simplicity - has probably been the existence of exotic, 
so-called Centauro events reported in the cosmic ray literature [6,7,8], where 
clusters consisting of almost exclusively charged pions and no neutrals have 
been observed. The DCC would indeed provide a simple explanation for such 
phenomenon. 

On the theory side, a plausible mechanism for DCC formation in the context 
of high-energy heavy- ion collisions has been identified [9,10]: Due to the large 
energy deposit in the collision zone, a hot, chirally symmetric state (quark- 
gluon plasma) is formed. Due to the fast expansion at early times, the system 
is suddenly quenched down to the low-temperature phase, where chiral sym- 
metry is spontaneously broken. The subsequent far-from-equilibrium evolution 
triggers an exponential growth of long-wavelength pion modes, resulting in the 
formation of a strong (semi-classical) pion field configuration. This "quench 
scenario" has become a paradigm for DCC formation in heavy-ion collisions 
and has been widely used to investigate the phenomenological aspects of DCC 
production. As a side effect, this has led people to think more about non- 
equilibrium dynamics in the context of high-energy physics. In particular, this 
has triggered a number of theoretical developments concerning the descrip- 
tion of far-from-equilibrium quantum fields, ^ in connection with other areas 
of research, such as condensed-matter physics, or early-time cosmology (see 
e.g. Ref. [12]). 



^ For a recent review see e.g. Ref. [11]. 
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Numerous experimental searches have been performed in parallel with the 
development of theoretical ideas. These include the analysis of various cos- 
mic ray experiments [6,7], nucleon-nucleon collisions at CERN [13,14,15] and 
Fermilab [16], with, in particular, the dedicated MiniMAX experiment [17], as 
well as nucleus-nucleus collisions at the CERN SPS [18,19,20,21] and presently 
at RHIC [22,23]. The search for DCCs and other exotic events is part of the 
heavy-ions physics program to be performed by the multi-purpose detector 
ALICE at the LHC [24,25]. These investigations have led, in particular, to 
the development of powerful experimental tools to search for non-statistical 
fluctuations and/or to detect non-trivial structures in high-multiplicity events. 
No clear positive signal has been reported so far and upper bounds have been 
put on the likelihood of DCC formation, in particular, in heavy-ion collisions 
at SPS energies [18,19,21]. However, at the same time, it has been understood 
that the present experimental limit actually is consistent with theoretical ex- 
pectations based on the quench scenario [26,27]. 



This paper presents a status report of this field of research. In the first part, we 
review the main theoretical ideas and developments relative to DCC physics, 
including DCC formation and evolution in the multi-particle environment as 
well as phenomenological aspects. The second part is devoted to the numerous 
experimental studies which have been performed to search for DCC signals in 
high-energy hadronic and/or nuclear collisions. We emphasize that this report 
is, by no means, an exhaustive review of the extensive literature on this topics. 
Our main aim is, instead, to provide a comprehensive synthesis of what we 
believe are the most relevant aspects of DCC physics, both theoretically and 
experimentally, in view of further investigations, in particular, with upcoming 
DCC searches at RHIC and LHC. We mention that earlier reviews can be 
found in [28,29,30,31] 



To end this introduction and open our review, let us cite Bj or ken's words in 
his "DCC trouble list" [32], which illustrate very well the spirit of the present 
report: 



Existence of DCC: 



Must it exist? NO 

Should it exist? MAYBE 

Might it exist? YES 

Does it exist? IT'S WORTH HAVE A LOOK 
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2 Theory 



2.1 The disoriented chiral condensate: basic ideas 

The theory of strong interactions exhibits an approximate chiral symmetry 
SUr{2) X SUl{2), which is spontaneously broken in the vacuum - or in thermal 
equilibrium at sufficiently low temperatures. The associated order parameter, 
namely the quark condensate {quqi)-, can be represented as a four- component 
vector ifa = (vf, a) transforming under the 0(4) subgroup of SUr{2) x SUl{2). 
In the vacuum, the order parameter points in the cr-direction. However, due 
to the approximate chiral symmetry, one might expect that under appropriate 
conditions, there could exist a region of space, separated from the physical 
vacuum for some period of time, where the order parameter develops a non- 
trivial pionic component: This is the disoriented chiral condensate. ^ 

2.1.1 Baked-Alaska 

Bjorken and collaborators [1,28,34] have put forward a very simple and intu- 
itive physical picture of the possible formation of a disoriented pseudo- vacuum 
state in hadronic coUisions. We start our report by reviewing this so-called 
baked-Alaska scenario, which provides a useful guide for physical intuition. 
Consider a high multiplicity collision event with a large transverse energy 
release but no high-pT jets. In this situation, the hadronization time can be 
rather long, as large as a few fm/ c. Prior to hadronization, the primary partons 
carry most of the released energy away from the collision point at essentially 
the speed of light. One imagines a thin, "hot" expanding shell which isolates 
the relatively "cold" interior from the outer vacuum. If the energy density left 
behind is low enough, the interior of the fireball should look very similar to the 
vacuum, with an associated quark condensate. However, if the time it takes to 
cook this baked-Alaska is short enough, the quark condensate in the interior 
might be rotated from its usual orientation since the energy density associated 
with the explicit breaking of chiral symmetry is small. ^ When the hot shell 
hadronizes, the disoriented interior comes into contact with the true vacuum 
and radiates away its pionic orientation, resulting in coherent emission of soft 

^ It is worth emphasizing that there is no contradiction with the well-known Vafa- 
Witten theorem [33], which states that no pionic component of the chiral condensate 
can develop in a stationary state: The DCC is, by essence, a transient phenomenon. 
Similarly, it should be stressed that the production of a DCC state does not contra- 
dict usual conservation laws, such as e.g. parity, isospin, or charge conservation, as 
one eventually has to average over all possible equivalent directions in isospin space 
to obtain physical results (see subsection Sec. 2.1.3 below). 
^ A rough order of magnitude is: e /^"^^ ~ 20 MeV/fm^. 
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pions, with strong isospin correlations. 



As a simple realization of these ideas [35] , one might represent the hot debris 
located on the surface of the fireball as a source for the long-wavelength pionic 
excitation associated with the disoriented interior. In the linear approximation, 
the dynamics of the latter can be described by the following equation: 

{n + mi)7r{x) = J{x), (1) 



where the arrows denote vectors in isospin-space. A simple example is, for 
instance, J{x) — J{t) Q{t — to) 5{t — r), corresponding to a spherically expand- 
ing shell with initial radius = to, where is the time where the expansion 
starts. After a typical decoupling (hadronization) time, the source J{t) van- 
ishes and the field excitation 7i{x) decays into freely propagating pions. In this 
simple model, pion emission is characterized by the coherent state: 

l^") = -A^ exp (j:^ J ^ J„(k) at (k)) |0) (2) 



where A/" is a normalization factor and the subscript a = 1, 2, 3 denotes Carte- 
sian isospin orientations. Here, (k) is the creation operator of a pion with 
momentum k and isospin component a and Ja(k) is related to the on-shell 
4-dimensional Fourier transform of the source component Ja{x) through: 



Ja(k) - 



I 



(fxe'^^'Jaix 



(3) 



with ujk 



nit 



Equivalently, it can be directly related to the spatial 



Fourier components of the asymptotic out-going field configuration 7r°"*(k, t) 
and its time derivative 7ra"*(k, t):^ 



Ja(k) = 



i7rr(k,i)+u;,7rr(k,i) 



(4) 



For a given realization of the source, particles are produced independently and 
follow a Poisson distribution characterized by the average: 

n„(k) = (J |al(k) a„(k)| J) = | J„(k)|=^ . (5) 



When the latter is large enough ~ which is required for the present classical 
description to make sense at all, one can neglect the quantum fiuctuations, of 

^ It is easy to checlc that the RHS of Eq. (4) does not depend on time by using the 
fact that 7r°"*(a;) satisfies Eq. (1) with vanishing sources. 
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relative order ~ 1/ ^Jn. The number of pions produced per unit phase space in 
the coUision event corresponding to the source J is then approximately given 
by: 



The magnitude and chiral orientation of the source Ja(k) for each mode k 

fluctuate from event to event. The above picture of a large region of space 
where the chiral condensate is coherently misaligned ideally corresponds to all 
relevant (soft) modes pointing in a given direction e in isospin space: 

/Dcc(k) = J(k) e . (7) 



Notice from Eq. (4) that this implies that both the field and its time-derivative 
be aligned in the same direction e. This can be viewed as an out-going wave 
linearly polarized in isospin space. Clearly, the source (7) induces non-trivial 
correlations between emitted pions with different isospin orientations. This 
can be nicely illustrated by means of the neutral fraction of emitted pions: 



where Na is the total number of pions with isospin component a in a given 
event. For the DCC state (7) one has: 

/Dcc = cos2^^ (9) 



where 9 is the angle between the unit vector e and the third axis in isospin 
space. Assuming that there is no privileged isospin direction, all possible ori- 
entations of the unit vector e are equally probable and one finds that the 
event-by-event distribution of the neutral fraction / is given by [1,2]: 



dPjf) 
df 



DCC 



2^7 



(10) 



The latter exhibits striking fluctuations around the mean value / = 1/3, which 
are a direct consequence of the coherence of the DCC state. Equation (10) is to 
be contrasted with the narrow Gaussian distribution predicted by statistical 
arguments for incoherent pion production. ^ For instance, the probability that 

^ For a binomial distribution, the typical fluctuations around the average value 
/ = 1/3 are ~ l/V-^totj where Ntot is the total multiplicity. 
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less than 10 % of the DCC pions be neiitral is predicted to be as large as 30 %. 
This is what makes the DCC an interesting candidate to explain the Centauro 
events in cosmic ray showers. The property (10) is also at the basis of most 
existing strategies for experimental searches. 

2.1.2 A dynamical perspective 

The idea that multiple-pion emission in high-energy hadronic and/or nuclear 
collisions might be associated with classical radiation is rather old (see e.g. 
Refs. [36,37,38,39]) and can actually be traced back to some old papers by 
Heisenberg [40]. This idea has, however, received only marginal attention un- 
til the early 1990 's, where it has been rediscovered and further developed 
in the modern context of low-energy effective theories [4,5,2,28]. In Ref. [2], 
Blaizot and Krzywicki have investigated the question of soft pion emission 
in high energy nuclear collisions by studying classical solutions of the non- 
linear a model, which describes the dynamics of low-energy pion fields. They 
adopted Heisenberg's ideal boundary conditions [40] (see also Ref. [41]) to 
model the expanding geometry of the collision. The picture is very close to 
the one described above, with the non-trivial pion field configuration gen- 
erated by classical sources localized on the light-cone, that is receding from 
each other at essentially the speed of light. Remarkably enough, there exist 
solutions which correspond to the DCC configuration (7), that is where the 
pion field oscillates in a given direction in isospin space. The analysis can be 
extended to the linear a model [42] and provides an instructive dynamical 
realization of the qualitative ideas described previously. 

Written in terms of the quadruplet of scalar fields ipa = (tt, a), the classical 
action of the linear u-model with the standard chiral-symmetry-breaking term 
reads: 

j d^X ^d^ipa d^ipa - ^ {'Pa'Pa " t^^)^ + H^lfia^ , (H) 

where Ha = (0, H) points in the a-direction in chiral space. The parameters 
V, A and H are related to physical quantities through: 

ml^X{f:-v') (12) 

Note that in the phenomenologically relevant limit, namely rricr ^ or, 
equivalently, A ^ 1, one has v ~ and rricr ~ V2XfT^. In the following, we 
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set the scale v = 1 for simplicity. The classical equations of motion read: 

U + X{^''-v%^a = H,- (13) 



In the original treatment of [2,42], the longitudinal expansion is modeled by 
viewing the colliding nuclei in the center of mass frame as two infinitesimally 
thin (Lorentz contracted) pancakes of infinite transverse extent (see also [41]). 
The symmetry of the problem then implies that the classical field ipa is a func- 
tion of the proper time r = \Jt'^ — only. The field equations (13) become 
ordinary differential equations and can be solved analytically. It is not diffi- 
cult to extend Blaizot and Krzywicki's treatment to the case of a symmetric 
d-dimensional expansion. Here, we present the main results of such an analysis 
and essentially follow the presentation of Ref. [30]. The relevant proper-time 
variable is r = \/t'^ — with — Yli=ixj, and the four-dimensional Lapla- 
cian becomes: 

n<pa{r) = r-''dr{r%ipa) = iPa + -<Pa , (14) 



where the dot denotes derivative with respect to r. Hence, the equations of 
motion involve a friction term ~ (f, which simply reflects the decrease of the 
energy density due to expansion. The initial conditions are to be specified on 
the surface t — tq. 

Prom Eqs. (13) and (14) one easily obtains that: 

— * 

nx^^^ (15) 



and 



TTfT — (TTT 



b H 



(16) 



The first of these equations is a consequence of the conservation of the iso- 

— * 

vector current — n x d^'-n and the second one reflects the partial con- 
servation of the corresponding axial-vector current A^^ — ifd'^a — ad'^if. The 

''' Notice that strictly speaking, this assumes that the expansion never stops, or, in 
terms of the previous baked-Alaska description, that the sources of the pion field 
never decouple. This reflects the fact that the initial energy density is infinite in the 
present idealization. We stress however that, due to expansion, the energy density 
decreases inside the Hght-cone, and the pion dynamics eventually freezes out at a 
time Tf (see Eq. (23) below). Therefore, the present boost-invariant idealization 
should provide a reasonable description if the decoupling (hadronization) time th > 
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orthogonal iso-vectors a and b are integration constants. Their lengths mea- 
sure the initial strength of the respective currents. We shall focus here on the 
regime where a <ti b, which is the relevant one for our present purpose. From 
(15), one sees that the motion is planar in isospin space: Tia = n ■ a = 0. At 
short enough time, the pion mass is irrelevant and one can neglect the second 
term on the RHS of Eq. (16). Following [42], one can show that, for times 
r < {b/m^f": 



TTfoRi— rsin^ (17) 

TTc ~ r cos 9 (18) 

b - 

rcosO, (19) 



where c = axb. The motion is approximately planar in the 4-dimensional chiral 
space. Notice that for a <^ b, the component tt^ is very small and the pion 
field oscillates along the (random) direction defined by the iso-vector b. This 
precisely corresponds to the linearly polarized DCC configuration described 
in the previous subsection. 

For times {b/m^Y^'^ ^ ^ ^ (b/m^y^''', the length of the chiral field r = 
Vo^ + TT^ undergoes rapid, damped oscillations around the approximately de- 
generate minimum of the potential. For instance, for large enough r, one has: 

r«l + ^^°;'%+'' . (20) 
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where C and 6 are integration constants. In this regime one can replace r 
by its time-averaged value f = 1 for all practical purposes and the motion 
essentially takes place near the minimum of the mexican hat potential. One 
therefore obtains that the angle 9 is approximately given by: 

e^bj"^, (21) 



corresponding to a circular motion in chiral space. In this regime the energy 
density is approximately given by: 



(22) 



and is still high enough so that the system does not feel the explicit symmetry 
breaking due to the pion mass, hence the circular motion. 
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At a time r ~ {b/mT^y^'^, the damping produces a cross-over from the circular 
to an oscillatory motion around the true minimum of the potential. This is also 
the time at which the energy density (22) becomes comparable to the mass of 
a pion, divided by the cube of its Compton wavelength ~ m^. Therefore, the 
non-linear dynamics freezes out and one is left with freely propagating pions. 
Indeed, for times r > one finds: 



Vb 



(23) 



and 



62 



(24) 



which describes free propagation of pions in the expanding geometry. These 
correspond to the DCC decay products. As already noticed, when the relative 
strengths of the initial vector versus axial- vector current is small (a -C 6) , the 
pion field is mainly polarized along the random direction b, which results in 
the event-by-event distribution (10) for the neutral fraction /. 

We see from Eq. (23) that the parameter b controls the amphtude of the 
out-going pion field, that is in turn, of the amount of radiated energy. It is 
interesting to compute the probability P{a, b) that a and b take particular 
values in a simple model characterizing the high energy, chirally symmetric 
initial state produced in the collision. Assuming that the values of the fields 
(tt , a) and their proper-time derivative at r = Tq are Gaussian random numbers 
of zero mean and of variance ai and (72 respectively, one easily obtains [42] : 

d^P{a,b) i^i (\/a2 + 6V«^o) e-'/'^" , , 

= — r\J (25) 

dadb 2n^Kl ^ ^ ' 



where K\{z) is the second modified Bessel function and = (Ti(T2To- The 
second term on the RHS shows the behavior at b ^ a, kq: The probability of a 
significant signal appears to be exponentially suppressed in this simple model. 



2.1.3 Coherent state descriptions 

It is instructive to see how the qualitative argument leading to the prediction 
(10) is modified when one takes the quantum nature of the emission process 
into account. This can be done in the previous coherent state picture, Eqs. 
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(2) and (7),^ where the DCC state is characterized by a given orientation 
e in isospin space and a complex source J. For each given source J, isospin 
symmetry is ensured by averaging over aU possible orientations with equal 
weight. This can be described by the following density matrix: 

p{J,r) = J^\J){J\, (26) 



where | J) denotes the DCC state, Eqs. (2) and (7). and the integral is over all 
possible orientations of the unit iso- vector e characterized by the solid angle 
fl. Alternatively, one can imagine [36] a coherent superposition of states |J) 
with all possible orientations, corresponding to the following zero-isospin pure 
state: ^ 

= j <fnY^{n)\j) , (27) 



with the spherical harmonic YQiVt) = l/\/47r. It is easy to show that the two 
descriptions (26) and (27) are actually equivalent in the limit of large particle 
numbers | J(k)p » 1. In that case, one can write: 

p{j,r)^\j){j\. (28) 



To completely characterize the emission process, one has to consider all pos- 
sible realizations of the complex source J, with appropriate weight V{J^ J*) 
(the latter reflects the detailed dynamics of DCC formation, to be discussed 
later in this review). This can be described by the following density matrix: 

pdcc = / vjvr v{j, r) p{j, r) (29) 



^ Squeezed quantum states have also been considered in the literature, see e.g. 

[43,44,45,46]. 

^ Notice that this state is not normaUzed. One has: {J\J) = nk-^(l'^(k)P)' with 
M{x) = 47re~^ sinhar/x. 

The precise meaning of this relation is that correlation functions - or multiplicity 
distributions considered below - computed with either the mixed state (26) or the 
pure state (27) agree in the (classical) limit of large particle numbers. This can be 
seen by computing the generating functional Z{K,K*) = {J\0{K,K*)\J), where 
0{K,K*) = exp / ^[K(k) • Sk + -^*(k) • 4]> ^om which one can obtain aU 

correlation functions by derivatives with respect to the sources K(k) and K*(k). The 
calculation involves an integration over two directions, f cPflcPil'. By employing a 
saddle point approximation in the limit of large particle numbers, one can show that 
Z{K,K*) c:, J ^{J\OiK,K*)\J) = Tt{pO{K,K*)}, with p given by Eq. (26). 
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where the integration runs over all modes in the DCC state: VJVJ* = 
nkciJ(k)ciJ*(k). 

The multiplicity distribution of emitted pions can be obtained as (we consider 
here a single mode k for simplicity and we omit the explicit k-dependence) : 

p(no,n+,n_) = j VJVrV{J,r) | ^ | K n+, n_| J) (30) 



where no and n± denote the actual multiplicity of neutral and positive/nega- 
tive pions respectively. The corresponding state is given by: 



|no,n+,n_) 



|0) 



(31) 



where the creation operators aj and a\ for neutral and charged pions re- 
spectively are related to the Cartesian-coordinates creation operators a\ as: 
aj = al and a^. = ^{a\ ± ia2)/y/2. Correspondingly, it is convenient to de- 
compose the unit iso- vector e by: 60 = 63 = cos 9 and e± = =F(ei ± ie2)/\/2 = 
^sm9e^'^'^/\/2, where 9 and (p denote respectively the azimuthal and polar 
angles in isospin-space. The probability distribution p{nch, no) for charged and 
neutral pions is easily obtained from (30) and can be written as {N = rich+no)' 

p(«ch, no) = P{N) pjv(nch, no) , (32) 



where P{N) is the total multiplicity distribution and PNinch,no) denotes the 
conditional probability that, among N produced pions, rich be charged and no 
be neutral. In the present model, one gets: 

r 1 7P^ 

P(N) = J VJVJ* P( J, J*) e-l-^l' (33) 



and 



PN{nch,no) = ^^ /^|eo|^"o(|e+r + |e_r)"- 
richlno! J 2V/ 

where we introduced the variable / = cos^ 9 in the second line. The last 
integral in Eq. (34) can be expressed in terms of the Euler gamma function 
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r{z) and one gets: 



r(iv + i)rK + |) 



For large mulitplicities no, rich ^ 1, one finally obtains (see also [39,47]): 

PN{nch,no) —^j^ (36) 



where / = Uq/N is the fraction of neutral pions. Thus one recovers the in- 
verse square- root law obtained previously by a simple geometrical argument, 
neglecting quantum fluctuations. As expected, this is a valid approximation 
when the number of produced particles is large. 

Another possibility for the DCC state has been extensively discussed in the 
literature [36,47,48,49,50,51], namely the zero-isospin state with fixed number 
N of pions (note that a state of total isospin zero can only contain an even 
number of pions) [36]: 

|7 = 0, = 2m) = , ^ Ualy - 2aiaLl™ |0) . (37) 
'(2m + 1)! ^ ^ 



To make link with the present considerations, we note that the isospin-avera- 
ged coherent state (27) can in fact be written as a coherent superposition of 
the states (37) with arbitrary (even) number of pions. Indeed, one easily gets, 
after some calculations: 

j2m 

|J) = v^e-^l-^l' ^ — |7 = 0,A^ = 2m). (38) 



m>0 



^(2m + l)! 



It is clearly seen on this expression that production of a DCC is consistent 
with usual conservation laws, such as charge (n+ = n_), isospin, or parity 
(even number of pions) conservation. 

It is interesting to characterize the correlations between the different charge 
states [49,50,51]. A convenient way is to introduce the following normalized 
correlation function — +,0, — ) [52,51]: 

{N{N-l)){n,){n,) ^''^ 
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where brackets denote an average with respect to the multiphcity distribution 
(30). Here, the factor {N)'^/{N{N — 1)) is introduced to correct for possible 
non-Poissonian multiphcity fluctuations. For the model described above, one 
obviously has: 

{ni) = l{N), (40) 



The correlators {niUj) are most conveniently calculated by means of the for- 
mula: {niTij) — Tr{pDCC ofoi where Pdcc is given by Eqs. (29) and (26). 
One finds, neglecting terms of relative order 1/{N): 

(n,n,)^r,,(iV(iV-l)), (41) 



where 

= / ^ |e.P|e,P. (42) 



Finally, one obtains the following correlations [49,51]: 



Roo-T (43) 
5 

Ro+^-l (44) 

R+^^l- (45) 

This is to be contrasted with the corresponding result Rij oc 6ij/{N), ex- 
pected from statistical arguments for incoherent pion emission. Clearly, the 
large correlations obtained here reflect the large event-by-event fluctuations 
of the neutral fraction of emitted pions, see (36). Notice, in particular, the 
negative neutral-charged correlation (44), characteristic of DCC emission. 



2.2 Dynamics of DCC formation: the out- of- equilibrium chiral phase tran- 
sition in heavy-ion collisions 



We have seen that the DCC field configuration is a solution of the low energy 
dynamics of strong interactions. The approaches described in the previous 
sections assume the existence of a classical pion field - or, equivalently, of a 
classical source for the latter. To go beyond this level of description, one needs 
to understand the dynamical origin of such a classical field. An important 
progress has been made in this respect in Ref. [10], where Rajagopal and 
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Wilczck have realized that a strong long-wavelength pion field configTiration 
could indeed be formed during the out- of- equilibrium chiral phase transition in 
the context of heavy-ion collisions: The rapid expansion of the system results in 
a sudden suppression of initial fluctuations (quenching) which in turn triggers 
a dramatic amplification of soft pion modes [10]. 



2.2.1 The quench scenario: amplification of long wavelength modes 

To model the dynamics of the far-from-equilibrium chiral phase transition, 
Rajagopal and Wilczek [10] consider the 0(4) linear sigma model for the chi- 
ral quadruplet = (tt, cr) with action given in Eq. (11) above, where the 
parameters (12) are chosen so that = 135 MeV, ma- ~ 600 MeV and 
/t^ — 92.5 MeV. Anticipating the fact that the relevant field configurations 
correspond to large field amplitudes (that is large occupation numbers), it is 
justified to employ the classical statistical field approximation. In practice, the 
possible initial field configurations at time t = 0, characterized by the values of 
the fields ^q{x) = ^a{^: t = 0) and their time-derivatives 0o(x) = 0a(x, t = 0) 
at each point of space, are sampled from a given statistical ensemble refiecting 
the initial density matrix of the corresponding quantum system. Each such 
field configuration is then evolved in time according to the classical equations 
of motion (13) - supplemented by appropriate boundary conditions. The time- 
evolution of a given physical observable O = 0[(p,i^] is obtained by averaging 
over all possible field configurations. This is summarized in the following for- 
mula, where we omit the explicit spatial dependence as well as chiral indices 
for simplicity: 

_ / dcpodipo V[cpo, ipo] 0[cp{t), cpjt)] 



where / dipodipo represents an integral over possible initial field configurations 
{(fio, (po} with appropriate weight Vl^po- V^o] and {v?(t), <p(t)} is the correspond- 
ing time-evolved configuration at time t. In the strong field regime, where 
quantum fiuctuations are suppressed compared to statistical fiuctuations, the 
previous procedure provides a good approximation to the full (quantum) dy- 
namics (see e.g. [53,54]). 

The physical picture of the Rajagopal and Wilczek scenario is as follows: One 
assumes that a large amount of energy has been deposited in the collision 
zone, corresponding to a very high temperature, presumably well above the 
chiral transition temperature. Due to the rapid expansion, this hot system 
experiences rapid cooling, which results in a strong suppression of the initial 
(thermal) fluctuations. To model this effect, Rajagopal and Wilczek assume 
an instantaneous quench from above to below the critical temperature at the 
initial time. The initial fleld conflgurations are sampled from a chirally sym- 
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Fig. 1. The squared amplitude of the pion field in Fourier space |7r3(k, t)p as 
a function of time (in lattice units) for a given initial field configuration in the 
quenched ensemble (see text), as in Ref. [10]. Low-momentum modes experience 
dramatic amplification, the softer the mode the larger the amplification, and exhibit 
coherent oscillations with period ~ 2-1: /m-n. 

metric probability distribution, characteristic of the high-temperature phase, 

but where the fluctuations are frozen by hand. In the simplest realization of 
this scenario [10], the values of the field V'o(x) and its time-derivative v^o(x) 
at the initial time are chosen as independent random variables on each site 
of a cubic lattice. They are sampled from a Gaussian distribution with 
respective means ((/Jq) = and (y?g) = 0, variances {^q^q) — erf Sab and 
(VoVo) = '^2 ^ab, and covariance {^Pq'^q) = 0, where af and are chosen to 
be smaller than what they would be in the high-temperature phase, thereby 
describing quenched fluctuations. The authors of Ref. [10] choose af = 
and (7| = v'^/AXl, where Ac = a is the correlation length in the initial high- 
temperature phase, here identified with the lattice spacing a. 

The main result of Ref. [10] is reproduced in Fig. 1, which shows the squared 
amplitude of different field modes as a function of time, starting from a given 
field configuration in the ensemble described above: Due to the far-from- 
equilibrium initial condition, one observes a dramatic amplification of low 
momentum modes of the pion field at intermediate times. This phenomenon 
is analogous to that of domain formation after quenching a ferromagnet be- 
low the critical temperature. It is also operative e.g. in the physics of particle 
creation in the early universe in so-called new inflationary scenarios (see e.g. 
Ref. [55]). 



The lattice spacing has therefore the physical meaning of the correlation length 
Ac in the high-temperature phase. 
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It is worth emphasizing that the dramatic growth of field amphtudes for low- 
momentum modes at intermediate times is independent of the particular initial 
configuration chosen here. It is a generic feature of the typical field configu- 
rations in the statistical ensemble described previously. To illustrate this, we 
define an amplification factor as: 

(47) 



where 

3 



P{k,t)=UkJ2MKt) (48) 



a=l 



is the pion power spectrum in mode k at time t, with the average occupation 
number in mode k defined by {cuk — ^Jl^^^rrrt^) : 



n„(k, t) 



i(fia{k,t) +UJk (fa{K t) 



(49) 



in analogy with Eqs. (4)- (5).^^ Here, (/9a(k, t) and (^^(k, t) are the spatial 
Fourier transforms of the field (/^^(x, t) and its time-derivative (^^(x, t) respec- 
tively. Figure 2 shows the amplification factor averaged over a large number 
of field configuration as a function of momentum at two different times [27] 
(see below). One observes that, starting from a quenched initial ensemble, 
the typical field configurations indeed experience dramatic amplification at 
low-momentum, the softer the mode the stronger the amplification. 

This can be qualitatively understood using the following mean- field type ap- 
proximation [10]: In the classical equation of motion (13), replace the quadratic 
term in the square brackets on the the LHS by its ensemble average value: 
(p^ {^P^). Assuming a spatially homogeneous ensemble for simplicity, the 
latter average is a function of t only and the equation of motion for the pion 
field can be written, in Fourier space: 

7r(k,t) = 0, (50) 



This is inspired by the fact that each classical field configuration in the statisti- 
cal ensemble might be viewed as a particular coherent state of the corresponding 
quantum system. 

The approximation described here corresponds to the so-called large- approxi- 
mation. 
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Fig. 2. Momentum dependence of the amplification factor (47) averaged over the 
initial ensemble in the quench scenario at t/ = 10a (21.0 x 10^ events) and tf = 56a 
(10.9 X 10^ events), where a is the lattice spacing. Error bars are statistical. The 
dashed lines are just guide for the eyes. Prom Ref. [27]. 
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Fig. 3. Effective mass squared m^^ (t) as a function of time (in units of the lattice 
spacing a) in the quench scenario [10]. 

where the effective time-dependent mass 



m, 



{t) = X {^'m-v' 



(51) 



The problem is analogous to that of an ensemble of oscillators with time depen- 
dent frequencies. The mass (51) can be viewed as the instantaneous curvature 
of the effective potential seen by the field modes at time t. One immediately 
sees that whenever the effective mass squared (51) becomes negative, soft 
modes with k'^ < |mgg(t)| undergo exponential growth. This is analogous to 
the so-called spinodal instability. This happens when the field fiuctuations get 
small enough: < v"^, which is the case in particular at early times in the 
quench scenario. Figure 3 shows the actual time evolution of the effective mass 
squared defined in Eq. (51), as obtained from the exact classical Monte-Carlo 
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simulation described previously [10]. We see that the amphfication of low- 
momentum modes at early times (t < 10a) is indeed due to the phenomenon 
of spinodal instability. The corresponding amplification factor is shown on the 
left panel of Fig. 2. 

It is clear from Fig. 3 that the spinodal instability alone cannot explain the 
amplification observed on Fig. 1 after times t ~ 10a, where the effective mass 
squared (51) is always positive. In fact, this intermediate-time amplification 
for 10a ^ t < 60a, can be understood as resulting from the quasi-periodic 
oscillations of mlff{t) around its asymptotic late time value [56,57]. This phe- 
nomenon is known as parametric resonance and plays an important role e.g. 
in describing the (p)reheating of the early universe in chaotic inflationary sce- 
narios [58,59,55]. Figure 2 shows the momentum dependence of the average 
amplification factor (47) at two different times, characterizing the two amplifi- 
cation mechanisms at work. We see that most of the observed amplification is 
due to the spinodal instability at early-times, leading to an amplification ~ 80 
for the most amplified mode, whereas the parametric resonance phenomenon 
at intermediate-times only gives an extra factor ~ 5 amplification. Moreover, 
it is important to emphasize that the regular oscillations of the effective mass 
are in fact suppressed when expansion is explicitly taken into account (see e.g. 
Fig. 5 below) and, therefore, so is the intermediate-time parametric amplifi- 
cation. 

Thus we see that, due to the very unstable initial state, which was assumed 
to be formed as a consequence of rapid expansion, typical pion field configu- 
rations undergo a dramatic amplification and a strong long-wavelength pion 
field develops rapidly. Expansion causes the energy density to drop and the 
dynamics eventually linearizes as the modes stop interacting. If this freeze- 
out happens short enough after the collision, the system may be left in such 
a strong field configuration and subsequently decay through coherent pion 
emission. This provides a microscopic scenario for the formation of a strong 
(semi-classical) pion field configuration in heavy-ion collisions. It is worth 
mentioning that this phenomenon has been demonstrated in other models as 



This justifies a posteriori the use of classical statistical field theory. 

One may wonder whether a similar phenomenon could lead to enhanced produc- 
tion of strange mesons due to the approximate restoration of the SUji{3) x 5^7^(3) 
symmetry at very high temperatures. This has been investigated in Rcf. [60] (see 
also [61] for a phenomenological study), with somewhat negative conclusions (which 
can be partly understood as being due to the too large mass of the strange quark). 

It is of common use in the literature to interpret each representative of the statis- 
tical ensemble described here as a possible realization of the pion field formed in a 
given collision event. It is worth emphasizing, however, that physical results, which 
require averaging over the statistical ensemble, do not depend on this interpretation. 
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well [62,63,64]. 
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2.2.2 Expansion and initial conditions 

A crucial assumption in the previous treatment concerns the role of expansion. 
Including expansion explicitly in the description is needed e.g. in order to 
assess whether the decoupling of modes (freeze-out) occurs soon enough before 
rescatterings destroy the large field amplitude and thermalize the system (see 
e.g. Ref. [67]). But the main aspect of implementing expansion is to relax 
the drastic quench assumption made in Ref. [10]. There are two points to this 
assumption: First, the time-scale of the initial energy drop must be sufficiently 
short compared to the typical time characterizing the interactions between 
long-wavelength modes; Second, the amplitude of the initial energy drop must 
be large enough so that, starting in the high- temperature chirally symmetric 
phase, the system is indeed cooled down to the unstable region. Clearly, the 
possible occurrence of an instability will depend both on the efficiency of the 
expansion in quenching initial fluctuations and on the actual initial state of 
the system. 

These aspects have been investigated in details by Randrup in Ref. [68] in the 
context of classical statistical field theory. We review this work here. The effect 
of expansion is modeled by adding a cooling term by hand in the equations of 
motion (13). Although not a rigorous treatment of an expanding system, this 
captures the physics of the cooling process in a simple way. ^® For a symmetric 
d-dimensional expansion, one replaces the 4-dimensional Laplacian by: 



The assumption of quenched initial fluctuations is relaxed and the initial fleld 
conflgurations are sampled from a thermal initial state at a given tempera- 
ture T. Employing a Hartree approximation at t = 0, the initial ensemble 
can be described by a thermal bath of non-interacting quasi-particles with 
self-consistently determined masses and the corresponding probability distri- 
bution V[ipo,^o], see Eq. (46), is simply given by a Gaussian [72]. Assuming 
a high-temperature initial state, the latter distribution is characterized by al- 
most vanishing average values for the chiral fleld and its time derivative and 

-"^^ We mention that a slightly different scenario for the growth of long- wavelength 
pion modes has been considered in Ref. [65]. In this so-called annealing scenario, 
soft pion modes evolve in an effective thermal potential with slowly decreasing 
temperature. This essentially differs from the quench scenario in that the system 
evolves close to equilibrium, which leads to less efficient amplification [66]. 
^® For descriptions in actual expanding geometries, see e.g. [67,69,70,71] (see also 
subsection 2.2.3.2 below). 



D = d^ -A 




(52) 
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Fig. 4. Dynamical trajectories in the plane (To-(5</'^)^/^, where ao = {<j)a4>aY^'^ 
for various expansion scenarios (left) and various initial temperatures (right). The 
dashed line connects the equilibrium states from T = to above 500 MeV and the 
shaded area corresponds to the instability region where the effective mass (51) is 
negative. Successive points on a given trajectory are separated by At = 0.2 fm/c. 
Prom Ref. [68]. 

by thermal two-point functions. Initial field configurations sampled in this 

Gaussian ensemble are then evolved in time according to the exact classical 
equations of motion with cooling, Eqs. (13) and (52). 

The result of the time evolution is most nicely illustrated by plotting the time 
evolution of the order parameter 0^(t) = (pa(t)(pa(t), where (j)a{t) = (v?a(x, t)) 
is the ensemble averaged value of the field at time t, versus that of the 
square fluctuations {Sip'^){t) = ((^^(x, t)) — 0^(t), where = (fia'^a- The corre- 
sponding trajectories are shown in Fig. 4 for various initial temperatures and 
various types of expansions, with d = D — 1,2,3. Due to the cooling term 
in the equations of motion, the initial thermal fluctuations are suppressed - 
the higher the value of d, the more efficient the cooling in (52) - and might 
enter the instability region, where the effective mass squared (51) is negative, 
represented by the shaded area in Fig. 4. When this happens, long- wavelength 
modes get amplified. In particular, the order parameter - which is nothing 
but the zero-mode - grows rapidly and the system is pushed away from the 
instability region. As expected, we see that the phenomenon is more pro- 
nounced for higher values of d, corresponding to more rapid expansions. Also, 
it appears clearly that the time spent in the instability region, which even- 
tually determines the amplitude of the amplification for low-momentum pion 
modes, strongly depends on the initial temperature: For T < 200 MeV, the 
initial order parameter is too large and the system never enters the instabil- 

Assuming a spatially homogeneous ensemble, the latter does not depend on x. 
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ity region. Similarly, one would expect that for very high initial temperatures 
T ^ 500 MeV, the initial fluctuations are too large for the system to ever 
become unstable But one observes that there is a large range of temperatures, 
which roughly correspond to those expected in high-energy nucleus-nucleus 
collisions, where the expansion produces an efficient quenching of initial fluc- 
tuations. 

However, a more detailed analysis of the results shown in Fig. 4 reveals that, 
even for the most favorable case d — 3, the time spent in the instability region 
is rather short and barely exceeds 1 — 1.5 fm/c. This indicates that the typical 
trajectories do not get dramatically amplified, in contrast to what happens 
when the initial fluctuations are artiflcially quenched deep in the instability 
region as discussed before. In fact, the ensemble average amplification factor 
for the most amplified mode k — was estimated in Ref. [68] to be roughly 
of order 10 in the most favorable case. This is very different from the large 
amplifications, typically ~ 100, obtained in the quench scenario. Thus, we 
can already conclude that, in a more realistic scenario with cooling, large 
amplifications seem to be a rather rare phenomenon. This conclusion has been 
corroborated by other studies (see e.g. [73]). A crucial question is, therefore, to 
know how likely these events are. We shall return to this later in this section 
(see subsection 2.2.4). We first discuss the issue of describing the quantum 
corrections to the dynamics in the next subsection. 



2.2.3 Quantum corrections to the dynamics 

A considerable amount of work has been devoted to include quantum fluctu- 
ations in the dynamical description of DCC formation. Initial- value problems 
in quantum field theory are, however, known to be of notorious difficulty. 
Existing studies in the context of DCC physics have been mainly limited to 
the use of mean-field approximations, such as the so-called Hartree approxima- 
tion, or the large- approximation. We describe these developments here, for 
the case of a spherically expanding geometry. We comment on recent progress 
in non-equilibrium quantum field theory beyond mean-field approximations at 
the end of this subsection. 



2.2.3.1 Mean-field approximation In the present context, mean-field 
approximations consist in replacing the fully non-linear problem by an effec- 
tive dynamics in a self-adjusting quadratic potential: The system of interact- 
ing quantum fields is replaced by an ensemble of independent, self-consistent ly 

'^^ Notice, however, that the use of the sigma model cannot be justified at such high 
temperatures. 

For reviews concerning recent developments in this field, see e.g. [11]. 
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dressed quasi-particles. There are various ways one can formulate such approx- 
imations, such as Gaussian ansatze for the density matrix in the Schrodinger 
picture [74], Gaussian-hke factorization of correlation functions in the Heisen- 
berg equations of motion [74,68], non-perturbative resummation of so-called 
daisy and super-daisy diagrams [75], l/A'"-expansion at leading order, where 
= 4 is the number of component of the chiral field [76,77], etc. Here, we 
employ the latter approach: the so-called large-A^ approximation. We shall 
not enter, however, in the details of the 1/A^-expansion and instead, we give 
a simple derivation of the equations relevant for our purposes. 

The exact Heisenberg equations for the chiral quantum field read (see (11)): 

n + \{^\x) - v^)\^a{x) = Ha. (53) 

As we did previously for the classical field equations (see Eqs. (50)-(51)), we 
replace the non-linear term in parenthesis by its average ip^ {'-p'^)) where, 
here, the brackets denote an average over the initial density matrix p, which 
encodes both statistical and quantum fiuctuations in the initial state: (•■■) = 
Tr{p ■ ■ - j/Trlp}. Writing the field as a sum of its average value 0^ = {ipa) 
and a fiuctuation field 5<^„, one gets: 

□ + A(0^(x) + {5ip\x)) - v'')] (M^) + Sip,{x)) = Ha . (54) 

Taking the average of this equation with respect to the initial density matrix 
p, one obtain the following non- linear equation for the condensate (pa'- 

D + M\x)](Pa{x) = Ha, (55) 



where 



(P\x) + {S^\x)) - . (56) 



We see that 4>a{x) evolves in what is essentially the classical potential l^(<^^) = 
A(0^ — v'^Y/A — Ha, modified by quantum and statistical fluctuations through 
a term 5V oc 0^. Finally, subtracting Eqs. (54) and (55), one obtains the 
following Klein- Gordon-hke equation for the fluctuation fleld: 

□ + M2(x)1 (5(^„(x) = 0. (57) 



The excitations of the fluctuation fleld 5ipa describe self-consistently dressed 
quasi-particles evolving in a x-dependent quadratic potential, whose curvature 
depends on the local value of the condensate and on the mean effect of the 
fluctuations themselves. 
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Of course, {5(p^{x)) diverges and the above equations require regularization, 
e.g. through an ultra-violet cut-off A, and renormalization. The quadratic 
divergence in Eq. (56) can be argued to be ^-independent by general power- 
counting arguments and is easily eliminated through the following vacuum 
subtraction [78]: 



M'^ix) - ml 



A 



(^)-/.' + (V)(^)-(o|V|o) 



(58) 



The remaining logarithmic divergence can be eliminated by introducing the 
renormalized coupling constant [76,78]: 

^ _ 1 r k^dk 

Xr~ Stt^J (A;2 + m2)i ^ ' 



It is important to notice that the cancellation of the logarithmic divergence 
at any time requires a suitable choice of initial conditions (see Eqs. (80)-(81) 
below): Not all initial conditions are physically acceptable. Finally, it should 
be kept in mind that the present model is at best an effective theory for low- 
momentum scales. The ultra-violet cut-off has therefore a physical meaning 
and should be kept finite. The aim of the renormalization procedure is to 
reduce the cut-off dependence of physical results. 



2.2.3.2 Expanding geometry We describe the mean-field equations in 
more details in the case of a spherically expanding system, for which we es- 
sentially follow the treatment of Ref. [70]. This requires one to quantify the 
theory on equal-proper-time hyper-surfaces r = const, see Eq. (67) below, 
where the proper-time r = \/t'^ — r^, with r the radial distance to the origin, 
is the time measured by a co-moving observer sitting on an expanding shell. 
Exploiting the spherical symmetry of the problem, the relevant coordinates 
are: 

r = , ^ = ^In (^) J,^, (60) 



where r] is the spatial radial rapidity, and 9 and the usual azimuthal and 
polar angles on the sphere. With this choice of coordinates, the relativistic 



22 Moreover, it is well- known that the theory becomes trivial as A — > oo [79]. The 
physical cut-off should be chosen smaller than the Landau pole A^, which, in the 
present approximation, is roughly given by: h-L/m-j^ « exp(87r2/A''Ai^). 
2^ For details concerning techniques of quantization in curved geometries, see [80]. 
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metric is of the Priedrich-Robertson- Walker type and reads, explicitly: 

ds^ = g^^ dx^dx" = dr^ - [hij dx'dx^) , (61) 



where hij is the metric on the 3-hypcrboloid, or pseudo-sphere, r = 1: hij — 
diag(l, sinh^ ?7, sinh^ ?7 sin^ The four- dimensional Laplacian reads: 

□ = A(^) = + -d, - 1a(3) , (62) 

T T 



where A*^^) is the corresponding three-dimensional Laplacian on the unit pseu- 
do-sphere: 

A^^^ ^ ^diVhh'^ di, (63) 
Vh 



with Vh — sinh^ rj sin 6 the determinant and k^^ the inverse of the three- 
dimensional metric hij. 

Spherical symmetry and radial-boost invariance ensures that both the con- 
densate (pa and the effective mass only depend on the proper-time r and 
Eq. (55) reads: 

Mr) + - Mt) + M\t) Mr) = Ha , (64) 



The canonical quantization of the fluctuation field is more conveniently for- 
mulated in terms of the dimensionless rescaled fields (the various components 
of the fluctuation field dip^ being all equivalent in the present approximation, 
Eq. (57), we drop the chiral indices for simplicity): 



(^(r, x)=r599(r, x) 

7r(T, x) = t(^(t, x) = (p'{t, x) 



(65) 
(66) 



where the dot represents the derivation with respect to r and the prime rep- 
resents the derivation with respect to the conformal-time variable u — Iut/tq. 
These operators satisfy the equal-time commutation relations: 

[(P{t, x), 7r(T, x')] = i 5x,x' , (67) 

The rescaled fluctuation field introduced here should not be confused with the 
original field ipa- No confusion being possible in the present section, we use the same 
notation to avoid a proliferation of symbols. 
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where we introduced the short-hand notations x = (77, 6, (f) and (5x,x' = S{v ~ 
7]')5{9 — 9')5{(p — (p')/\/h. Exploiting the symmetry of the problem, one projects 
the fluctuation fleld ip on eigenfunctions of the three-dimensional Laplacian 
(63): 

<Ps{r) = / d^x^/hy:{^) <p{r, x) (68) 
7r,(r) = J d'xVhy:{^) 7r(r, x) , (69) 

where the so-called hyperbohc harmonics are defined as [81,82]: 

A(3)3;,(x) = -(s2 + i)3;,(x). (70) 

Here we introduced the notation s = (s, m), where s is a positive real variable 
and I and m take integer values: < Z < -|-oo and —l<m<l. The hyperbolic 
harmonics form an orthonormal basis, spanning the space of functions on the 
unit pseudo-sphere. One has, in particular: 



/ 



d'a:x/^X,(x)3;,(x) = V, (71) 



where 5s,s' = 5{s - s') 5mm'- 

In the mean-field approximation, the equations of motion (57) for the fiuctu- 
ation field are essentially linear and the most general solution can, therefore, 
be written as: 



(^,(t) = V'.(t) a, + ^^^:{T) {-ir (72) 
7rs(r) = i^'^r) a, + rj (r) (-1)- aU , (73) 

where, as before, the prime denotes differentiation with respect to u = Inr/ro. 
The annihilation and creation operators Ug and for quasi-particle excitations 
satisfy the following commutation relations: 

[as,a^,] = , (74) 



all other commutators being zero. The corresponding mode functions V's(t) 
satisfy the following differential equation: 

Mr) + -Mr) + cOs{r)Mr) = , (75) 

T 

Here, we have used the fact that the field 9j(r, x) is real, which imphes, using 
3^*(x) = (-l)™3;_s(x), with -s = {s,l,-m), that (^i(r) = (-l)"V-s(r). 
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which describes a set of damped oscillators with time-dependent frequencies: 



u^sir) ^ J + M^{t) . (76) 



The physical momentum of the corresponding quasi-particle, k = s/r, is red- 
shifted as time increases, as a simple consequence of expansion. In terms of 
the variable u — Iut/to, one has: 

<(r)+(^,2(r)^,(r) = 0, (77) 



where a;s(r) = tu!s{t) is a dimensionless frequency. The problem therefore 
reduces to the familiar one of a set of parametrically excited oscillators. The 
non-linear character of the mean-field dynamics enters through the calculation 
of the self consistent mass, Eq. (56). For instance, assuming that in the initial 
state: 

(aiosO = (5s,s' , (asa,,) = 0, (78) 



one has, using the property J2im |3^s(x)P = s^/27r^ [81]: 

= 4>\t) -v' + -J ^ |^,(r)|2 (2n, + 1) , (79) 





where N — A is the number of components of the chiral field. The integral on 
the RHS is limited to physical momenta k — s/r < A. Finally, the quantization 
is completed by choosing appropriate initial conditions for the mode functions. 
The cancellation of divergences required by the renormalization of the mass 
gap is ensured if one adopts the following adiabatic condition at r = tq: 



V's(to) = 



V^2a;3(ro) 



(80) 



With this choice, the solution of Eq. (77) for high momentum modes, s^/r^ S> 
|M^(r)| is approximately given by the second-order adiabatic function: 



Mr) 



exp < 



u 

l J du' ujs{u') 



where u = In t/tq. Using this solution, it is easy to check that the divergent part 
of the momentum integral in Eq. (79) is exactly canceled at all times by the renor- 
malization procedure described previously, Eqs. (58)-(59). 
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2.2.3.3 Physical particles, interpolating field and amplification fac- 
tor Assume that Heisenberg and Schrodinger representations coincide at 
T = To- With the choice (80)-(81) of initial mode functions, one can regard 
as the Schrodinger representation operator creating a particle with frequency 
ci;s(ro). These are the physical excitations appropriate for the description of the 
initial state. Similarly, one can introduce the particles with frequencies uJs{Tf) 
appropriate to the final state at t — Tf. Let bl denote the corresponding 
creation operator in the Schrodinger representation. The Heisenberg represen- 
tation field operators, given by (72)-(73), can then also be written as: 



Mr) = qshir) + q:i-irbUr) 



(82) 
(83) 



with 



Ps = -t\ 



1 - z^^" 



(84) 

(85) 



and 



b,{r)^U{r,ro)b,U\r,ro), 



(86) 



where U (r, tq) denotes the unitary timc-cvolution operator, which connects the 
Heisenberg and Schrodinger representations. Using (72)-(73) and (82)-(83), 
one easily obtains the Bogolyubov transformation connecting the operators 
as, al and b^ir), bKr): 

6s(r) = a.(r) a, + /5,(r) (-I)'" al, (87) 
(-I)'" bUr) - P:{t) a, + alir) (-I)'" aU , (88) 

where 



ia*(r) = qs i^* (r) - ^^(r) 
iP*s (^) = Qs i^'s (^) - Ps i^s (r) . 



(89) 
(90) 
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Fig. 5. Proper-time evolution of the effective mass squared M^(t) for the initial 
condition of Ref. [70], namely: (/^^(to) = (0.3,0) [fm"!] and <j)a{To) = (-1,0) [fm^^], 
at To = 1 fm. The initial temperature is T = 200 MeV. The employed parameters 
are: A = 800 MeV, = 139.5 MeV, U = 92.5 MeV and = 7.3. The dashed 
line shows the value of the pion mass squared = 0.5 fm~^. 
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Fig. 6. The amplification factor As{Tf) corresponding to the event of Fig. 5, for 
various values of the final time r/. 

Assuming that the initial state is characterized by the correlations (78), one 
finally obtains the number of 6-particles at the final time t/ for a given number 
of a-particles at the initial time Tq as: 

nf\rf) + \ = MTf)(nf^ + l) , (91) 

The derivation given here differs from that of Ref. [70] , where the so-called adia- 
batic basis is unnecessarily used [26]. 
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where nf^r) = (6J(t)6s('?')). The amphfication factor is given by: 



(92) 



Figure 5 shows the time evolution of the effective mass squared M^(r) for a 
particular set of initial conditions, taken from Ref. [70]. Here, one assumes 
a local thermal equilibrium at temperature T in the initial state and sets, 
accordingly: 

= , ^ . (93) 



One sees that, for the considered initial conditions (see the caption of Fig. 5), 
the system is driven into the unstable region and the effective mass squared 
becomes negative for some period of time. This leads to a dramatic amplifi- 
cation of low-momentum modes, as shown in Fig. 6. Notice that due to the 
freeze-out of interactions at large enough Tf, one has M'^^Tf) ^ and the 
particles in the final state indeed correspond to physical pions. These results 
are the quantum analog of those obtained previously using classical statistical 
field theory techniques and demonstrates that significant amplifications can be 
obtained for some appropriate choice of initial conditions {(f)aiTo),<j>a{To))-^^ 



2.2.3.4 Non-equilibrium dynamics beyond the mean-field approxi- 
mation Mean-field approximations, such as the large- approximation em- 
ployed above, neglect scatterings between quasi-particles. The latter tend to 
redistribute the energy stored in the amplified modes among other modes and 
to thermalize the system, hence washing out any signal in the final state. 
It is a crucial issue for the quench scenario to be successful, that freeze-out 
occurs before the early-time amplification get washed out. 

The description of the far-from-equilibrium dynamics of quantum fields be- 
yond mean-field-like approximations has long been a major difficulty in prac- 
tice: Similar to perturbation theory, standard approximation schemes, such 

In the present set-up, the initial conditions are completely determined by the 

choice of the eight real numbers {(f>aiTo), (f'aiTo))- From these, one computes the 
initial mass squared M^(ro) and its time derivative ((iM^/(ir)(ro), which are then 
used to initialize the mode functions (V's('ro)) '^s('^o)) through Eqs. (80) and (81). 

This effect is visible on Fig. 1 at late times. In contrast to mean-field approxi- 
mations, the classical field approximation employed there treats the full non-linear 
character of the problem exactly. Notice, however, the processes of classical versus 
quantum thermalization are very different and, in particular, they are characterized 
by different time-scales (see e.g. Ref. [54] ) . 



30 



as e.g. a coupling or a 1/iV-cxpansion of the one-particle-irrcduciblc (IPl) ef- 
fective action, are plagued by the problem that a secular (unbounded) time 
evolution prevents the study of the late-time behavior. Important progress in 
this field has been made in recent years with the use of efficient functional 
integration techniques, so-called n-particle-irreducible (nPI) effective actions 
[83,84,85], for which practicable non-perturbativc expansion schemes are avail- 
able. In particular, it has been demonstrated that approximations based on 
a systematic coupling or 1/A^-expansion of the 2P1 effective action allow for 
practical, first principle calculations of the far-from-equilibrium dynamics as 
well as late-time quantum thermalization [86] . ^° 

Clearly, it would be of definite interest to apply these methods to the problem 
of DCC formation. For instance, the large- iV calculation described in the pre- 
vious subsection can be systematically improved using the 2PI l/A'"-expansion 
[88,89]. In particular, this provides a valid description of the dynamics at non- 
perturbatively high occupation numbers [90], as is the case in the context of 
the out-of-equilibrium chiral phase transition. 

2.2.4 The likelihood of DCC formation 

The quench scenario provides a plausible mechanism for DCC formation in 
high-energy heavy-ion collisions. A proper description of the cooling due to 
expansion reveals, however, that significantly large amplifications are not very 
likely. It is therefore crucial for phenomenology to estimate the probability 
that this happens. It is true that the present models are not realistic enough 
to be trusted at the quantitative level. Nevertheless, even within the existing 
framework it is legitimate to seek for a crude estimate of the probability in 
question. This problem has been addressed in Ref. [26] (see also [73]), which 
we review in the present subsection. 



2.2.4.1 Physical picture and sampling strategy Consider the evolu- 
tion of a spherical droplet of DCC in its rest frame. One starts with a small ball 
of radius Rq, filled with hot matter in local thermal equilibrium at a tempera- 
ture T and assumes that the ball undergoes a radial expansion at the speed of 
light. Due to the time dilation the temperature stays approximately constant 
within a layer near the boundary of the ball. The equations of motion of the 
sigma model are supposed to describe what happens in the interior of the ball, 
the cooling observed as one moves away from the surface towards the center. 
As a first approximation, it is simpler to assume that the expanding bubble 
is connected forever to a heat bath kept at constant temperature T, so that 
the process never stops. Of course, the hot shell becomes thinner with increas- 

^° Recent reviews on this topics as well as relevant literature can be found in [11,87]. 
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ing time and eventually disappears. We shall discuss the effect of switching 
the heat bath off at the end of this subsection. The relevant boundary con- 
ditions for the pion field inside the bubble are specified on the hyper-surface 
T — To = Ro, where the proper-time r = y/W^^r^, with r the radial coordi- 
nate. A given field configuration completely determines the amplitude of the 
resulting pionic amplification. The point is that the probability of occurrence 
of a particular initial field configuration can be completely determined once 
one has assumed that the initial droplet is in local thermal equilibrium. 

The quantum evolution of the system is treated in the large- approximation 
discussed previously. In that case, one can convince oneself (see also footnote 
28 above) that the only relevant initial conditions are those concerning the 
average- value of the field (pai^o) and its proper-time derivative 4>a{To)- For the 
sake of the argument, consider first the case of a classical field. The key obser- 
vation is to realize that an observer living within the ball cannot distinguish 
the initial values (pai^o) and <pa{To) from the spatial averages (pa and Tta of the 
field and of its time derivative respectively, calculated over the volume Vq of 
the initial ball. The point is that the latter variables fiuctuate in a predictable 
manner: Assuming that the initial ball, of fixed volume Vq, is in local thermal 
equilibrium at temperature T means that the field fluctuates as if the ball 
were part of a larger system, with volume V much larger than Vq, in equi- 
librium at the same temperature T. In this large system the variances of the 
spatial averages of the field and of its time derivative are very small, of order 
1/V (since long range correlations are absent). The corresponding variances 
for spatial averages over the volume of the ball are just the same, enhanced 
by a factor of order V/Vq (assuming that the radius of the ball is at least of 
the order of the correlation length at temperature T). 

In the quantum theory, the c-numbers 4>a and 7fa correspond to the possible 
measured values of the observables: 



They can be sampled from the probabihty distribution characterized by the 
means: 




(94) 



(95) 



E[7r„] = (n„)^ = 0, 




2 "i*^ ' 



(96) 



(97) 
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and by the variances (no sum over a): 



Var 



{^l)r - {^a)l = ^ J d'xd'y {Sip,{^)SipM)T (98) 



° Vb 



and 

Var [7f„] = (n^), - {Ua)l = ^ I d'xd'y (50„(x)50„(y)), , (99) 



Vo 



where thermal averages {■ ■ correspond to the equihbrium density ma- 
trix p oc exp(— if/T), with H the Hamiltonian. The thermal correlators in 
Eqs. (98) and (99) can be written as: 

{S^^{^)ScpM)r - / (0 GUk) e'^<-y^ , (100) 

and similarly for (50<j(x)5(^a(y))j, with G^{k) K^{k). In the large-A^ ap- 
proximation, o ne haa: G ^j/c) = coth(^fc/2T)/2Efe and Kf^{k) = ElG^{k), 

where Ek — ^Jk"^ -\- M^. Here, the equilibrium mass-gap at temperature T, 
Mt, is the solution of the gap-equation (56) written in thermal equilibrium: 



where 0t — H/M'^ is the thermal condensate (cf. Eq. (96)). 

In the large-A^ approximation the interacting system is effectively replaced 
by an ensemble of independent excitations with mass Mt and the probability 
distribution we are interested in is, therefore, simply given by a Gaussian. The 
latter is completely characterized by the parameters given in (96)- (99). The 
variances (98)- (99) can be estimated analytically when the radius Rq of the 
ball is much larger than the correlation length At = 1/Mt- In that case, it 
is sufficient to calculate the variances for the quasi-infinite volume V and to 
multiply the result by V^/Vq, since the V/Vq small cells fluctuate independently. 
Obviously, only the modes k <l/Ro contribute to the integrals in Eqs. (98)- 
(99). Therefore, for large volumes, there is only the zero mode contribution 
and one immediately gets: 



Var 



0a 



1 



Var [TTa] 



cm 

Km _ 

Vo 2Vo 



2VoMt 
Mt 



coth 



Mt 
2T 



coth 



Mt 
2T 



(102) 
(103) 
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The dispersion of (pa calculated exactly [91] is smaller by a factor of 2 (1.5) 
for Rq = Xt (i?o = 2At) and T = 200 to 400 MeV. The dispersion of 7f„ 
obtained from (99) differs from the estimation (103) by 20% (8%), respectively. 
For i?o < At the discrepancy between the analytical formulas and the exact 
results increases rapidly [91]. Finally, as expected, the fluctuations within the 
ball depend rather weakly on the environment provided Rq > Xt- 

The formalism of Ref. [70], reviewed in subsection 2.2.3.2 above, together with 
the sampling method proposed in Ref. [26], and described above, enable one 
to estimate the likelihood of a coherent amplification of the pion field. More 
precisely, one can calculate the probability that the amplification factor Ag 
given by Eq. (92) takes a given value. In such a calculation the size Rq of 
the initial ball appears as a free parameter. Remember, however, that one 
has to set tq = Rq and that the friction force responsible for the quench is 
proportional to l/r. Thus the hkelihood of DCC formation decreases with 
increasing Rq. This parameter should be assigned the smallest possible value 
in order to get an upper bound for the probability we are looking for. 

At this point, it is important to realize that the model described here only 
makes sense for Rq > Xt- Indeed, the concept of local thermal equihbrium is 
meaningful for a cell whose degrees of freedom fluctuate more or less indepen- 
dently from what happens in the neighbouring cells. Also, the validity of the 
mean-fleld approximation requires the size of the cell to be larger than the 
Compton wavelength of a typical excitation. With these arguments in mind, 
one focuses on the values of Rq in the range of one or two correlation length. 



2.2.4.2 Results The main result of Ref. [26] is reproduced in Fig. 7, which 
shows the histograms of the probability P{A) that the ampliflcation factor of 
the 5 = 0"^ mode takes the value Aq = A. The relevant parameters are the 
same as in Ref. [70], sec Figs. 5 and 6. The ampliflcation is calculated at 
Tf = 10 fm/c, where the system is in the stationary regime. Clearly, large 
amphflcations occurs in a small fraction of events only. It is interesting to 
notice that one condition for a large ampliflcation appears to be the smallness 

— * 

of the absolute strength of the initial classical iso- vector current = tt o x 

— * 

S^TTo as compared to that of the corresponding axial current A^ = Ti^d^aQ — 
(^od^Tio, where (7fo,(To) = 0a(To) and (vfcfjo) = 0a(To). This is iUustrated by 
the histograms of Fig. 8, which shows the the respective currents versus the 
ampliflcation factor A. Recall that a similar condition was required in order 
that the classical solution of Blaizot and Krzywicki, described in Sec. 2.1, 
actually corresponds to a DCC, that is to a linearly polarized in isospin space. 

In order to judge what ampliflcation should be considered as large for phe- 
nomenological purposes, one can estimate the multiplicity of produced pions. 
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Fig. 7. Probability distribution of the amplification factor A of the softest mode for 
various choices of the initial temperature T and size of the ball Rq (in units of the 
correlation length). All histogram are made from at least 10^ events. From Ref. [26]. 
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Fig. 8. The strengths of the initial classical iso-vector and iso-axial-vector currents, 
Vn ■ and • A^^ respectively (in fm^^), versus the amplification factor A. The 
parameters are Rq = Xt and T = 200 MeV. Prom Ref. [26]. 

The invariant one-pion inclusive spectrum can be obtained as [92]: 



(104) 
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where g is the determinant of the metric tensor, see (61), f{x,p) is the in- 
variant phase-space density of produced particles and u^^ = dx^^/dr = x^/t 
is the 4-velocity of the co-moving volume centered at x^. It is a unit 4- vector 
orthogonal to the hyper-surface t = Tf, where the particles are counted. In the 
present case, one has [26]: f{x,p) — N{2T:)~^nf^{Tf), where nf'^Tf) is given by 
Eq. (91) above. The momentum s/r/ in the expanding frame is related to the 
4-momentum p'* in the rest frame by the obvious relation [70] (see Eq. (76)): 



^ + Mf, (105) 



\ 



where M| = M'^{Tf) ^ is the mass gap at time Tf. In the present case, the 
integrand in (104) depends on a single external 4- vector, i.e. p^, and therefore 
the integral, being a Lorentz scalar, only depends on — Mj: The invariant 
spectrum is flat. Using hyperbohc coordinates, see Eq. (60), one finally obtains 
[91]: 




(106) 



Notice that the RHS of Eq. (106) does not depend on the choice of Tf, provided 
the latter is large enough (see also Fig. 6). 

The flatness of the spectrum is, of course, an artifact of the unreahstic assump- 
tion that the boost invariant expansion continues forever. In a real collision 
process the expansion would last a finite time and the resulting spectrum 
would be cut, the value of the cut-off reflecting the behavior of the envi- 
ronment. Obviously, the predicted total multiplicity depends strongly on this 
cut-off and cannot be estimated in a reliable manner within the present model. 
However, the RHS of (106) is presumably a reasonable estimate of the invari- 
ant momentum-space density of soft pion radiation (see, however, [35]). The 
latter can be compared to the corresponding density of incoherently produced 
pions. A simple example is instructive: The one-particle spectrum in the cen- 
tral rapidity region of a heavy-ion collision can be parametrized as follows: 



E— 
d^p 



dn 



dyd^Pt 




(107) 



where h denotes the height of the rapidity plateau. In a central Pb-Pb collision 
at the CERN SPS one observes [93] about 200 vr^ per unit rapidity, i.e. for 
all pions h ~ 600. Thus the invariant momentum space density at very small 
transverse momentum is roughly 1900 GeV~^, using (pf) = 0.1 GeV^. The 
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corresponding density fluctuation is expected to be ~ ~ 75 GeV~^. 

For the DCC signal to be detectable, the RHS of (106) should be significantly 
larger than the expected fluctuation. One finds [26] that the signal is more 
than three times the fluctuation for events where the ampliflcation factor for 
the zero mode is ^4 > 45. Setting Rq — At, the corresponding probabihty is 
roughly 4 x lO'^ both for T = 200 MeV and T = 400 MeV. Of course, this 
is a conditional probability, as it has been assumed that the initial plasma 
droplet was formed in the collision. Moreover, one should keep in mind that 
the present estimate has been obtained in the most favorable scenario and 
should, therefore, be considered as an upper-bound. 

Thus, the probability of a potentially observable DCC signal appears small. 
The present crude estimate indicates that in a central Pb-Pb collision at CERN 
SPS this probability is at best: 

Probability ("observable" DCC) < 10"^ . (108) 

This should be taken into account by experimenters designing DCC hunt 
strategies as well as theorists studying possible DCC signatures. It is worth 
emphasizing that the above prediction agrees with the present experimental 
upper limits for DCC formation obtained at CERN by both the WA98 and 
the NA49 collaborations [18,19,21].=^^ 

2.2.5 The unpolarized DCC 

The quench scenario has been widely accepted as a microscopic description 
of DCC formation in heavy-ion collisions and, as such, has been extensively 
used to study DCC phenomenology [94,95,71], see Sec. 2.4. However, if the 
rapid suppression of initial fluctuations indeed provides a robust mechanism 
for the formation of a strong coherent pion fleld, it is not clear whether it 
can explain the hypothetical polarization in iso-space originally proposed (see 
e.g. Eq. (7)). In fact, the expected large event-by-event fluctuations of the 
neutral ratio, see Eq. (10), have never been observed in actual simulations of 
the out-of-equilibrium phase transition [96,94,95]. Of course, deviations from 
the ideal law (10) are to be expected in a realistic situation, but the complete 
absence of large fluctuations in the quench scenario suggests that the original 
picture of linearly polarized waves coherently aligned in isospin space might 

We stress that the present estimate is pertinent to the CERN SPS experimental 
conditions and would certainly lead to a different result for RHIC or LHC energies. 

Naively, one would expect that the likelihood of a potentially observable DCC fluc- 
tuation decreases at higher collision energies, where both the initial temperature 
and the average multiplicity of normal pions are increased compared to the present 
case. 
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Fig. 9. Neutral fraction distribution in a single mode, here k = 0, at initial time 
(dashed line) and at final time, here = 10 fm (histogram). Prom Ref. [97]. 

actually not be realized in the context of the far-from-equilibrium chiral phase 
transition. 

It is a question of great phenomenological relevance to inquire to what extent 
is the original picture realized in a realistic microscopic model. The point is 
that the usual assumption of a (locally) thcrmalizcd initial state implies that 
the field modes are completely uncorrelated at initial time. In order to actually 
generate a DCC configuration, the microscopic mechanism at work needs not 
only to be efficient in amplifying the modes amplitudes, but it should also build 
correlations between amplified modes. This question has been investigated in 
Ref. [97], where a detailed statistical analysis of the pion field configurations 
produced after a quench has been performed, with particular emphasis on 
their isospin structure. For this purpose, it is sufficient to consider the original 
model of Rajagopal and Wilczek [10], see subsection 2.2.1, which includes all 
the relevant physical features. In particular, the classical field approximation 
allows one to take into account the full non-linear character of the problem, 
which is crucial to study correlations. Moreover, in the quench approximation, 
the system is artificially prepared in an unstable situation, all configurations 
undergo a dramatic amplification, as in Fig. 1. This automatically selects the 
events which are relevant for the present analysis. 

To analyze the isospin orientations of distinct field modes, a sensitive observ- 
able is the neutral fraction of pions in each mode k, defined as: 



m 



n3(k) 



(109) 



ni(k) +n2(k) +n3(k) ' 



For a more detailed analysis, see Ref. [97]. 
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Fig. 10. Normalized correlation between the fluctuations of the neutral fraction in 
the zero mode and in mode k as a function of momentum. Prom Ref. [97]. 
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Fig. 11. Event-by-event distribution of the neutral fraction in a bin in momentum 
space containing 10 amplified modes. The expected DCC signal is already consid- 
erably reduced, due to the absence of correlation between modes. From Ref. [97]. 

where ria(k) = na(k,tf) represent the averaged occupation numbers corre- 
sponding to the classical field configuration at the final time tf, see Eq. (49). 
The event-by-event distribution of the neutral fraction (109) for the most am- 
plified mode, k = 0, is shown in Fig. 9 for = 10 fm/c, corresponding roughly 
to the end of the spinodal instability period, cf. Fig. 2. The corresponding dis- 
tribution at time tf — 56 fm/c, characteristic of the parametric amplification 
mechanism, looks exactly the same [97]. In both cases, all the amplified modes 
exhibit similar distributions. In fact, the neutral fraction distributions for am- 
plified modes are found to be essentially the same as the corresponding ones in 
the initial state. The latter can be computed exactly [97,91] and is represented 
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by the dashed hne in Fig. 9. For a given mode k, it reads: 



Pkif) 



Fa,{f) + F.nM) 



(110) 



where 



n + 1 



3/2 



Faif) /)] 



O - (1 - 2/) 



(111) 



and 




(112) 



with af and o"2 the variances of the field and its time derivative respectively 
(see subsection 2.2.1). Although not exactly l/\/7, the distribution is very 
broad, exhibiting large fluctuations around the mean value / = 1/3, which is 
the relevant point for phenomenology. As emphasized previously, these large 
fluctuations are a direct manifestation of the classical nature of the fleld modes. 
The deviations from the ideal law (10) come from the fact that the latter are 
not strictly linearly polarized waves [97]. 

However, one finds that distinct modes have completely independent polariza- 
tions in isospin space, as can be seen on Fig. 10, which shows the statistical 
correlation between the neutral fractions in different amplified modes: Their 
directions of oscillation in isospin space are completely random. In other words, 
different modes act as independent DCCs. This has the important phenomeno- 
logical consequence that the large event-by-event ffuctuations of the neutral 
fraction are rapidly washed out when the contributions of several modes are 
added in a momentum bin, even when one limits one's attention to soft modes 
only. This is illustrated in Fig. 11, which shows the neutral fraction distribu- 
tion in a bin containing 10 modes. The expected signal is considerably reduced, 
already for such a small bin. This explains the absence of large ffuctuations 
reported in previous studies [96,94,95], where the authors typically considered 
the contribution of a large number of modes. 

In conclusion, the non-linear dynamics of the linear sigma model does not 
build the required correlation between modes: The state produced in the sim- 
plest form of the quench scenario, where no correlations are present initially. 

The following formula is obtained for Neumann boundary conditions (where the 
Fourier components of the field are real numbers) which are convenient for discussing 
the question of polarization. 

Similar conclusions have been reached in a slightly different context in Ref. [98] . 
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Fig. 12. Schematic representation of the ideal DCC configuration (left), where all 
field modes oscillate in the same direction in isospin space (see Eq. (10)), and of the 
unpolarized DCC formed after a quench (right), where the modes have independent 
directions of oscillation. 

is not identical to the originally proposed DCC. Instead, a more realistic pic- 
ture is that of a superposition of waves having independent orientations in 
isospin space: an "unpolarized" DCC configuration, as depicted on Fig. 12. 
Of course, one cannot exclude the possibility that the required correlations 
between modes be indeed (partially) formed in actual nuclear or hadronic 
collisions by means of some other mechanism. Reality most probably lies 
between the two extreme pictures represented on Fig. 12. The point is that 
the only presently existing microscopic scenario for DCC formation predicts 
an unpolarized state. This information should be taken into account in further 
theoretical as well as experimental investigations. 



2.3 Lifetime of the DCC 

So far we have been mainly concerned with the description of the intrinsic 
DCC dynamics, that is the possible mechanism responsible for its formation. 
We now discuss extrinsic aspects of DCC dynamics, namely those concerning 
the interactions of the DCC bubble with its hadronic environment. If at all, the 
DCC is produced in a high-multiplicity environment and must be considered 
as an open system interacting with the surrounding degrees of freedom [100]. 
In the context of ultra-relativistic heavy-ion collisions, one typically imagines 
a thermal bath of pions and nucleons. The interactions with the latter tend to 
destroy the coherence of the DCC excitation and, if it were not for expansion 
- which eventually causes the interactions to freeze out, the produced DCC 
state would simply melt in the external bath: The DCC pions would thermalize 
before being emitted. The possibility to detect any hadronic signal from a DCC 
therefore crucially depends on the lifetime of the latter in the environmental 
bath. 



See, for instance, Ref. [99]. 
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2.3.1 Interactions with the hadronic environment: effective dissipation 

The interactions of the DCC with the debris of the coUisions have been studied 
by many authors in various contexts [100,73,64,101,102]. Although a quanti- 
tative description is clearly out of reach within the present theoretical under- 
standing, a qualitative discussion is already useful. A widely used approach 
in the literature consists in integrating out the degrees of freedom of the bath 
to obtain an effective dynamics for the DCC state [73,103,104,105,106]. As a 
simple picture, one imagines a DCC excitation, characterized by a classical 
pion field tTq, in contact with a thermal bath of pions and nucleons at a given 
temperature T. The DCC being essentially a soft excitation, it is reasonable 
to assume that the dynamics of the bath is characterized by comparatively 
short time scales. Tracing over the rapid degrees of freedom in analogy with 
the standard description of Brownian motion [107,108], one essentially obtains 
an effectively dissipative dynamics for the DCC pion field. In general, the lat- 
ter is characterized by a non-local memory integral and an associated colored 
random noise, reflecting the presence of the bath [109]. For the present ar- 
gument, however, it is sufficient to assume the so-called Markov limit [109], 
where the dissipation kernel becomes local, thus leading to a simple friction 
term ~ rjkTraQs., t) in the equations of motion for the DCC pion field excitation 
in momentum space, and where the associated fluctuating noise field ^^(t) is 
a white noise (see Eq. (114) below). The relevant equations of motion for a 
given realization of the noise field ii^{t) can be written as: 

7f„(k, t) + rik 7r„(k, t) + ulnaiK t) + ■ ■ ■ ^ ^^(t) , (113) 



where ook = yk'^ + and were, the dots on the LHS represent non-linear 
terms in the field, which we neglect for the present discussion. Physical 
results are obtained after averaging over all possible realizations of the noise, 
with an appropriate weight. Denoting by (. . the corresponding average 
over the degrees of freedom of the bath, one has, for instance, (^k(^))B ~ 0- 
Moreover, for a bath in thermal equilibrium the correlation function of the 
random noise and the corresponding damping coefficient rjk are related by: 

mtm')k = *-.ooth (^) '"'"''■f . (114) 

as a consequence of the fluctuation-dissipation theorem. Here, V is the volume 
of the system and V^5k,k' = (27r)^5*^^)(k — k'). Equation (114) ensures that 

In general, the presence of the bath modifies the effective mass and couplings of 
the soft DCC modes and induces higher-order couplings as well [109,103]. However, 
these effects do not play a major role in the present qualitative discussion and we 
neglect them for simplicity. 
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the DCC eventually thermalizes with the heat bath at temperature T. For 
instance, one finds that, for large times: 

{^kW^k'W)B ^ (> ('k.k' ^—^ , (115) 



which indeed corresponds to the required thermal correlator. Finally, the fric- 
tion coefficient 7]^, which determines the rate at which equilibrium is ap- 
proached, is simply given by the on-shell damping rate at temperature T [109]: 
?7fc = 7fc = cr^ — cr^, where and denote the in-medium absorption and 
production rates for on-shell pionic excitations of momentum k respectively. 



2.3.2 The lifetime of a DCC excitation 

Averaging Eq. (113) over the noise, one obtains that the amphtude of the DCC 
field excitation decays as^^ (7ra(k, i))^ ~ Q-Vkt/2_ 'pj^g corresponding number 
of pions therefore decays with a rate 77^: 

UaiKt) = \l{7ra{Kt))^+UJk{7Ta{Kt))j' ~ e'^'^' . (116) 



Assuming a thermal bath of pions - which are the most abundantly produced 
particles in heavy-ion collisions, one can compute the relevant friction term 
?7fc. As a first estimate, the authors of Ref. [73] have computed the latter at 
lowest order in perturbation theory in the context of the linear sigma model 
in the high temperature, chirally symmetric phase, resulting in a rather short 
lifetime. For instance, they obtain, for the zero mode: tj^Iq ~ 1 fm/c, which 
is a characteristic scale of typical hadronic processes. However, as pointed out 
in [103], this estimate is too crude as it neglects the fact that, after formation, 
the DCC is supposed to evolve in the phase of spontaneously broken symme- 
try, where the dynamics is strongly affected by the presence long-wavelength 
Goldstone modes. In particular, the interactions of the latter are suppressed 
at low energies and the associated time-scales are expected to be correspond- 
ingly longer. Various possible contributions to the DCC decay rate have been 
studied in the literature [103,101,106]. We review the main results below. 

In Ref. [103], Rischke has investigated the damping of a DCC excitation in 
a thermal bath of pions in the context of the linear sigma model. In the 
broken phase, damping of DCC pions may arise from the absorption of a 
thermal pion, TTDccTTbath CTbath, together with the reverse process, cxbath 
TTDCcTi'bath, namely the decay of a thermal a. These processes, however, are 

This assumes that co^ > rjk/2, which is reasonable if the interactions with the 
heat bath are weak enough. 
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strongly suppressed due to the restricted available phase space [103]. The 
corresponding contribution to the damping rate reads: 



Vk 

"'k 



1 r ~ ~ 

{n;(l + n^)-n;(l + n;)}, (117) 



where we introduced four- momenta as e.g. /cf = (a;^., k) and we used the 
notation dp^ = ^2^)Auji ' ^ ~ Here uj\ — ^A;^ + mf and n\ — 

[expa;j./T — 1]~^ is the Bose- Einstein distribution for the respective species. 
The terms 1 + n in brackets correspond to Bose- enhancement factors in the 
final state. Finally, I^Vf^j^^^p denotes the relevant matrix element. The rate 
rjl'^'^ has been estimated in [103] using perturbation theory as a qualitative 
guide. At lowest order, one has IA^tttt^o-P = i'^d^aY, where gnia- — A/tt is 
the dimensionful ttttct coupling, and the momentum integrals in (117) can be 
easily computed [103]. Neglecting terms of relative order m^/m^, one obtains 
the following approximate expression: 



As expected, this results in rather long lifetimes for temperature of interest 
in heavy-ion collisions. For instance, one obtains l/ry^Z^J^'^ ~ 20 (9) fm/c for 
T = 150 (170) MeV. Of course, the result (118) should be regarded only as 
a rough qualitative estimate, being based on a perturbative expansion (recall 
that g ~ ma/2fTj k, 3). Still, this indicates that the DCC lifetime in the broken 
phase is considerably longer than typical hadronic time-scales. 

For a more quantitative estimate, it is important to include other processes, 
which are formally of higher-order in perturbation theory, but which can 
give significant contribution to the damping rates. The first such contribu- 
tion comes from two-body elastic scatterings tttt — > tttt. A non-perturbative 
estimate of the latter has been performed in Ref. [101], where the authors 
used experimental data to constrain the relevant matrix element. ^ The cor- 
responding contribution to the damping rate reads: 

An exact expression for rj'^Zo*'^ ^^^^ order of perturbation theory can be found 
in [103]. We mention that a more general expression, including non- vanishing mo- 
mentum, can be found in Ref. [106]. 
^^For a perturbative calculation, see Ref. [106]. 

^ We mention that the authors of Ref. [101] employ a slightly different approach 

to describe the damping of the DCC excitation. Indeed, they derive an expression 
for the DCC decay rate by generalizing the LSZ reduction formula to include the 
DCC coherent state in the initial and final states. The expression they obtain for 
the damping rate of a given mode k is identical to the one used here, that is the 
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Fig. 13. Various contributions to the decay rate r]k=o = —N/N of an homogeneous 
{k = 0) DCC excitation as a function of the temperature T of the surrounding heat 
bath. From Ref. [101] 



1 /■ _ _ _ 

77j^ = — - / dpdqdr {27ry6'-^\k + p- q-r)\M 



^7r7r<->7r7r 



TTTT — ^TTTT 



X 

2 



Un;{l + 0(1 + <) - «(! + n;)}, (119) 



where the factor 1/2 in the second Hne accounts for identical particles in the 
initial (final) state. The isospin averaged matrix element can be written as 
jX^^^^^p = I XI/ (2/ + l)|Mj^^p, with the standard partial- wave decompo- 
sition of the amphtude for scattering in a state of total isospin / = 0, 1, 2 
[101]: 

M^, = 327r ^(2^ + l)P^(cos 9) ^e'^i sin 5l , (120) 



where ^/s, q and 6 denote the center of mass energy, three-momentum transfer 
and scattering angle of the binary collision respectively. Using a parametriza- 
tion of the data for the three dominant phase shifts 5q, Si and 5q [110], and 
assuming, as above, that the DCC pions have essentially zero-momentum, 
the authors of Ref. [101] obtain the estimate l/rjll^^'^ fti 8.9 (5.6) fm/c for 
T = 150 (170) MeV. Again, as expected from low-energy theorems, this is 
much longer than typical hadronic time scales and is instead characteristic 
of Goldstone dynamics. Note also that the contribution to the damping rate 
from two-body scatterings is larger than the previous one from pion absorp- 
tion and sigma decay. Taking both contributions into account one obtains 
{rilZo'^ + ilk=o'^^)~^ ~ 6-1 (3.5) fm/c at the respective temperatures. 

usual on-shell damping rate (see Eq. (119) below). 
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Similarly, the authors of [101] have estimated the contributions from irp scat- 
tering, assuming that the latter is dominated by the formation of a ai(1240) 
meson in the s-channel, as well as from scatterings of DCC pions with baryons. 
While the former is found to give a negligible contribution due to kinematics, 
it was argued that the latter can play a substantial role since a number of 
baryons are produced even in the central rapidity region. The contribution 
of pion-nucleon scattering to the damping rate has a similar expression as 
Eq. (119), but replacing the combination of statistical factor on the second 
line (including the overall 1/2) by: 

<(1 + <)(1 - n^f) - <nf (1 - O , (121) 



where rip — [exp(£'p — ijl)/T + 1] ^ is the Fermi-Dirac distribution for nu- 

cleons, with Ep — ^Jp^ + m^. The baryon chemical potential ji is introduced 
to account for the finite baryon density in the central rapidity region and 
is chosen to reproduce the observed ratio between pions and nucleons. The 
authors of [101] use the value ii = 260 MeV at T = 150 MeV in order to 
enforce a 5 : 1 pion-to-nucleon ratio as observed at SPS energies [111].*^^ Us- 
ing, again, a parametrization of experimental data to describe the relevant 
scattering matrix elements, they find that the effect of nucleons is about a 
20% reduction of the DCC lifetime at the SPS. '^^ This is shown in Fig. 13, 
where the various contributions to the DCC decay rate are shown as a func- 
tion of temperature. The chemical potential is adjusted in order to maintain 
a fixed pion-to-nucleon ratio at all temperatures. The figure also shows the 
contribution from irp scattering. 

Finally, it is interesting to estimate the influence of the finite size of the DCC, 
that is to take into account the spread of the momentum distribution of the 
DCC pions. Adapting the treatment of Ref. [101], we define a momentum 
averaged lifetime f as: 

1_ ^(^^ IkVknjKt) 

f- N{t)^ /,n(k,t) ' ^^''^ 



where we used the notation = J Here, N{t) = n{k,t) is the total 
number of DCC pions, with n(k, t) = J2a^a{k,t) and na{k,t) the number of 
DCC pions of isospin a in mode k, as defined in (116). We used the fact that 

The pion-to-nucleon ratio quoted in Ref. [Ill] is actually of 6 to 1, corresponding 
to = 230 MeV at T = 150 MeV. However, this small difference has little impact 
on the results described below [101]. 

At RHIC energies, the pion-to-nucleon ratio tt/N is increased by about a factor 
2, but this is partly compensated by the increase of anti-nucleon production (see 
e.g. [112,113]). One has, roughly, Tr/{N + N)^7. 
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Fig. 14. The momentum-averaged DCC lifetime (122)-(123) as a function of the 
momentum spread Ak = cj of the DCC excitation, assuming a Gaussian shape (see 
text). The sohd hne shows the contribution from tttt scattering alone, whereas the 
dashed curve is obtained by including np and ttN scatterings as well. Prom Ref. [101] 



n(k, t) ^ —rikn(k,t) in the limit r)k ^ uj^ (see Eq. (113)) to approximate 
the RHS of Eq. (122). For a rough estimate, we write n(k, t) ^ A^(t)jF(k), 
where ^(k) describes the momentum shape of the DCC spectrum. Inserting 
this ansatz in (122), we get: 

^^-jr^kH^). (123) 

k 



This expression is similar to the one used by the authors of Ref. [101]. Their 
result is reproduced on Fig. 14, which shows the momentum averaged life- 
time (123) as a function of the momentum spread a = Ak of the DCC. 
This assumes a Gaussian shape, ^(k) oc exp[— A;^/2A/c^], and includes the 
contributions from tttt, vrp as well as nN binary scatterings (with baryonic 
chemical potential adjusted to maintain a fixed pion-to-nucleon ratio). The 
corresponding matrix elements are parametrized as above. Obviously, increas- 
ing the width Ak decreases the lifetime: Smaller DCC domains decay faster. 
The typical momentum spread predicted by the quench scenario discussed in 
previous sections is given by the symmetry breaking scale Ak < /j^ ~ 90 MeV, 
corresponding to an estimated average lifetime r > 5 fm/c at SPS energies. 
Although not very large, this is reasonably long for a hadronic DCC signal to 
be observable. Notice that larger DCCs, which have a better chance to survive 
until freeze-out, are also the rarest. 
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2-4 Phenomenology 



The possible detection of a DCC in high-energy hadronic or nuclear colli- 
sions^^ would bring valuable information concerning the chiral structure of 
the QCD vacuum and/or concerning the chiral phase transition. Although, as 
wc have seen in the previous sections, the present theoretical understanding 
of DCC formation is far from being accurate enough to produce reliable num- 
bers, it is important to characterize even qualitatively possible experimental 
signatures of the phenomenon. The DCC pions are expected to be mostly con- 
centrated at low momentum in the DCC rest frame and the DCC emission 
should, therefore, be characterized by a cluster of pions with low relative mo- 
menta. Ideally, the latter would present characteristic isospin fluctuations. In 
the present section, we describe various hadronic signatures which have been 
discussed in the literature as well as experimental tools used in actual DCC 
searches (see Sec. 3). We also discuss possible electromagnetic signatures. 

As already emphasized, the formation of a DCC is expected to be a rare 
phenomenon and the main difficulty is to isolate any DCC signal from the 
enormous background of "incoherent" pions, produced by standard mecha- 
nisms. In particular, the presence of DCC pions would hardly be detectable 
in global observables, such as e.g. the single-pion inclusive spectrum and 
event-by-event analysis seem preferable. We focus on the latter in the follow- 
ing. 



2.4-1 Multiplicity fluctuations 

A simple hadronic observable is the pion multiplicity in a given bin in mo- 
mentum-space. In particular, multiplicity fluctuations can reveal the dynamics 
underlying pion production. They can be characterized by means of the so- 
called factorial moments, defined as: 

TV' 

/, = {N{N - 1) • • • (A^ - i + 1)) = E P(^^W—)\_ ' (124) 



where P{N) denotes the multiplicity distribution in the considered bin of 
momentum-space. Clearly, the i-th factorial moment probes events where 

"^'^ We mention that the possibility of DCC production at a photon collider has been 
advocated in Ref. [114]. 

Multi-pion correlations have also been proposed as possible probes of the forma- 
tion of a DCC [49,50,51,115]. It is indeed well-known that the presence of a coherent 
source reduces the so-called Hanburry-Brown-Twiss (HBT) effect on identical pion 
correlations [52]. However, the interpretation of HBT measurement is highly non- 
trivial [116,52] and such effects are very difficult to observe in practice. 
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Fig. 15. Left: Reduced factorial moments for the multiplicity of positively charged 
pious emitted at mid rapidity iu central Au + Au collisions at the full RHIC energy, 
as obtained from semi-classical solutions of the linear sigma model in an expand- 
ing cylindrical geometry (triangles) [71]. Also shown are the results generated by 
HIJING with (squares) or without (circles) jet quenching and by UrQMD (dia- 
monds). Solid and open symbols correspond to soft {pT < 200 MeV) and hard 
{pt ^ 200 MeV) pious respectively. Right: The corresponding ratios of reduced 
factorial moments calculated for soft and hard pious. From Ref. [117]. 

the multiplicity N > i. Purely statistical fluctuations are characterized by a 
Poisson distribution, for which fi = (iV)*. It is therefore useful to introduce 
reduced factorial moments: 



Deviations from = 1 reflect the presence of dynamical (non-random) fluc- 
tuations. 

It is interesting to compute the reduced factorial moments in a specific model 
of DCC production. In Ref. [117], Bleicher et al. have performed such an anal- 
ysis in the context of a semi-classical treatment of the linear sigma model 
dynamics, in an expanding geometry with cylindrical symmetry [71]:^^ Us- 
ing the classical field theory methods described in Sec. 2.2, one computes 
the average occupation number na(k) in a given momentum mode k from an 
expression analogous to Eq. (49), suitably adapted to the expanding geome- 
try. Identifying each classical field configuration with a quantum-mechanical 
coherent state, the actual multiplicity is then sampled from the associated 
Poisson distribution with average given by na(k). The latter fluctuates from 
event to event and one finally performs a statistical average over initial field 

This is mainly relevant to the case of nuclear collisions. For analysis of multiplicity 
fluctuations in models based on coherent and/or squeezed states descriptions of the 
DCC, see e.g. Refs. [43,44,45,46]. 
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configurations in a thermal ensemble at a given initial temperature Tq. The 
corresponding reduced factorial moments are shown in the left panel of Fig. 15, 
for both soft (px < 200 MeV) and hard {px > 200 MeV) pions. Also shown 
are the corresponding results obtained from the event generators HIJING [118] 
and UrQMD [119,120], based on standard production mechanisms. The lat- 
ter clearly exhibit non-Poissonian fluctuations, with a gentle increase of the 
reduced factorial moments Fi with the order i, both for soft and hard pions. 
In contrast, the dynamics of the linear sigma model leads to reduced factorial 
moments that remain close to unity for the hard pions while increasing rapidly 
for the soft pions, indicating enhanced multiplicity fluctuations. 

The qualitative difference with standard production mechanisms is better em- 
phasized by considering the ratios between the reduced moments for soft and 
hard pions, as illustrated in the right panel of Fig. 15. These ratios may pro- 
vide a useful observable to search for dynamical fluctuations beyond what has 
been included in standard event generators. 

2.4-2 Charged-neutral fluctuations and robust observables 

The most striking feature of the coherent (semi-classical) nature of DCC emis- 
sion is the anomalously large event-by-event fluctuations of the neutral frac- 
tion of produced pions, see Eq. (10). This is at the basis of most existing DCC 
hunt strategies. Similar to the previous analysis of multiplicity, such charged- 
neutral multiplicity fluctuations can be conveniently studied by means of the 
bivariate factorial moments 

{nchirich - 1) • • • (rich - i + 1) no{no -l)---{no- j + 1)) , (126) 

where rich and no denote the charged and neutral-pion multiplicity respec- 
tively. Here we review the discussion of charged-neutral fluctuations presented 
in Ref. [121] for the analysis of the MiniMAX experiment at the Fermilab 
Tevatron. In particular, these authors have introduced specific observables 
(sec Eq. (143) below) which arc rather insensitive to the fiuctuations of the 
total multiplicity and, therefore, essentially probe the charged-to-neutral ra- 
tio fiuctuations for a wide class of production models. These so-called robust 
observables are also largely insensitive to detector efficiencies. 



2.4.2.1 Generating functionals for bivariate factorial moments Let 

p{nch, no) denote the probability for the occurrence of rich charged and no neu- 

Similar results have been obtained in a slightly different context in Ref. [104]. Such 
large multiplicity fluctuations are also already visible from the broad distributions 
of the amplification coefficient reported in Ref. [26] (see Fig. 7 of the present report). 
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tral pions in a multi-particle event, in a given region of phase-space. Equiva- 
lently, one introduces the generating function: 



G(^ch,^o)= E PKh,no);^r^o% (127) 

nch,no>0 



which encodes the bivariate factorial moments (126): 
rich! no! \ d'^^G{zch,zo 



(rich -i)! (no -j)!/ dzi^dz^ 



(128) 

2ch=20 = l 



Denoting by P{N) the total multiplicity distribution in the considered region 
of phase-space, one writes 

p(nch, no) = P{N)pN{nch, no) , (129) 



where j97v(nch,no) denotes the conditional probability that rich of the pro- 
duced pions be charged and Hq = N — rich be neutral. A wide class of models 
can be described by the following "binomial transform" [121]: 

Piv(neh,no) = ^^ / df p{f) ni - fr^K (130) 
nch'no! J 



For instance, the case of generic pion production can be described by a stan- 
dard binomial distribution, which corresponds to p{f) = 6{f — /), where 
/ = 1/3 is the mean fraction of neutral pions. Equation (130) also includes 
the case of ideal DCC production, which corresponds to p{f) — 1/2-^//, see 
Eq. (34). Using Eqs. (127), (129) and (130), one can write: 

1 

^(^ch, ^o) = / dfpif) Gifzo + (1 - /)^ch) (131) 



where 

G{z) = E P{N) (132) 

N>0 



is the generating function for the factorial moments of the total multiplicity 
fluctuations. 
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2.4.2.2 Charged-pions— photons fluctuations Neutral pions are detec- 
ted through their decay in a pair of photons vr'^ — > 77. It is, therefore, useful 
to measure directly the bivariate factorial moments for charged-pions-photons 
fluctuations. The latter can be defined as in Eq. (126) with the replacement 
no n-y-i the number of detected photons in the region of momentum-space 
under consideration, and d/dzo — > d/dz^. For ideal detection efficiency, the 
corresponding generating functional can be obtained from Eq. (131) by re- 
placing zq — > z'i^. More realistically, there are probabilities £9,1,2 that 0, 1, or 2 
photons from the 7r° decay be actually detected, with eo + ei + £2 = 1. These 
three possibilities can be taken into account by introducing the following gen- 
erating function: 

go{z^) = eo + eiz^ + e2Zch ■ (133) 



Similarly, taking into account the detection efficiency for charged pions, there 
is a probability that a charged pion be actually observed in the detector 
and a probability 1 — e^h. for not observing it. As before, these two possibilities 
can be described by the following generating function: 

5'ch(^ch) = (1 - Cch) + ech^^ch • (134) 



Thus the relevant generating functional for observed charged-pions-photons 
fluctuations, taking into account flnite detection efficiencies, can be obtained 
from Eq. (131) with the replacement Zch. — 5'ch(-2ch) and zq goiz-y) [121]: 

Gobs(^ch, %) = G{gc\y{zch)igo{z^)) ■ (135) 



One has, for the observed bivariate moments: 



II- I 



{rich - «)! (% -i)!/°bs 



■9^+-^'Gob.(.-ch.%) 
dzi^dz^^ 



^ch — ^'Y 



(136) 



For instance, one easily obtains, for the first moments: 



(^7)obs 


= (61 + 262) (/) (N) 


(137) 


(^ch)obs 


= ech(l-/) (N) 


(138) 


('^7^ch)obs 


= ech(ei + 2e2) (/(I - /)) (7V(7V - 1)) 


(139) 


{n^(n^ - l))obs 


^(e, + 2e,r{f){N{N -1)) + 2e,{f){N) 


(140) 


(rich (rich - l))obs 




(141) 



where the brackets on the RHS denote an average with respect either to the 
multiphcity distribution P{N) or to the distribution p{f) (cf. Eq. (130)). 
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Prom the above expressions, one observes that the ratio 
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^'^ (nch(nch - l))obs ("'7)obs ((1 - fy) 



is independent of the efficiencies introduced above as well as of the multiplicity 
distribution and only depends on the distribution p{f). More generally, it is 
easy to check that the following ratios: 

ri,i = ^, (143) 



of reduced bivariate factorial moments: 



("di("di -!)••• (/'di - / + -!)••• ("-;, - J + l))c 



\?T'ch/obs\%/obs 



are robust observables in the above sense: 

^ ^ (/(I - /)■> 



Moreover, these ratios are sensitive to the difference between generic and DCC 
production. For generic production, (/') = (/)' = (1/3)* and one immediately 
gets: 

ri,i (generic) = 1 (146) 



for all i's. In contrast, for ideal DCC production, one has: 

if (1 - /) )dcc - 2r(z + j + |) ' ^^^^^ 



from which one obtains [121]: 

r,i(DCC) = ^. (148) 



The increasing difference between the generic and ideal DCC production mech- 
anisms with increasing order i of the ratios rj^i clearly reflects the broadness 
of the charged-to-neutral distribution in the latter case. 



Here and in the following, we assume exact isospin symmetry, that is (/) = 
Jo dfp{f)f = 1/3, for simplicity. 
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2.4.2.3 Mixed events The above discussion concerns the cases of pure 
generic versus pure DCC production. More reahstic descriptions should include 
intermediate situations, where both generic and DCC production mechanisms 
contribute. Various types of such mixed events have been studied in [121]. 
Here, we consider a model where a DCC is formed with probability a and, 
when this happens, the detected pions are either of generic or of DCC origin. 
The corresponding generating function is given by the product: 

Gmix(^ch, Zq) = Goccl^ch, ^o) G^gen(^ch, ^o) ■ (149) 



where Goccl-^ch) -^o) is given by the generic form Eqs. (131)-(132) with multi- 
phcity distribution Pdcc(^) and withpDcc(/) = 1/2a/7. Similarly, G'gcn(-2ch, zq) 
corresponds to a given multiplicity distribution Pgen(-/V) and Pgen(/) = ^{f — 
1/3). In the cases where no DCC is formed, which occur with probability 1 — a, 
the production is purely of generic type. The total generating function takes 
the form: ^® 

G{Zch, Zq) = aGmi^{Zch, Zo) + (1 - a)Ggen{Zch, Zq) . (150) 



From the definitions (143) and (144), and using (128) and (150) above, one 
can compute the robust ratios r^^i in the present model. For instance, one 
finds: 

r,,^l±4^, (151) 



l + a{A + 2B) ' 



where 



{N{N - l))ge„ ^ ^ 

^^3(iV(iV-l)W_ (153) 



5 {N{N-1)) 



gen 



Here, the brackets denote averages with respect to the multiplicity distribu- 
tions corresponding to either the DCC or the generic production mechanisms, 
with obvious notations. The case of pure generic production is obviously re- 
covered for a — 0. Assuming that the probabihty of DCC formation is small. 



Taking into account the detection efficiencies as in the previous discussion, the 
actual generating function is still given by Eq. (135), with the RHS corresponding 
to Eq. (150). It is easy to check that the ratios (143) remain independent of the 
efficiencies. 
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Fig. 16. The robust ratios rj^i for the model described in the text (cf. Eq. (155)) 
with, from top to bottom: i = 1,...,5, as a function of the fraction /3 of DCC 

pions in the considered region of phase space. The probabihty that a DCC forms is 
10^^. We have assumed that the total multiplicity distributions Pgen{N) and 



a 



Pdcc{N) are both Poissonian for simplicity. 

a <^ 1, one obtains the approximate expression: 

3a {N{N - 1))dcc 



5 {N{N-1)) 



(154) 



gen 



Similarly, one can compute the 0{a) deviation from pure generic production 
due to the presence of a DCC for higher-order ratios. One obtains, after some 
calculations: 



i< k (3/3)^« Fff^^Ffn 



where j3 = (A^)Dcc/(^)gen is the fraction of DCC pions in the region of 
phase-space considered and the F^'s denote the reduced factorial moments, 
cf. Eq. (125), corresponding either to the DCC or to the generic multiplicity 
distributions. As expected, the deviation from i = 1 is larger for higher- 
order ratios and increases as the probability of DCC formation and/or as the 
fraction /3 of DCC pions increase. This is illustrated in Fig. 16, where we 
show the ratios rj^i for i = 1, . . . , 5 as a function of [3 for a given value of the 
probability a. For illustration, we have taken both the multiplicity distribu- 
tions Pgcn(^) and -Pdcc(-^) to be Poissonian, that is all the reduced factorial 
moments = pj^cc — i (155). This constitute a lower bound on the 
expected deviation from 1 for the robust ratios rj^i. 
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The results described here assume the general form (130) for the multiplicity 
distribution. This includes ideal DCC production as well as simple modeliza- 
tions of generic production. It would be interesting to extend the analysis of 
robust ratios to dynamical scenarios of DCC formation, such as e.g. semi- 
classical simulations of the out-of-equilibrium chiral phase transition, as de- 
scribed in previous sections. 



2.4-3 Multi-resolution discrete wavelet analysis 

As described above, a possible DCC signal would appear as a localized struc- 
ture in momentum space. To search for such structure in the lego plot, the au- 
thors of Ref. [122] have proposed to use the so-called multi-resolution wavelet 
analysis. This provides a useful tool to identify possible structures simulta- 
neously in term of their size and location in momentum space. This generalizes 
the standard Fourier analysis. The technique has been applied to classical sim- 
ulations of the linear sigma model [122,124]. It has also been used by various 
experimental collaborations (see Sec. 3). Here, we review the basics of the 
discrete wavelet analysis. 

To illustrate the method, we consider a one- dimensional phase-space, described 
by the dimensionless variable x in the interval [0, 1]. Any observable can be 
represented by a function f{x) on this interval. Let us cut the accessible phase- 
space in 2^ bins of size Ax — 1/2^ , where j is a positive integer j < jmax, with 
jmax corresponding to the finest accessible resolution. To each collision event 
corresponds a given sample of the function f{x). Call fl = f{x = kAx) the 
value of the sampled function in the kth bin, the complete sample can be 
represented by the following function: 

fe=0 



where 0j.(x) is zero everywhere but in the kth bin: 

lk/2^ <x<{k + l)/2^ 
(Piix) = { . (157) 

otherwise 

Notice that the "bin functions" (/)l.{x) can be obtained as translations and 
dilation of a single function (/){x), called the mother function: 

4(x) = <f>(2^x - k) , (158) 
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For a review see e.g. [123]. 
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where, in the present case, = ^(1 — x)9{x). 

Equation (156) describes the sampled observable at the resolution scale j. To 
detect structures at lower resolution scales, one can apply the following coarse- 
graining procedure: Replace two adjacent bins 2k and 2A; + 1 by a single bin 
of size 2Ax — l/2^~^, with corresponding bin function given by: 



and define the value fl ^ of the function in the new bin k as the average of 
the values in the previous smaller bins: 

ft'-kfL + fk+i)- (160) 



The resulting coarse-grained sample at scale j — 1 is simply given by: 

f'-'\x) J'i:' ft' <i>i-\x) . (161) 



fe=0 



Equation (161) defines the mother function representation of the distribution 
f{x) at scale j — 1. The /^^'s are accordingly called the mother function 
coefficients (MFCs). Clearly the procedure described above eliminates all fluc- 
tuations at shorter scales than the resolution scale j — I. These steps can be 
repeated to search for possible structures from the finest experimental resolu- 
tion scale jmax to the lowest one j = 0. 

The information lost at each step in the coarse-graining procedure is encoded 
in the difference /^•'"^^(x) = f^\x) — f'^^'^^x). It is easy to show that the 
latter can be given a similar representation: 

f^'-'\x)^'i:'fi-'^i-\x), (162) 

k=0 

where the functions V'fc"^ ^-^^ given by: 

V'r'(^) = 02.(^)-0Wi(^)- (163) 



They can be obtained as translations and dilations of the so-called father 
function: ipl{x) = ip{2^x — k), where, in the present case, ip{x) = (f){2x) — 
(f){2x — 1). The corresponding father function coefficients (FFCs) in (162) are 
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Fig. 17. Scale dependence of the discrete wavelet power spectrum for various dy- 
namically generated data. Prom Ref. [122]. 

related to the MFCs at the previous resolution scale j through: 

fl-' - lifL - ■ (164) 



Equation (162) defines the father function representation of the sampled dis- 
tribution f{x) at the scale j — 1. Clearly, the FFCs (164) at a given scale 
j measure the variation of the sampled distribution / between two adjacent 
bins. The information about the size of the fluctuations at each scale can be 
conveniently encoded in the following power spectrum: 

^ fe=0 

Similar to the Parseval theorem for the Fourier transform, one has: 

1 

^P,= fdx\f{x)f. (166) 

j>0 



To demonstrate the method, the authors of [122] have applied the above analy- 
sis to various dynamically generated data. They have considered the evolution 

The sample at scale j can be fully reconstructed from the mother function rep- 
resentation at lower resolution scale j — k with the help of all father function rep- 
resentations at the intermediate scales j — k, . . . ,j: 

k 

m=0 
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of the pion field as described by the hnear sigma model in the Hartree approx- 
imation, for a boost invariant 1 + 1 expanding geometry. The initial field con- 
figurations are sampled from the quenched ensemble described in Sec. 2.2.1. 
The discrete wavelet analysis is applied to the ratio: 



where vri.2,3(?7) the isospin Cartesian components of the pion field as a 
function of the spatial rapidity i] = ^ \n{t + z)/(t — z), where z is the collision 
axis. The ratio (167) is related to the neutral pion fraction discussed in previ- 
ous sections. The rapidity-interval —rj^ax < i] < rjmax, with rjmax = 5, is bined 
with a bin size Ar) — 0.08, corresponding to a finest resolution scale j^ax — 7, 
such that 2r/jnax/^'7 = 2^"'°-'^. Figure 17 shows the corresponding power spec- 
trum (166) at various scales, for both the initial quenched ensemble (random 
sample) and the final time evolved ensemble (DCC sample). Also shown are 
the results corresponding to a thermal ensemble as well as the power spectrum 
obtained from HIJING Monte Carlo data. The DCC sample clearly exhibit 
a non-trivial structure, refiecting clustering in spatial rapidity. Similar results 
have been obtained in Ref. [124] for three-dimensional classical simulations of 
the linear sigma model in a static geometry with a cooling term, as described 
in subsection 2.2.2 above. 

It is important to mention that the analysis of Refs. [122,124] concerns data 
sampled in coordinate space. Although this is sufficient to demonstrate the 
wavelet analysis technique, it would be of interest, for phenomenological pur- 
poses, to perform similar analysis in momentum space. 



2.4-4 Electromagnetic signatures 

Because they interact strongly with the background of other hadrons, hadronic 
signals may rapidly loose the information they carry before they leave the 
interaction region. In contrast, electromagnetic signals do not undergo multiple 
rescatterings and, therefore, provide a sensitive probe of the various stages of 
the collision. If the DCC lives long enough, they might probe its formation 
even if the latter is diluted in the hadronic environment before freeze-out. 

It has been proposed that anomalous pion production from a DCC would result 
in anomalous radiation of low- momentum photon and/or low-mass dileptons 
[125,126,127,128,129]. Here we illustrate these ideas in the case of dilepton 
production. At lowest order in the electromagnetic coupling constant aem, the 

In that case the sampled observable is the fraction of neutral pions as a function 
of rapidity. 
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dilepton production rate reads: 



{q^^q- -q^g^^-)W^,{q), (168) 



where q^ is the four-momentum of the lepton pair, such that q^ > 4m|, with 
rrii the lepton mass, -B(g^) = (1 + 2m^/g^)(l — Aml/q^Y^'^ and W^^^'^q) is the 
4-dimensional Fourier transform of the electromagnetic current correlator: 

lyM^(g) = / d^xd^ye-''i-^''-y\j^'{x).r{y)) , (169) 



where the brackets represent the average in the state under consideration, and 
J^{x) is the electromagnetic current operator: 



Ji'ix) = ^[7rt(x) d^" 7r(x) - 7r(x) t:\x)\ , (170) 



where fd^g = fd^g — {dnf)g and where tt^x) and ti\x) denote the charged 
pion field operators. 

To illustrate the anomalous electromagnetic emission due to a DCC, we es- 
sentially follow the schematic model presented in [126,127]. The presence of 
a DCC may be modeled by splitting the pion field operators into a classical 
part 7rci(a;), which represent the coherent DCC excitation, and the rest 7f(a;), 
which describes the background pions: 

'k{x) = -Kclix) + Tx{x) (171) 

The electromagnetic current (170) can be written as a sum of three contribu- 
tions: 

r{x) = J^,{x) + J^{x) + J^{x) , (172) 



where J^iix) is the purely classical contribution, obtained from Eq. (170) with 
the replacement vr tTci, and Jn=i 2(2^) contains n = 1, 2 powers of the field tt 
respectively. Inserting this decomposition in Eq. (169) and assuming that the 
system of background pions is electrically neutral, one obtains: 

W>^^{q) = WS'iq) + Wr{q) + Wr{q) , (173) 

where the first term on the RHS, 

Wi:^{q) = Ji:,{q)J^{{<D (174) 
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describes the purely classical contribution and 



(175) 



with J^(g) and J^(g) the four-dimensional Fourier transform of J^(a;) and 
J^{x) respectively. 

The contribution to the rate (168) from the first term on the RHS on Eq. (173) 
corresponds to the anomalous bremsstrahlung from the DCC itself. As an 
illustrative example [126], it can be estimated using the Blaizot-Krzywicki 
classical solution [2], described in Sec. 2.1. For this purpose, it is useful to 
notice that the electromagnetic current (170) actually coincides with the third 

— * 

component of the iso- vector current — if x d^^Tf: 

J^(x-) = TTi{x)d^Tr2{x) - 7r2(a;)a^7ri(a;) = Vs^{x) . (176) 



For the one-dimensional boost-invariant solution of Refs. [2,42], one has 7i{x) = 
7f(T), with proper-time r = y/t^ — z^, and the classical iso- vector current can 
be written as (cf. (15)): 

KM=K'^^^(r-ro)^(x^) (177) 



where u'^ — (coshr^, sinhr^, 0, 0), with spatial rapidity r) defined as tanhr^ = 
z/t. Here, following the authors of Ref. [126], we have multiplied the boost- 
invariant solution of Blaizot and Krzywicki by a slowly varying function of 
the transverse coordinates g{'x.±) to account for the finite transverse extent of 
the system. The Fourier transform of the current (177) is easily evaluated and 

one obtains, after averaging over the possible orientations of the iso-vector a 
[126]: ^2 

- ""^4) [JoiM^ro) + N^oiM^ro)] , (178) 



dydMd?q^ 247r^ ^' MlM^ 



where M = is the invariant mass of the lepton pair, y and its longitu- 
dinal rapidity and transverse momentum respectively, and = \J q\ + M^. 
Here, Jq and A^o are usual Bessel functions and ^(q±) is the Fourier transform 
in transverse space of the shape function g{yi±). High transverse momentum 
dilepton are suppressed due to the finite transverse extent of the source. For 
instance, for a Gaussian shape 5'(xx) = exp(— x^/i?^), one has a strong expo- 
nential suppression: |^(q±)P oc exp(— g^7?^/2). Notice that, as a consequence 

A similar expression can be obtained for the spectrum of real photon produced 
by the DCC classical field, see Ref. [126]. 
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of boost-invariance in the present model, the rate (178) docs not depend on 
the dilepton rapidity y. Finally, {a\) measures the fluctuations of the strength 
of the iso-vector current in the initial state. Neglecting the pion mass. The 
latter can be related to the mean initial energy density (eo) using Eq. (22). 
Assuming a high temperature symmetric initial state where the values of the 
fields (tt, a) and their proper-time derivatives at r = tq are Gaussian random 
numbers of zero mean and of variance ui and (T2 respectively, as in (25), one 
obtains: 

(al) = 2«:^ = |(eo). (179) 



with — (7i(T2To. 



The dependence of the yield (178) on the dilepton mass at fixed transverse 
momentum is shown by the dotted curve in Fig. 18. Also shown are the con- 
tributions from the annihilations of on-shell background pions, assuming that 
the latter form a thermal bath (dashed curve), as well as the contribution 
arising from the interactions between thermal and DCC pions as described 
below (dot-dashed curve). One observes that the bremsstrahlung contribution 
(178) from off-shell DCC pions have no threshold and is actually mainly con- 
centrated at low invariant mass M < 2m,r- In this region, however, the yield is 
extremely sensitive to the initial time tq due to the singular behavior of Bessel 
Neumann function Nq{z). Moreover, this signal would probably be very diffi- 
cult to discriminate in a heavy-ion collision due to the large background from 
Dalitz decays 7r° — > e"'"e~7 in this kinematical region (see e.g. [130]). 

Another possibility, suggested in Ref. [127], is that non-coherent pions from 
the surrounding heat bath may annihilate on the coherent DCC state, forming 
dilepton pairs of typical invariant masses in a rather narrow window around 
M ~ 2m^, thus avoiding the 7r° Dalitz-decay region. The corresponding con- 
tribution to the dilepton rate (168) comes from the second term on the RHS 
side of Eq. (173). Assuming a thermal bath of non-interacting pions at tem- 
perature T, one obtains, after some calculations, the following contribution 
[127]: 



(1) 



2al 



-[ 
M J 



dydMd'^q± Stt^ MJ 2u;fe(27r)3 

{(l + n^)|7rei(r)|2 + (l + n+) |7rci(-r 



+ % Fci( 



+ nZ kci( 



-or}- 



(180) 
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The energy density seems to be over-estimated by a factor 4 in [126]. 
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Fig. 18. The e+e^ yield (168) as a function of the dilepton invariant mass. The 
dotted line corresponds to the anomalous bremsstrahlung from the DCC classical 
field configuration; The dashed line represents the contribution from the surround- 
ing thermal bath at temperature T: Finally, the dot-dashed curve represents the 
contribution arising from the interactions between DCC and thermal pions. The 
relevant parameters (see text) are: T = 145 MeV, (eo) = 58 MeV/fm^, tq = 1 fm. 
The absolute normalization of the DCC contribution depends on the details of 
the shape function g{xj_) for the DCC and is given here for a Gaussian shape 
g{^±) = exp(— x^/i?^) with R± = 2 fm. The relative normalization of the various 
contributions are not realistic (see text). Prom Ref. [127]. 

where Uk = \Jk'^ + ^1, = (t^fc,k), r'^ = (/c + q)'^, = {k — q)^ and 
denotes the equilibrium Bose-Einstein distribution at temperature T for 
positive and negative pions. The first two terms in curved brackets on the 
RHS of (180) correspond to the decay of the classical DCC field into an on- 
shell thermal pion and a dilepton pair: tTci — > Tti'^i~, whereas the last two 
terms represent the annihilation of a DCC pion against a thermal one: tTcitt — > 
The latter contribution has a threshold at M = 2m^, and produces 
an enhancement of the dilepton yield around M ~ 2m7r, as is clearly visible 
on Fig. 18. A more detailed analysis has been performed in Ref. [127], using 
classical field dynamics, which indicates that this enhancement is restricted 
to dilpeton momenta |q| < 300 — 500 MeV. 

We would like to stress that Fig. 18 is only meant here to illustrate the depen- 
dence of the various contributions to the dilepton yield with the dilepton mass. 
The relative normalization of the various contributions should not be taken 
too seriously as it is most probably far from realistic. First, because the present 
model completely neglects the small probability that a DCC actually forms. 
In a more realistic description, the rates (178) and (180) should be accordingly 
rescaled. Second, because the parameters used here correspond to an unex- 
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pectedly large contribution from the DCC. Indeed, as a rough estimate of the 
energy density of the background thermal pion gas, one may use the Stefan- 
Boltzmann law for a gas of non-interacting massless pions: e^r = g{n^/30)T'^ 
with g — 3 the pion degeneracy. For T — 145 MeV used in Fig. 18, one obtains 
~ 55 MeV/fm^ which is similar in magnitude as the initial energy density of 
the DCC state in the above caknilation: (eo) = 58 MeV/fm^. This is certainly 
unrealistic, at least in the context of the non-equilibrium chiral phase transi- 
tion scenario discussed in earlier sections. In a more realistic situations, one 
expects (eo) <C e,r and the various contributions to the dilepton rate shown in 
Fig. 18 should be appropriately rescaled. 

Finally, we mention that the generic dilepton spectrum in this region is dom- 
inated by T] Dalitz-decays, thus making the measurement of dileptons from 
a DCC a very difficult task. The possibihty of such a measurement would 
certainly require making appropriate cut on the data, probably in correlation 
with other hadronic observables. This certainly deserves further investigations. 

2.4-5 Further speculations 

It is usually believed that, in the context of high-energy nuclear collisions, the 
best place for DCC searches is the central rapidity region of the most central 
collisions, because this is where one expects the deposited energy density to 
be the highest, thereby providing the most favorable conditions for chiral 
symmetry restoration. However, as we have seen previously, this is also a place 
where DCC signatures might be difficult to observe. Here, we briefly mention 
alternative suggestions. 

Rajagopal has argued in Ref. [131] that semi-peripheral collisions at RHIC en- 
ergies might provide more favorable conditions for efficient quenching as well 
as sooner freeze-out, hence a better chance that DCC pions be visible in the 
detectors. Another interesting motivation for looking at such collision events 
has been proposed by the authors of Ref. [132,99] , who argued that strong elec- 
tromagnetic fields could result in a coherent excitation of the ttq field through 
the chiral U{1) anomaly. Such an excitation could be subsequently amplified 
by the quench mechanism, as discussed at length in the present report. More 
peripheral collisions might also be interesting for DCC searches because they 
essentially consists in the collision of the pion clouds of the respective nuclei. 
The relevant dynamics is, therefore, that of low-energy pion fields and the 
corresponding physical picture might be closer to the original baked-Alaska 
scenario (see Sec. 2.1). One could further imagine to favor a particular isospin 
direction by colliding asymmetric (e.g. neutron rich) nuclei. 

Another place of special interest to look for DCC signatures is the projectile 
This suggestion is due to D. Vautherin (private communication). 
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or target fragmentation region of nuclear collisions, as advocated by Bjorken 
[133,134]. The basic idea is that, viewing the collision in the rest frame of one 
of the ions, the projectile traverses the target at essentially the speed of light 
and ejects all the resting nuclear matter away almost instantaneously. If the 
smoking ruins left behind relax sufficiently slowly toward ordinary vacuum, 
the relevant dynamics should be that of soft pion fields and relaxation process 
should be characterized by low-momentum pion emission. 

It is also important to stress that present theoretical calculations in the con- 
text of the quench scenario in heavy- ion collisions indicate that the DCC might 
actually be more difficult to observe than originally thought [96,94,95,97,27] 
(see subsection 2.2.5). In particular, the results of Ref. [97] suggest that one 
should try to isolate as much as possible the various Fourier modes of the DCC 
pion field. Multi-resolution analysis techniques such as the discrete wavelets 
described above provide very efficient tools for such task. Moreover, in re- 
gards to the rather small expected probability that a DCC forms, it would be 
desirable to perform appropriate cuts on the data sample. For instance, one 
could restrict the analysis to events with a particularly large pion multiplicity, 
say larger than the average, in a given region of momentum space. We shall 
discuss such type of analysis in the experimental part of this review. 

Finally, we mention that a number of other observables, which we have not 
described here, have been discussed in the literature. This includes e.g. non- 
trivial charged-pion fluctuations [135] - which present the advantage that one 
does not have to reconstruct the ttq's, or anomalous strangeness production 
[61,136], etc. 



3 Experiment 

3.1 DCC-search experiments 

Can a DCC domain be produced in actual hadron or heavy-ion collisions 
and, if so, are we capable of detecting them? In the search for the answers to 
these questions, several experiments have been carried out and some have been 
planned for future. Some of them are based on the study of cosmic ray showers 
and others are accelerator-based. The list of all these experiments is given in 
Table 1. Below we discuss some of these experiments and their capability to 
look for DCC in terms of measurement of right observables. The results from 
these experiments will be discussed in detail in Sec. 3.4. 
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Table 1 

Experiments related to search for disoriented chiral condensates. 



Year 


Type 


Experiment 


V 


Observable 


1972-80 


Cosmic Ray, Balloons 


JACEE [6] 


> 1.7 TeV 


No. of 7 &: c/i in 5 < T/ < 9 


1980 


ClosTTiir Rriv 
Mt. Chacaltaya & Fuji 


RrRzil-.Tanan 


> 1.7 TeV 


No of ^ ^ ch in 
forward rj 


1977-91 


Cosmic Ray, Mt. Pamir 


PAMIR [7] 


> 1.7 TeV 


No. of 7 &; c/i in forward rj 


1982 




TJA5 \1A 151 


540 GeV 


IVn of 'V Ht' rh in i-^j A'k 

l.y\J, KJi. J txj V 111 ' ^ ^ / I 


and 


nucleon-nucleon 




and 


& asymmetry in electro- 


1986 


collision at CERN 




900 GeV 


magnetic & hadronic energy 


1983 


nucleon-nucleon 
collision at CERN 


TJAl fl3l 

KJ Xi. J. J. O 


540 GeV 


^mvTTiTTipfrv iTi olfpfrnTTi^i CTlPflP 

& hadronic energy 
in -3 < r? < 3 


1996 


Accelerator based, 


DO 


1.8 TeV 


asymmetry in electromagnetic 




nucleon-nucleon 


and 




& hadronic energy 




collision at FERMILAB 


CDF fl6l 

\_/J— /-I- 1 -I. 1 




in -4 < 77 < 4 


1997 


- do- 


MiniMAX [17] 


1.8 TeV 


No. of 7 &; c/i, in 
3 2 < rj < 4 2 


1998 


Accelerator based, 


WA98 [18,20,21] 


17.3 GeV/n 


No. of 7 & c/i in 


and 


nucleus- nucleus 






2.9 < ?7 < 3.75 


2001 


collision at CERN 








1999 


- do- 


NA49 [19] 


17.3 GeV/n 


Pt fluctuation in charged 
particles in 4.0 < r/ < 5.5 


2001 


Accelerator based. 


PHENIX [22] 


200 GeV/n 


No. of 7 & c/i, in 


and 


nucleus-nucleus 






-0.35 < 77 < 0.35 


2002 


collision at RHIC 









3.1.1 Cosmic ray experiments 

Cosmic ray experiments were performed by exposing stacks of nuclear emul- 
sion chambers at high mountain tops or sending them to high altitudes in 
balloons. They recorded secondary particles (mostly pions and gamma rays) 
originating in the successive interactions of a particle of primary cosmic radia- 
tion (such as proton) with atmospheric nuclei. The primary cosmic ray particle 
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Fig. 19. The rj — (p phase space distribution of an event from the JACEE collabo- 
ration [6]. The region indicated by the circle shows a region with anomalously high 
ratio of gammas to charged secondaries, called as anti-Centauro event. 



energies are in the range of 100 TeV to 10000 TeV [6]. A cosmic ray event can 
therefore be thought of as a result of collision between a projectile nucleon 
of an estimated energy of about 1500 TeV (equivalent center of mass energy 
^/s > 1700 GeV) and an atmospheric nuclei [15]. The mountain top experi- 
ments, carried out at Mt. Chacaltaya in Bolivia, at Pamir in Russia and Mt. 
Fuji in Japan, reported events having large number of charged hadrons (about 
63 - 90) and very few or no photons. These events are known as "Centauro- 
type" events. The first Centauro-type event called "Centauro-I" was reported 
in 1972. The JACEE collaboration [6], which had sent emulsion chambers 
in balloons to high altitudes also observed events having anomalously high 
number of gammas in comparison to charged secondaries. Such events are 
called "anti-Centauro-type" events. One such interesting event reported by 
the JACEE collaboration [6] is shown in Fig. 19, where the distributions of 
charged particles and gammas are presented in rj — (p space. One of the pos- 
sible explanations for these peculiar events with large isospin fluctuations in 
cosmic ray experiments is the possibile production of a DCC in high-energy 
collisions. 

There has been critical reviews of the exotic characteristics of Centauro-I. In 
Refs. [137,138] the authors found that there is no original correspondence 
between the clusters in upper and lower chambers of the emulsion chamber 
experiment. The authors of Ref. [137] claim that such an event no longer 
required any new physics to describe the event structure. In Ref. [138] there is 
a suggestion that the event structure may have arisen due a bundle of target 
interactions. They however claim that the characteristics, like hadrons not 
accompained by any 7 rays, can still not be understood and are very different 
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from commonly observed cosmic ray events. 

Recent results from a systematic study of the neutral pion fraction distri- 
bution (using sensitive techniques developed for DCC search) in high-energy 
cosmic ray families detected at high altitudes at Pamir [7] do not exclude such 
a possibility. The Pamir experiment took data by exposing thick lead cham- 
bers of homogeneous structure at Pamir in 1977-1991. The chamber thickness 
was between 40 - 110 cm Pb and, hence, had a detection probability of 7- 
hadron families close to 1. One of the limitations of cosmic ray experiments 
is that, because of the high energy detection threshold in experiments with 
emulsion chambers, the investigation of high-energy interaction is limited to 
the projectile particle fragmentation region. The situation is opposite to that 
in collider experiments where secondary particles emitted with small angles 
in far forward regions of the coUision cannot be detected because of the beam 
pipe. 



3.1.2 Nudeon-Nudeon collisions 

The observation of Centauro-type events motivated experimentalists to look 
for such events in accelerator-based experiments, even before the concept of the 
DCC came into picture. The UAl experiment [13] searched for events having 
characteristics of cosmic ray Ccntauros in p—p collisions at ^/s = 540 GeV, us- 
ing information on charged particle multiplicities, transverse momenta and the 
criteria of large hadronic and low electromagnetic energy content in calorime- 
ters. The central image chamber detector of the UAl apparatus was used to 
get charged particle tracks. Momenta were obtained from the curvature of the 
tracks in uniform magnetic field of 0.56 Tesla. The energy information was 
obtained from layers of finely divided sheets of lead and plastic scintillators 
which formed the electromagnetic calorimeter. The calorimeter surrounded the 
central detector. The electromagnetic energy deposition was measured from 
the first 4 radiation length of material. The hadronic energy contribution was 
estimated from the rest of the electromagnetic calorimeter and from energy de- 
position in the iron scintillator calorimeter behind it. No Centauro-like events 
were seen out of a total data sample of 48K events. Based on this the UAl 
collaboration reported an upper limit for the production cross section of the 
order of a micro-barn for a Centauro events at collider energies equivalent to 
155 TeV incident energy on stationary nuclei. 

The UA5 experiment [14,15], through an analysis of 3600 minimum bias events 
for correlation between photon and charged-particle multiplicities, reported no 

events having features corresponding to the Centauro events as observed in 
cosmic-ray emulsion chambers. The detector used for charged particle mea- 
surement was a streamer chamber, with a coverage of almost the entire 47r 
solid angle. The photon measurements were based on electromagnetic showers 
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resulting from photons converting throTigh the walls of the vacuum chamber 
(0.05 to 2.6 radiation length of material) situated between the two streamer 
chambers. The search were carried out at \/s energies of 540 GeV and 900 
GeV without any success. The search was extended to high ^/s energy of 1.8 
TeV in the DO and CDF [16] experiments at Fermilab. These experiments 
used the technique of the asymmetry in hadronic to electromagnetic energies 
for DCC search. They also reported null results. 

Subsequently a dedicated experiment at Fermilab Tevatron, called MiniMAX 
[17], was designed by Bjorken, Kowalsky, Taylor and collaborators to look 
for DCC signals in the far forward rapidity region oi p — p collisions at 
^/s = 1.8 TeV. This experiment, located at the CO area of the Fermilab 
Tevatron was designed to measure the ratio of charged-to-neutral pions. The 
acceptance in the lego space of pseudo-rapidity r] and azimuthal angle (p is 
roughly a circle of radius 0.65 centered at = 4.1. The MiniMAX apparatus 
included Multi-Wire Proportional Chambers (MWPC) to measure charged 
particles, movable 1 radiation length (Xq) Pb converter to improve the resolu- 
tion of charged-particles, a trigger scintillator, an array of 28 lead-scintillators 
acting as electromagnetic calorimeter located behind the MWPC to measure 
photons and upstream tagging detectors. The analysis of the data with about 
8 X 10^ events was done using the technique of robust observables, discussed 
earlier in Sec. 2.4, to take care of the small acceptance of the detector and 
various efficiency factors. The analysis yielded results consistent with no DCC 
production. 



3.1.3 Nucleus- Nucleus collisions 

Various DCC searches have also been pursued in high-energy heavy-ion colli- 
sion experiments. This is motivated by the high energy density and temper- 
ature achieved in such collisions, which provide favorable conditions for the 
DCC production scenario based on the non-equilibrium chiral phase transition 
[10], as discussed in Sec. 2.2. Among the various heavy- ion experiments, WA98 
was the first one which had the capability of measuring photon and charged 
particle multiplicities over a common region of phase-space, which is a key 
ingredient for DCC search through neutral fraction fluctuations [20,18]. The 
data analyzed by the WA98 collaboration were taken with the 158 ACeV Pb 
beam of the CERN SPS on a Pb target of 213 fim thickness during a period 
of WA98 operation without magnetic field. The analysis makes use of a subset 
of detectors of the WA98 experiment which were used to measure the mul- 
tiplicities of charged particles and photons. Charged-particle hits (Nch) were 
counted using a circular Silicon Pad Multiplicity Detector (SPMD) located 
32.8 cm downstream from the target. It provided uniform pseudo-rapidity 
coverage in the region 2.35 < r] < 3.75. The photon multiplicity was measured 
using a preshower Photon Multiphcity Detector (PMD) placed 21.5 m down- 
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stream of the target and covering the pseiido-rapidity range 2.9 < rj < 4.2. 
It consisted of an array of 53, 200 plastic scintillator pads placed behind 3Xq 
thick lead converter plates. Clusters of hit pads having total energy deposit 
above a hadron rejection threshold are identified as "photon-hke" . For the 
DCC analysis the pseudo-rapidity region of common coverage of the SPMD 
and PMD was selected (2.9 < 77 < 3.75) by WA98 collaboration. The ac- 
ceptance in terms of transverse momentum extends down to Pt = 30 MeV/c, 
although no explicit pr-selection was applied. As far as analysis goes, there has 
been a global search of DCC formation [18] in the WA98 experiment. Asym- 
metry in A^^ to Nch ratio over the full common phase space of the photon and 
charged particle detectors was used to get some idea regarding DCC forma- 
tion in Pb-Pb collisions at 158 AGeV. Following the theoretical expectations 
[1,2,10] that the isospin fluctuations, caused by formation of DCC, may pro- 
duce clusters of coherent pions in phase space forming localized domains, the 
WA98 collaboration has also searched for localized domains of DCC [20,21]. 
This was done by trying to detect large and localized fluctuations in the ratio 
of photons-to-charged-particles numbers. 



An analysis using the momentum information of charged particles without 
looking at photons was carried out by the NA49 experiment [19] at the CERN 
SPS to look for DCC. There, particles were selected that had a measured track 
length of more than 2 m in one of the two Main Time Projection Chambers 
(MTPC) outside the magnetic field and were also observed in at least one of 
the Vertex TPCs inside the superconducting magnets. The NA49 collabora- 
tion studied particles in the region: 0.005 < Pt < 1-5 GeV/c in transverse 
momentum and 4 < y^r < 5.5 in rapidity. Both the WA98 and NA49 ex- 
periments reported results in terms of upper limit on the frequency of DCC 
formation. 



At the RHIC, the STAR experiment [23] plans to look for DCC formation 
using the combination of photon-multiplicity detector (PMD) and the elec- 
tromagnetic calorimeter (EMCAL) for photon detection along with the time 
projection chamber (TPC) and forward-time projection chamber (FTPC) 
for charged-particle detection [23]. The PHENIX collaboration [22] has also 
demonstrated its capability to study isospin fiuctuations, essentially through 
charged-particle measurements using drift and pad chambers and photon mea- 
surements from the electromagnetic calorimeter of the PHENIX detector [22] . 
Finally, we mention that the ALICE experiment [24] at the upcoming LHC at 
CERN will have the combination of a PMD and a charged-particle multiplicity 
detector (FMD) in the forward rapidity region as well as an electromagnetic 
calorimeter (PHOS) and a time projection chamber (TPC) near mid-rapidity, 
which allow for DCC search in different rapidity regions. 
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3.2 Simulated DCC events 

It is important to model DCC pion production in realistic heavy ion collision 
event generators in order to study the sensitivity of the various techniques 
developed to look for DCC-like fluctTiations in experimental data. This is also 
useful to get a clear idea of the limitations of both the experiments and the 
analysis tools at hand. Finally this provides a basis to interpret actual data, 
as done e.g. by the WA98 [20,18] and NA49 [19] collaborations. Here, we 
present a simple model proposed in Refs. [18,139], where the generic particle 
production is governed by the VENUS event generator [140] and the DCC 
pions are introduced by hand at the freeze-out stage. The goal is to observe the 
effect of DCC domains on measured quantities, assuming DCC pions survive 
till the freeze-out time. As discussed previously in Sec. 2.3, this assumption is 
justified by theoretical calculations [101]. The main features of the model are 
the following: 

• DCC-type fluctuations: 

For DCC simulation in a given domain, the identity of charged pions taken 
pairwise (tt+tt") is changed to neutral ones (7r°7r°) and vice-versa [18,139], 
according to the ideal DCC inverse square- root law, Eq. (10). 

• Domain size: 

The size of a domain is defined in terms of its extent in pseudo-rapidity (r/) 
and azimuthal angle (0). In most analysis described in the following, one 
choscs a domain-size of one unit in rapidity: Ary = 1 and a variable extent 
A(f) in azimuthal angle. This is not a bad assumption as it is expected 
that DCC doamins will occur in small regions and will mostly modify the 
production of low momentum pions. The influence of DCC on the charge- 
neutral pion ratio may then be limited to localized ?7-0 regions due to the 
motion of the DCC domain within the overall collective motion or strong 
radial flow in heavy ion collisions. 

• Number and Pt spectrum of DCC pions: 

For simplicity all the pions inside a chosen domain are considered to be 
DCC-type. As already emphasized, theoretical calculations suggest that 
DCC pions have low px- In order to study the usefulness of pt informa- 
tion of particles, as discussed in Sec. 3.5 below, one typically consider two 
cases [141]: (i) all pions with pt < 150 MeV are taken as DCC pions in 
a given domain; (ii) a variable number (depending on the domain size) of 
DCC-pions with px < 150 MeV are added by hand to the generic sample 
of pions in a given domain. The number of these additional pions depends 
on the size and energy density of the DCC domain. Domains of typical 
radius i? ~ 3 — 4 fm are expected from numerical simulations of the non- 
equilibrium chiral phase transition in the context of heavy-ion collisions 
[96,65,142]. Assuming the energy density eDcc of the DCC domain to be 
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about 50 MeV/fm^, the number of DCC pious is estimated from the re- 
lation |7ri?^(eDcc/'"^7r) to be between 40 and 100. This probably gives a 
reasonable upper bound on actual numbers to be expected. 

• 7r° decay and photon fraction: 

After introducing DCC- type fluctuations the 7r°'s are allowed to decay. The 
neutral pion fraction / then gets approximately modified to: 



where is the photon multiplicity and A^ch is the total multiplicity of 
charged pions. 

• Percentage of events being DCC-type: 

The sets of events with DCC will be referred to as "nDCC" events. The 
ensemble of nDCC events may or may not have all the events of DCC type 
and the percentage a of DCC-type events in a given ensemble can be varied 
to study the sensitivity of various physical and detector related effect to the 
observation of DCC. 

In order to understand the various techniques used for DCC search, one can 
apply them to "realistic" simulated DCC events as described above. For this, 
one assumes [143] a single DCC domain lying within the common pseudo- 
rapidity coverage of two hypothetical multiplicity detectors for photon and 
charged particles. To simulate what happens in a true experimental situation, 
it is necessary to include detector related effects. For photon and charged- 
particle detection the detection efficiencies are taken to be about 70 ± 5% and 
95 ± 2% respectively. It is also known that charged particles sometimes lead 
to photon-hke signals. Such a contamination can be as high as 25% [144] . For 
simulations, one therefore includes a 25% charged particle contamination in 
the photon signal. These parameters are taken considering the WA98 PMD 
[144] and SPMD [145] detectors in view. Finally, we have the size of DCC 
domain in 0, i.e., A0, varied. To study the sensitivity one also prepares several 
sets of data with different DCC fractions a from 1% to 100%, mixing normal 
events and DCC-type events in an appropriate manner. The /-distribution in a 
localized DCC region of Arj = 1 and A0 = 90° within the common acceptance 
of the detectors is shown in Fig. 20. Also shown is the /-distribution of normal 
events. With the model for DCC event simulation in hand, we now discuss on 
analysis techniques developed specifically for DCC search. 
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Fig. 20. Neutral pion fraction (/) distribution from a DCC-domain as incorporated 
in the model for simulated DCC events (see text). The corresponding distribution 
for normal (non-DCC) events is also shown for comparison. Prom Ref. [143]. 

3.3 Experimental techniques developed for DCC search 



In this section we discuss the various techniques developed for looking at 
charged-neutral fluctuations or DCC-type fluctuations. We shall apply some 
of these techniques to simulated DCC events as discussed above and study 
their sensitivity to look for anomalous fluctuations. We follow the several such 
analysis methods proposed in the literature, among which: 

(1) N^. - iVch correlation [18,20,21]; 

(2) discrete wavelet analysis [122,139]; 

(3) robust observables [17]; 

(4) ^-measure [146]; 

(5) event-shape analysis [147]. 

The techniques listed in the flrst three items have been used in actual experi- 
mental analysis, so we describe them in detail in what follows. The technique 
based on event shape analysis is discussed only briefly. In later section 3.4 we 
present the experimental results based on these analysis techniques. 



3.3.1 Photon-to-charged-particle multiplicity correlation 

In order to look for any possible fluctuation in photon and charged particle 
multiplicities, which may have non-statistical origin, it is necessary to look at 
the correlation between iV^-ukc and N^i as a function of resolution in — 
phase-space [18,20,21]. We flrst discuss the actual method of analysis and then 
show the sensitivity of the method. 
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Fig. 21. A^-y-iikc vs. Nch correlation plot for various bins in (f), obtained from normal 
events. Also shown are the correlation line and its parameters. From Ref. [143]. 

In this analysis method event-by-event correlation between A^-y-nke and A^ch 
can be studied in different 0-segments, by dividing the entire 0-space into 2, 
4, 8 and 16 bins of equal size. The correlation plot for each of the 0-segments, 
starting with the case of one single bin - no segmentation - are shown in 
Fig. 21 for a sample of normal events. A common correlation axis {Z) is then 
obtained for the full distribution by fitting the A^-y-uke and A^ch correlation with 
a second order polynomial. The correlation axis with fit parameters are shown 
in the figure. The distances of separation (Dz) between the data points and 
the correlation axis is then calculated with the convention that Dz is positive 
for points below the Z-axis. The distribution of Dz represents the relative 
fluctuations of A^-y-nke and A^ch from the correlation axis at any given (p bin. 

One then applies this technique to simulated DCC events. The DCC events 
considered corresponds to a domain size of 90° in 0. In order to compare the 
fluctuations for different bins in the same level, one introduces a scaled 
variable, Sz — Dz/s{Dz), where s{Dz) represents the root-mean-squared 
(rms) deviation of the Dz distribution for normal events. The presence of 
events with localized fluctuations in A^^_iike and Nch at a given bin is expected 
to result in a broader distribution of Sz compared to those for normal events 
at that particular bin. The typical Sz distributions for both class of events for 
4 bins in cf) are shown in Fig. 22. One clearly observes a difference in width of 
the Sz distribution between the two class of events. So the rms deviation of 
the Sz distribution quantifies the amount of fiuctuation in such an analysis 
method. 

Figure 23 shows the sensitivity of the method for picking up DCC-type fluctu- 
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Fig. 22. Sz distribution for normal events (open circles) and DCC-type of events 
(solid histogram) with A(f) = 90° DCC domain. From Ref. [143]. 



2.25 



2 

c 
o 

1.75 

■> 
-o 

e 1.5 

1.25 



1 10 10" 

% of events being DCC type 

Fig. 23. RMS deviation of Sz distribution for DCC-type of events with 90° DCC 
domain in (/> as a function of percentage of events being DCC-type. This plot shows 
the sensitivity of the method. From Ref. [143]. 

ations for a case of 90° domain size. The width of the Sz distribution has been 
plotted as a function of the number of events being DCC-type. As expected 
the rms deviation of Sz distribution increases with increase in % of DCC-type 
events in the sample. 

This technique has been applied to WA98 experimental charged-particle and 
photon multiplicity data [18,20,21] to look for extra DCC-type fluctuations as 
compared to normal events. 
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3.3.2 Discrete wavelet analysis 



The discrete wavelet technique (DWT), described in Sec. 2.4 above, has the 
beauty of analyzing a distribution of particles at different length scales with 
the ability of finally picking up the right scale at which there is a fiuctuation. 
The method has been shown to be quite powerful to search for DCC-type 
of fiuctuation [122,124,139]. Here we discuss how this method is actually ap- 
plied in analysing experimental data. As before, we analyze the simulated 
normal and DCC events using the multi-resolution DWT technique to search 
for bin-to-bin fiuctuations in the charged-particle and photon multiplicity dis- 
tributions. The analysis has been carried out in Ref. [139] by making 2^ bins 
in where j is the resolution scale. The input to the DWT analysis is a sam- 
ple distribution function at the highest resolution scale jmax- In most analysis 
presented below, the sample function corresponds to the photon fraction: 



where A^-y and N^h are the multiplicities of photons and charged particles re- 
spectively. The 0-space is binned as discussed above and the analysis is carried 
out using the so-called D — A wavelet basis [148]. The output of the DWT anal- 
ysis consists of a set of wavelet or father function coefficients (FFCs) at each 
scale, from j = 0, . . . , jmax — 1- The FFCs at a given scale carry information 
about the degree of fluctuation at higher scales. Due to the completeness and 
orthogonality of the DWT basis, there is no information loss at any scale. For 
the present study the analysis has been carried out by taking jmax = 5. 

Typical FFC distributions at the scale j = 1 are shown in Fig. 24. The dashed 
histogram represents the FFC distribution for normal events. The solid line 
histogram shows the corresponding distribution when a DCC domain of size 
given by A(f> = 90° and 3 < < 4 is introduced in all the events. This 
constitutes an ensemble of DCC events having probabihty o; = 1 of occurrence 
of a DCC. It is observed that the distributions for normal events is Gaussian, 
whereas there is a broadening of the distribution in the presence of DCC. The 
presence of a DCC domain modifies the phase space distributions of charged 
particles and photons, resulting in an increase in the rms deviations of the FFC 
distribution. So in this analysis the rms deviation, s, of the FFC distribution 
quantifies the amount of fiuctuation. 



55 Neglecting detector effects, the photon fraction (182) can be approximately re- 
lated to the neutral pion fraction through (see Eq. (181)): «i 2//(l -|-/). Assum- 
ing ideal DCC production, the corresponding probability distribution is given by: 

dPDCc{fy)/df^ « 1/^/7(2 -/7)^. 

The distribution for normal events is a perfect Gaussian with a = 0.0908 and 
xVndf close to 1 [139]. 
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Fig. 24. FFC distribution for simulated normal events and pure DCC-like events. 
The solid and the dotted lines correspond to the DCC-like and normal events re- 
spectively. From Ref. [147]. 

To quantify the DCC-type effect further, we introduce a strength parameter 
C, given by: 



2 2 
DCC ~ '^normal/' 



^normal 



where Snormai is the FFC rms deviation for an ensemble of normal events and 
•sdcc is the same for an ensemble of events with DCC. Figure 25 shows the 
sensitivity of the method for picking up DCC-type fluctuations for of 
90° domain size. The strength parameter derived from the rms deviation of 
the FFC distributions of normal and DCC type events has been plotted as 
a function of the fraction of events being DCC-type. As expected ( increases 
with the latter. 

This technique has been adopted by the WA98 Collaboration, which measures 
multiplicities of charged particles and photons in overlapping parts of the 
detector phase-space [20,21]. 



3.3.3 Robust observables 

This analysis uses ratios of factorial moments of the charged-neutral pion 
distribution, as described at length in Sec. 2.4. It has been developed for a 
dedicated DCC search at the Fermilab TEVATRON by the MiniMAX Col- 
laboration [17]. The set of robust observables is defined by Eqs. (143)- (144), 
where and are given by the photon and charged-pion multiplicities re- 
spectively in the studied region of phase-space. For an inclusive analysis, the 
n's denote the particle multiplicities in an event and averages are taken over 
a large number of events. For event-by-event analysis, when the multiplicity 
is high, the phase space region is divided into bins and the n's denote the 
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Fig. 25. Variation of as a function of percentage of events being DCC type. DCC 
events have 90° DCC domain in 0. This plot shows the sensitivity of the method. 
From Ref. [143]. 

multiplicities in a bin. In this case, averages are taken over the total number 
of bins. As emphasized previously, the variables r^,! are robust observables in 
the sense that the detection efficiencies, often difficult to estimate - especially 
for photons, mostly cancel out and are thus expected to have little influence 
on the results. The main point [121] is that rj i = 1 for all 2 > 1 in the case 
of generic pion production, whereas rj i = + 1) (sec Eq. (148)) for pure 
ideal DCC production. Application of this technique to simulated normal and 
DCC-type events are shown in Fig. 26. 

This technique has been used by the MiniMAX Collaboration [17] and also to 
analyze cosmic ray experimental data at Pamir. 

3.3.4 ^-measure 

The <l>-mcasure for a system of particles is defined as [149,150]: 



where z — x — x, where x is the value of a given single-particle observable 
in a given event and x its average over all particles and all events. The event 
variable Z is a multi-particle analog of z defined as Z = — x), where 

the summation runs over all particles from a given event. The brackets in 
(184) represent averaging over all events. Note that, by construction, z = 
and (Z) — 0. As described in Refs. [151,152], the ^-measure can be used 




(184) 
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Fig. 26. Results of inclusive DCC analysis using robust observables. Normal events 
(VENUS event generator) do not show any correlation whereas pure DCC events 
follow the relationship given in Eq. (148). The simulated DCC events generated by 
taking detector effects into account follow a similar behavior to that of pure DCC 
events. Prom Ref. [143]. 

to study the fluctuations in particle species. In that case, the single-particle 
variable x = 1 if the particle is of a given sort (say a neutral pion) and x = 
otherwise (charged pion). This has been applied to study fluctuations of the 
neutral pion fraction / from a DCC in Ref. [146]. 

Denoting by N^^ = N^^o + Nch the total pion multiplicity in a given event, one 
can write A^^o = JNt^ and A^^± = (1 — f)N.^, where / is the fraction of neutral 
pion in the event. Assuming that the relative fluctuations in A'^ are small, 
that is {6N^) < (A/'^)^ one obtains [146]: 



^c^^{K){5P)-^{f){l-f). (185) 

It is instructive to consider the properties of the ^-measure in the following 
simple cases: 

• Non-DCC case: 

It is known from pp experiments [153] that the produced pions have their 
charge states partitioned binomially with (/) = 1/3. The fluctuation in / 
is inversely proportional to the total number of pions and is approximately 
given by {5p) ~ — f)/{Nn)- In that case, one obtains, for the 

measure: 

'^'non-DCC, uncorr. = 0. (186) 
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DCC pion fraction (P) 



Fig. 27. Variation of $ — ^normal as a function of the fraction (3 of DCC pions, 
obtained from simulated data using VENUS. Prom Ref. [146]. 

• No DCC, but Nt^o and N^h have non-trivial correlations: 

If A'^o and A^^± are assumed to be correlated in such a way that there are no 
(DCC-type) event-by-event fluctuations of the neutral fraction /. Moreover, 
the latter is assumed to be strictly independent of the event multiplicity. 
Then, f = S, where 5 is a constant smaller than unity, and (Sp) — 0. One 
has: 



Notice that The ^-measure is negative in that case. 
• Ideal DCC case: 

For ideal DCC production, the probabihty distribution of the neutral frac- 
tion is given by the inverse square-root law, Eq. (10). One has (/) — 1/3 
and {5p) — 4:/ 4:5, so the ^-measure is given as: 



It is positive for (A^^r) > 5/2. For a typical case where the total pion multi- 
plicity ~ 300, one finds $dcc ~ 4.7. It is worth emphasizing that this 
assumes that all the observed pions are of DCC origin. A more realistic case 
would include pions from non-DCC sources as well [143]. This, along with 
the decay of neutral pions to photons as well as detector effects, have been 
discussed in Ref. [146]. 




(187) 




(188) 
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DCC event fraction (a) 

Fig. 28. Variation of $ — ^normal ^ a function of the DCC event fraction (a) for a 
fixed DCC pion fraction (/? = 0.25), obtained from simulated data using VENUS. 
From Ref. [146]. 



The $- measure can be appHed to simulated data (both normal and DCC- type) 
in order to study the sensitivity of the measure to detection of DCC-type 
fluctuations. First wc study the effect of varying the fraction f3 of DCC pions 
in each event. To this aim, one calculates the quantity $ — $normab where $ 
denotes the measure of fluctuation for a given DCC pion fraction /3 and $normai 
is the measure of fluctuation for normal simulated events (detector and decay 
effects included). The variation of $ — $normai as a function of f3 is shown in 
Fig. 27. One observes that the size of anomalous ffuctuations decreases with 
decreasing DCC pion fraction, which is as per expectation. The statistical error 
on $ was calculated to be 0.006, from which one finds that such observable 
is a sensible probe of possible DCC-type fiuctuations if the fraction of DCC 
pions in a typical event is above ~ 3%. 

This analysis assumes that all events have DCC pions. Including the effect 
of the probability of DCC formation being a ^ 1 has been studied for a 
fixed fraction of DCC pions in each event in Ref. [141]. There, the authors 
assume that for DCC-type events the percentage of DCC pions is P — 0.25 
and vary the fraction a of events being DCC-type. The results are presented 
as $ — $normai VS. a iu Fig. 28. As expected, the value of $ — $normai decreases 
with decreasing a. Keeping in mind the above mentioned statistical error on $- 
measure, we find that for (3 = 0.25, the measure is sensitive to DCC formation 
if the latter occurs in more than 1% of the events. 
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3.3.5 Event shape analysis 



This method has been proposed in the context of DCC search in Ref. [147]. 
It is based on a technique that has been used successfully for flow analysis 
of heavy-ion data [154]. It uses the fact that localized DCC formation should 
lead to an event-shape anisotropy, which is expected to be out of phase for one 
detector (e.g. charged pions) compared to the other (e.g. photons). In other 
words, whenever there is large number of charged particles from a DCC region 
recorded in a given region of the first detector, there should be a depletion 
in the number of neutral particles recorded in the corresponding zone of the 
other detector. Therefore, from an event shape analysis using both detectors 
one can, in principle, look for DCC signals. In terms of the terminology used 
in flow analysis, this method is based on the fact that a simple redistribution 
of particles, with two detectors to detect charged particles and photons, would 
result in the same flow direction with the flow angle difference peaking at zero. 
However, events where the neutral pion fraction has been modified according 
to the ideal DCC probability distribution, will show the flow angles in two 
detectors to be almost 90° apart [147]. 

In Ref. [147], events with localized DCC domains were simulated along the 
lines described earlier in Sec. 3.2 and event anisotropy (indicating flow) was 
introduced according to the following simple procedure. First, a flow direc- 
tion is selected at random, uniformly distributed between 0° and 360°. The 
directions of particles are then distributed about this direction following a 
Gaussian distribution with a width of 10°. Results were analyzed using the 
standard second order Fourier analysis of elliptic flow [154]. The particle hits 
in a given detectors can be labeled by their positions in the rj — 4> plane. For 
each event, one constructs the two sums: 



X = ^cos(20i) (189) 

i 

y = ^sin(2(/),), (190) 

i 

where 0j is the azimuthal angle of the particle i and the sums run over all 
particles in a given event. The flow angle S is defined as: 

S = ^ arctan(y/X) . (191) 

In the absence of any detector imperfections and other geometrical effects, 
the distribution of S, taken over a large number of events, is expected to 
be completely flat over 0° to 180°. This is because the flow direction varies 
randomly from event to event. However, when the events are realigned, with 
respect to the flow angle in each event, as described above, one can see the 



82 



characteristic peaks (at 0° and 180°) in the azimuthal distribution of particles. 

In case of two detectors with the same phase space (77 — (f)) coverage, one 
detecting photons and other detecting charged particles, the situation is very 
interesting. If there is genuine flow in a particular event both detectors would 
show the effect in terms of their respective S-angles getting ahgned in the 
same direction. Therefore the distribution of ^ = S^^, — Sch, the difference 
between the flow angles of the two detectors, is expected to be peaking at 
zero. However, in case of DCC being prominent in a particular region, there 
will be more photons detected in one detector. The other detector is expected 
to show less charged particles in the same region of phase space. Therefore an 
event shape analysis is expected to show two flow angles for both detectors 
which will be out of phase (for the present second order Fourier analysis the 
angular difference is expected to peak at 90°). 

The results of Ref. [147] are presented in Fig. 29. The three left panels (a)- 
(c) correspond to the case with only DCC in which one can notice (panel 
(a)) a clear peak at 90° for the angle ^ between the event planes for the 
two detectors. Panels (b) and (c) show the angle between the event planes as 
obtained for two sub-events in each of the detectors separately. It is important 
to notice that both detectors show flow, but one with respect to the other 
clearly shows an anti-flow type behavior. In the same flgure is also presented 
the data corresponding to pure VENUS events for comparison. One notices 
that the latter do not present neither any signature of flow nor DCC-like 
fluctuation. Finally, results shown in panels (d)-(f) show the same plots for 
flow type events which have no DCC-like fluctuation. Here, the individual 
detectors are seen to show the same effect as one with respect to the other. 

A more detailed analysis Ref. [147] reveals that this method is sensitive if at 
least 10% of the events are of DCC-type. It has been argued that one can 
improve upon this by combining with the DWT technique Ref. [147] . The 
event-shape analysis has not yet been applied to experimental data. 



3.4 Experimental results 



We have discussed in the above section, the experimental set-up for various 
experiments looking for DCC and the experimental techniques used to look 
for DCC-type fluctuations in experimental data. Below we present the results 
from various experiments which have looked for possible DCC formation in 
high energy collisions. 
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Fig. 29. The distribution of for simulated DCC (a)-(c) and flow (d)-(f) 

events. Panel (a) shows anti-correlation between event planes determined for 
charged-particle and photon detectors. Panels (b) and (c) show correlations be- 
tween event planes obtained considering two sub-events for the photon and 
charged-particle detectors respectively. The dotted lines correspond to generic 
VENUS events. In case of flowy events, panels (d) shows correlation between event 
planes obtained from the two detectors, whereas panels (e) and (f) are analog to 
(b) and (c). From Ref. [147]. 

3.4-1 Results from cosmic ray experiments 

There have been many cosmic rays experiments [6,7]. Here we discuss the 
most recent results from a systematic study of the large asymmetries in the 
neutral pion fraction distribution in high-energy cosmic-ray families (100 TeV 
< Eyisibie < 700 TeV) detected at high mountain altitudes at Pamir [7]. The 
experimental data has been analyzed using the technique of robust observables 
discussed earlier. The results for the robust ratios rj^i is presented as a function 
of leading jet energy Eiead jet in Fig- 30. The pattern of energy flow in cosmic- 
ray interactions is understood through the picture of jets: The stream of energy 
in the air cascade are represented by jets produced in the nuclear interactions. 
In this picture the leading jet represents the pattern of atmospheric nuclear 
and electromagnetic cascade having the largest energy flow. 

Prom the figures one observes that the absolute values of r,^i are decreasing 
with increasing energy of the jets and the order i of the moments. The results 
from data have been compared to simulated generic events. One observes that 
with the increasing order of the robust observables the difference between 
experiment and simulation becomes more important. Recall that values of rj^i 
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Fig. 30. The robust observable ri^i, r2,i and r^^i as a function of visible leading jet 
energy: experimental points are represented by triangles, whereas squares are used 
for simulated data. Prom Ref. [7]. 



far below 1 can be an indication of events with DCC formation overlaying 
generic events. The experimental average values of r3 i that are most sensitive 
to DCC formation are clearly below 1 but by less than a standard deviation. 
The averaged experimental values for rs^i are presented in Table 2, together 
with the 95% confidence interval. The latter does not exceed 1. 



Table 2 

Experimental statistics and 95% confidence interval of the average rs^i in Pamir 
experiment. From Ref. [7]. 



Jet energy (TeV) 




Number of events 


95% confidence interval 


50 


0.48 ± 0.44 


78 


(0.38,0.58) 


100 


0.45 ± 0.41 


61 


(0.38,0.55) 


200 


0.43 ± 0.40 


30 


(0.28,0.58) 


300 


0.47 ± 0.49 


10 


(0.12,0.82) 



The conclusion from Pamir experiment [7] is that, due to the limited large error 
bars, the experimental data neither confirm nor rule out the possibility of a 
DCC formation mechanism in high-energy cosmic-ray interactions. A further 
increase in experimental statistics is needed to establish firm conclusions. 
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Fig. 31. Number of events plotted against multiplicities of charged particles observed 
in the pseudo-rapidities interval — 2 < 77 < 2. The distributions are drawn separately 
for each value of observed photon multiplicity in the same range of pseudo-rapidity. 
The region which is populated by Ccntauro-likc events if produced at the collider 
and observed by the detector is indicated. Prom Ref. [14]. 

3.4-2 Results from nucleon-nucleon collision experiments 

Here we discuss results on DCC search from high-energy nucleon-nucleon ex- 
periments UAl and UA5 at CERN [13,14,15] and MiniMAX at Fermilab [17]. 



3.4.2.1 The UA5 Experiment Figure 31 shows the distribution of events 
as a function of observed charged-hadron multiplicity in the pseudo-rapidity 
interval — 2 < 77 < 2 for each value of observed photon multiplicity in the 
same interval, for pp collisions at center of mass energy of 540 GeV [14]. This 
correlation plot does not support the existence of Centauro type of events, 
based on the characteristics of Ccntauros observed in cosmic ray interactions. 
One observes that there is no event with charged hadron multiplicity remotely 
comparable to the expected value of 30 — 40 with associated photon mul- 
tiplicity being close to zero. The shaded portion represents the region where 
one would expect some events in pp collisions at y/s — 540 GeV, had they 
similar characteristics as Centauro events in cosmic-ray experiments. 

The UA5 Collaboration [15] also looked for Centauro-hke events in the 900 GeV 
center of mass energy data, by studying the correlation of produced charged 
particles A'^ch versus photons (A^) and correlation between the electromagnetic 
(£^em) and hadronic energies (£^had) deposited in the calorimeters. Figure 32 



^' This number has been arrived at by making a comparative study of the collider 
data with the data obtained from the cosmic ray experiments [14] . 
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Fig. 32. The asymmetry in electromagnetic and hadronic energy (R) is plotted 
along the Y-axis and the transverse energy is plotted along the X-axis. Comparison 
to various model calculations are shown. From Ref. [15] 



shows the average value of the ratio 



R 



had 



Eha.d + Eg 



(192) 



as a function of the transverse energy E^, for experimental data, simulated 
normal and Centauro-like Monte-Carlo events, as well as a mixture of normal 
and Centauro events, the percentage of Centauro events being 0.24% of all 
events. The data clearly seems to follow the trend of normal events. Further 
the correlation of Nch with respect to A^^ for data and Monte-Carlo simulation 
look similar, thereby ruling out the possibility of Centauro events or DCC-type 
events in the experimental data. 



3.4.2.2 The UAl Experiment The UAl Collaboration [13] have looked 
for Centauro-type events in pp collisions at center of mass energy of 540 GeV 
in 48, 000 low bias events. They have measured the correlation between de- 
posited transverse electromagnetic and hadronic energies, where the energy 
deposited in the first four radiation lengths of material is termed as electro- 
magnetic energy and that deposited beyond twelve radiation lengths is termed 
as hadronic energy. By comparing experimental data with Monte-Carlo simu- 
lations, they have excluded possible Centauros or global DCC-type events in 
their data sample. 
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Tabic 3 

Values of rj^i for the whole event sample (about 13.8 x 10^ events) of the MiniMAX 
experiment [17]. 



i 


ri,i 


1 


1.026 ± 0.004 


2 


1.035 ±0.010 


3 


1.059 ± 0.027 


4 


1.118 ±0.065 


5 


1.310 ±0.151 


6 


1.904 ± 0.382 



3.4.2.3 The MINIMAX Experiment The analysis of the pp colhsions 
data at center of mass energy of 1.8 TeV for the MiniMAX experiment [17], 
in the forward direction 3.2 < r) < 4.2, was carried out using the technique of 

robust observablcs. which, as already emphasized, had been actually developed 
for this experiment. A total of 1.3 x 10^ events was analyzed. A detailed 
Monte-Carlo simulation specific to the MiniMAX experiment, with input being 
minimum bias PYTHIA events passed through full detector simulation and 
analysis chain, yielded a value of ri,i ~ 1.02 — 1.13 for the lowest robust ratio. 
A similar simulation with pure ideal DCC introduced by hand yielded a value 
of Ti l ~ 0.6 — 0.7. This demonstrated the sensitivity of the experiment to 
detect a single large ideal DCC. The analysis of the data reveals that the 
ratios rj^i are not smaller than one as would be expected for pion production 
from a DCC. The measured values of the ratios r^^i are given in Table 3. 

The MiniMAX Collaboration was able to put some limits on DCC production 
in terms of the size of the DCC - that is the number of DCC pions produced 
- and the hkelihood of DCC production. To this aim, two types of model for 
DCC production have been considered: 



(i) Associated production, where the number of produced DCC pions in each 
event is proportional to the number of generically produced pions; 

(ii) Exclusive production, where a given event is described by either pure DCC 
or pure generic production with an given probability. 

Knowing the most probable lower limit of the ratios Vi^i in data, they were 
able to extract an upper limit on DCC production for the above two scenarios. 
For the associated production model, the upper limit on DCC pion production 
has been set at 

{N)dcc_ < Q_21 . (193) 



(A) 



generic 
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Fig. 33. This is the scatter plot showing the correlation between N^h and A?^_iike- 
The solid outline shows the trend of the minimum bias data. The central sample 
(with Et > 300 GeV) is shown as points for the data and as a hatched region for 
VENUS (with much lower statistics). Overlaid on the plot are the Z axis and the 
Dz axis at a particular value of Z as explained in the text. Prom Ref. [18]. 

For the exclusive production model, the upper limit on DCC production is set 
in terms of probability that an event is DCC type. This value must be less 
than 0.05. 

3.4-3 Results from nucleus-nucleus collision experiments 

Here we discuss results on DCC searches from the high-energy heavy-ion colli- 
sion experiments at the CERN SPS (WA98 [20] and NA49 [19] Collaborations) 
and at the RHIC (PHENIX Collaboration [22]). 

3.4.3.1 The WA98 Experiment DCC search in WA98 experiment in- 
volved three kinds of analysis. One was to look for global DCCs [18], that 
is single large DCC domains, by measuring photon versus charged-particle 
correlation in the full common r] — (j) phase-space of the photon (PMD) and 
charged-particle (SPMD) detectors. Second was to look for smaller domains 
of DCC in common localized regions oi t] — (p phase space of the PMD and 
SPMD [20,21]. And the third was to look for anti-Centauro-type events in the 
data by looking at localized regions in the event displays [155], as for example 
in the JACEE cosmic-ray experiment. Below we discuss the results on DCC 
search from the above three analysis carried out by the WA98 Collaboration. 

- Global DCC search: 
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Fig. 34. Sz distribution for the experimental data is shown, overlaid with VENUS 
simulations incorporating 0%, 25% and 60% DCC-typc events. All of the distribu- 
tions are normalized to the total number of data events. Prom Ref. [18]. 

The WA98 Experiment attempted to look for DCC-type fluctuations by look- 
ing at event-by-evcnt correlation in the number of charged particles and pho- 
tons detected over the full common — coverage of the PMD and the SPMD 
[18]. The strong correlation between charged and neutral multiplicities shown 
in Fig. 33 suggests to employ a coordinate system with one axis being the 
measured correlation axis and the other being perpendicular to it. If all de- 
tected particles were pions and the detectors were perfect and had identical 
pseudo-rapidity acceptance, then the correlation axis would be a straight line. 
The details of the correlation analysis has already been described in Sec. 3.3.1. 
Here we present the results of such analysis apphed to the WA98 data. 

The Sz distribution for both data and simulated events are shown in Fig. 34. 
The simulated events have been generated using the VENUS 4.12 event gen- 
erator, with output processed through a detector simulation package in the 
GEANT 3.21 framework [156]. This incorporates the full WA98 experimental 
setup. The centrality selection for simulated data has been made in an iden- 
tical manner as in the data, determined from the simulated total transverse 
energy of mid-rapidity calorimeter used in the experiment. The discrepancy 
between VENUS and the data can be seen more clearly by measuring the 
width of the Sz distribution. The simulated and measured distributions are 
both Gaussian with respective widths crgimui. — 0.998 ± 0.002 (fit error only) 
and (Tmeas. = 1-13 ± 0.07 (error from relative scale uncertainties included). 

One expects that DCC events would show up as non-statistical tails on the Sz 
axis. No such events were seen by the WA98 Collaboration in their global anal- 
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Fig. 35. The 90% C.L. upper limit on DCC production per central event as a 
function of the fraction of DCC pions under two assumptions. The thick line gives 
the upper limit obtained by assuming the a^z is completely given by the 

VENUS calculation requiring to make a cut at 6ct. The dashed line shows a less 
conservative limit obtained by using the ctdz measured in the data itself. This 
allows one to make a tighter cut at 5fT, increasing the DCC detection efficiency. 
From Ref. [18]. 

ysis of the data sample. They concluded that a possible single DCC domain is 
either very rare, very small, or both. To check which hypothese is consistent 
with their data, they determined upper hmits on the frequency of DCC produc- 
tion as a function of its size, as represented by the fraction (3 = C, = Ny^qq/Nt, 
of DCC pions. For this they have computed Sz distributions for several values 
of C, ranging from 15% to 90%, and they have considered two possible scenar- 
ios. The first is based upon the conservative assumption that VENUS should 
describe the data perfectly in the absence of a DCC signal. The corresponding 
90% C.L. excluded region for possible DCC formation is presented in Fig. 35 
as a solid line. The second scenario assumes that the difference between the 
data and VENUS seen in Fig. 34 is due to detector effects and that the widths 
should in fact be the same. The resulting excluded region is shown in Fig. 34 as 
a dashed line. The boarders of the two regions are quite different at C = 15%, 
but get closer at ^ > 30%. 

- Search for localized domains of DCC: 

Photon and charged-particle multiplicity correlation studies have also been 
carried out in localized regions of the rj ~ (p phase-space to look for the pos- 
sible existence of small localized DCC domains. This analysis uses both the 
techniques of — N^h correlation and DWT in various bins in azimuthal 
angle [20,21], discussed previously under sections 3.3.1 and 3.3.2. The WA98 
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Collaboration have also studied the centrality dependence of this analysis. 
For this, four different centrality bins were considered: top 5% (henceforth 
referred to as centrality-1); 5% — 10% (centrality-2); 15% — 30% (centrality-3); 
and 45% — 55% (centrahty-4) of the minimum bias cross-section. 

In looking for non-statistical DCC-type fluctuations in experimental data, it 
is necessary to understand all detector related effects and to try to filter out 
the various components contributing to the total observed fluctuations. For 
this, results from experimental data need to be compared to some base line 
results to draw proper conclusions. Due to the inherent uncertainties in the 
description of "normal" physics and detector response in simulations, the ob- 
servation of an experimental result which differs from the case with statistical 
fluctuation, given by simulations, cannot be taken alone as evidence of pres- 
ence of non-statistical fluctuations. This problem can be solved by generating 
different sets of mixed events from data, keeping specific physics goals and 
detector effects in mind. This remarkable ability of mixed events to probe dif- 
ferent physics goals is possible through construction of suitable ensembles of 
mixed events by removing various correlations in a controlled manner while 
preserving the characteristics of the measured distributions as accurately as 
possible. Fluctuations in the ratio of to A^ch can have three different ori- 
gins: They can arise due to independent fluctuation in either only, A'ch 
only, or both. Further, in the latter case the fluctuations in A^-.^ and Nch may 
be correlated event by event. In particular, such a correlation is expected to 
arise in the case of DCC production. By intelligently constructing different 
types of mixed events one can probe possible non-trivial fluctuations in the 
N^-toNch ratio and their origin. 

A first set of mixed events, referred to as Ml-type, is constructed to provide 
a base line for comparison to real event. These are generated by mixing hits 
in both photon (PMD) and charged-particle (SPMD) detectors separately, 
satisfying the A^^-uke — -^ch correlation as in the real event. This is done by 
producing one mixed event for each real event with the same multiplicity of 
A'^-iikc and A'ch pair as in the VENUS event. The idea is to globally keep the 
distributions similar to real events and look for event-by-event fiuctuations in 
localized regions of phase-space. In the construction of such mixed events, care 
is taken such that no two hits come from the same real event. Moreover, hits 
within a detector in the mixed events are not allowed to lie within the two 
track resolution of the detector. Results from such mixed events when com- 
pared to those of real events were expected to reveal the presence of localized 
fiuctuations. However, they do not give any information neither regarding the 
origin of such fiuctuations nor concerning possible — N^h correlations. 

A second set of mixed events, referred to as M2-type events, has been con- 
structed to provide information regarding the possible presence of localized 
event-by-event correlated fiuctuation in A^-y-uke and A^ch- They are generated 
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Tabic 4 

Sources of fluctuations preserved by various mixed events. 



Type of 




Mixed 


Event 




fluctuation 


Ml 


M2 


M3-7 


M3-ch 


A'^^ only 


No 


Yes 


No 


Yes 


Nch only 


No 


Yes 


Yes 


No 


correlated Nj — Nch 


No 


No 


No 


No 



by mixing photon events with difTerent charged-particle events so that no two 
events were repeated. Here also, care is taken that the original A^^_iike — -^ch 
correlation is maintained. The original event-by-event hit structure in each 
individual detector being left unchanged, such mixed event set keeps identi- 
cal the individual localized fluctuations due to iV^_iike and A'ch, but removes 
possible event-by-event correlations. Therefore, by comparing the results from 
such mixed events to those from real events reveals the presence of localized 
event-by-event correlated fluctuations in A^^-uke and A^ch- 

Finally, two other sets of mixed events, referred to as M3-7 and M3-ch type, 
have also been considered to provide information regarding contribution to 
localized fluctuations in the Ny-to-Nct ratio from either A^^ fluctuations or 
A'ch fluctuations separately. Such events are generated from real events by 
mixing hits in one of the detectors (following the procedure of construction 
of Ml-type events) and keeping the hit structure of the event in the other 
detector intact. M3-7 events corresponds to the case where the hits within the 
photon detector are unaltered while the charged-particle hits from different 
events are mixed, whereas M3-ch events correspond to the opposite case. The 
total number of mixed events generated is the same as the number of real 
events. Results from such mixed events, when compared to those of real and 
Ml-type events, are expected to reveal the presence of event-by-event localized 
fluctuations in A'^-uke and A^ch individually. The various sources of fluctuation 
in the A^^-to-A'ch ratio which are preserved by the different types of mixed 
events described above are summarized in Table 4. 

Mixed events constructed from experimental data were compared to simulated 
events, generated as for the global DCC case discussed previously, and were 
analyzed using both the correlation and the DWT analysis techniques. Since 
the widths (rms deviations) of the Sz and FFC distributions carry the infor- 
mation concerning fluctuations, below we compare these widths as obtained 
from data as well as mixed and simulated events. 

The rms deviations of the Sz distributions calculated at different bins are 
shown in Fig. 36 for data, mixed (Ml, M2, M3-7 and M3-ch) and VENUS 
events, in various centrality bins. The statistical errors on the values are small 
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Fig. 36. RMS deviations of the Sz distributions for data, mixed and simulated 
events for the four centrahty bins: top-left for centrality-1, top-right for centrality-2, 
bottom-left for centrality-3 and bottom-right for centrality-4. Prom Ref. [21]. 

and are within the sizes of the symbol. The bars represent statistical and 
systematic errors added in quadrature. The total error is not more than 3% for 
the most central events and goes to not more than 10% for the most peripheral 
events. From Fig. 36, one observes that the nature of the rms deviations of 
the mixed events closely follow those of the data. The mixed events have 
been constructed such that the A^-y_iike vs. A^ch correlations are maintained 
for full azimuth - that is one single bin, so the rms deviations of data and 
mixed events are identical by construction in that case. The rms deviations of 
M2-type events are found to agree with those of the experimental data within 
errors for all four centrality classes, thereby suggesting the absence of localized 
correlations in A^-^-uke and N^h- The same for Ml-type events are found to be 
lower than those obtained for data for 2, 4 and 8 bins in (j) for centrahty bins 1, 
2 and 3. The results form M3-type events are found in between those obtained 
from data and Ml-type events. These results indicate the presence of localized 
fluctuations in data, due to both photons and charged particles. However for 
the case of centrahty-4, the rms deviations of -S'^ distributions from data and 
the mixed events are found to more or less agree with each other within the 
quoted errors. The rms deviations for the simulated events are 1 for all bins 
in azimuth by definition of Sz- 

The rms deviations of the FFC distributions calculated at different bins are 



94 



0.13 



0.12 



E 

O 0.11 



0.1 



0.15 



0.2 r 



S "'IS 

■o 



0.16 - 



(b) 



Data 
MixM2 
Mix M3-ch 
Mix M3-Y 
Mix Ml 
Venus+Geant 



4 8 16 32 

Number of bins in ({> 

• Data 

□ MixM2 



Mix M3-ch 
Mix M3-y 
Mix Ml 

Venus+Geant 



4 8 16 32 

Number of bins 



0.13 

e 

u 

u 

E 0.12 



u 
ft 



0.24 



• Data 

□ MixM2 

T MixM3-ch 

A MixM3-Y 

A Mix Ml 

Venus+Geant 



4 8 16 32 

Number of bins \n& 





• 


Data 




□ 


Mix M2 




T 


Mix M3-ch 




A 


Mix M3-Y 






Mix Ml 







Venus+Geant 












1 




1 



Nurnberofbinsin* 



32 



Fig. 37. RMS deviations of the FFC distributions for data, mixed and simulated 
events for the four centrahty bins: top-left for centrality-1, top-right for centrality-2, 
bottom-left for centrality-3 and bottom-right for centrality-4. From Ref. [21]. 



shown in Fig. 37 for data, all four types of mixed events, as well as generic 
VENUS events. We observe that the rms deviations of data, VENUS and 
mixed events match with each other (within quoted errors) for 32 bins in for 
all the the four centrality classes. The VENUS results are close to the Ml-type 
events, for centrality classes 1, 2 and 3 and are slightly higher for centrahty-4. 
The rms deviations of the FFC distributions for M2-type events are found to 
agree closely with those for data for all centrality classes and all bins in (p. The 
rms deviations for M3-type events lie in between those obtained from data and 
Ml-type mixed events. These results are consistent with those obtained from 
the Sz distributions. They indicate the presence of individual localized fluctu- 
ations in photons and charged particles for certain bins in azimuth (resulting 
in fluctuations of the ratio, as seen). However, they show no clear evidence for 
localized correlated fluctuations (DCC-type) between iV^-uke and N^h- 

Comparing these results with those obtained from simulated DCC events de- 
scribed in Sec. 3.2, one can extract an upper limit on the occurrence probability 
of DCC-like fluctuations at the 90% confldence level, following the standard 
procedure [157] as described below. This has been done, in particular, for cen- 
tral Pb-Pb events with Et > 341 GeV at 158 AGeV collision energy, as DCC 
formation is expected to be more likely in these events. First, a quantity x is 
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Fig. 38. 90% confidence level upper limit on DCC production per central event at 

SPS as a function of DCC domain size in azimuthal angle within the acceptance 
of the detectors. The solid line is for top 5% central events and dashed line for 
5% - 10% centrality class. From Ref. [21]. 

defined as: 



where s and Si correspond to the rms deviations of the FFC (or Sz) distribu- 
tions of M2-type mixed events and real data respectively. It is assumed that 
errors are Gaussian-distributed. Since errors are asymmetric one considers the 
largest error for the calculation to be on the conservative side. The errors on 
the rms deviations are propagated to get errors on x- Upper limit contours 
with 90% CL. have been calculated by taking the limit at % + 1.28e^ where 
is the error in % from Sz and FFC analysis. The two upper limit contours 
set by the two different analysis were found to be consistent with each other. 
It must be mentioned that, for calculating the upper limits, it has been as- 
sumed that the total difference in rms values of data and M2-type events is 
due to DCC-like fluctuations. Now, to relate them to DCC domain size and 
frequency of occurrence of DCC, one proceeds as follows: Obtain the rms de- 
viations of FFC distributions for various domain sizes (from 15° to 180° in 
(f) with 15° increments) for different values of the frequency of occurrence of 
DCC (varying from 0% to 100%), where simulated DCC events are generated 
from the simple model discussed earlier. Then carry out the entire analysis 
as done for data to get x values for each set of events. One then finds for 
which frequency of occurrence for a fixed DCC domain size, the x value from 
simulation matches with that of the x+ l-28e^ for data. This is done for each 
of the above domain sizes. This is then used to set the upper limits in terms 
of domain size and frequency of occurrence of DCC, as shown on Fig. 38. It 
may be mentioned that results from both the Sz and FFC analysis methods 
gives similar upper limits. 
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Fig. 39. Plot showing photon hits (PMD) and charged particle hits (SPMD) in an 
azimuthal plane. The marked 90° patch corresponds to fmax = 0.77. Prom Ref. [155]. 

Thus, within the context of this simple DCC model, an upper limit on the 
presence of localized non-statistical DCC-hke fluctuations is found to be of 
the order of 10~^ for an azimuthal domain size 45° < < 90° and aboTit 
3 X 10~^ for 90° < A(f) < 135°. This upper limit is the most strict limit set 
so far at SPS. Compared to earlier limits set in global analysis of WA98 and 
NA49, this analysis is more sensitive to possible presence of small localized 
domains of DCC in rj — (j) phase-space. 

-Search for anti-Centauro events: Finally, a new method to look for anti- 
Centauro events in WA98 data has been proposed [155], where one studies 
events having photon excess in azimuthal patches within the overlap zone of 
the PMD and the SPMD. One of the advantage of searching for anti-Centauro 
type events is that the purity of the photon sample observed by the PMD will 
be higher than those for normal patches due to the depleted fiux of charged 
particles. The WA98 Collaboration have analyzed 196K events corresponding 
to the top 15% of the minimum bias cross-section. One looks for fluctuations 
in the neutral pion fraction in locahzed regions of the rj — (/) phase-space on 
an event-by-event basis. A particular azimuthal window is selected in the 
pseudo-rapidity range 2.9 < rj < 3.75 and the entire azimuthal range from 0° 
to 360° is scanned in order to find the patch having the maximum value of /, 
referred to as fmax- This scan is performed by successive 2° rotations of the 
A77 — A0 patch in the r) — (j) plane. To minimize the statistical fluctuations, 
a patch with maximum / value in an event was required to have at least 40 
photons corresponding to a 15% statistical error. 

Results from such an analysis have been compared with simulated events 



97 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Neutral pion fraction (f) 



Fig. 40. fmax distributions for 60° patch for data (solid histogram), mixed Events 
(dashed histogram) and V+G (dotted histogram). Prom Ref. [155]. 

obtained using the VENUS 4.12 event generator and processed through the 
WA98 detector setup using GEANT 3.21. These are referred to as V+G events. 
The statistical significance of these results is obtained by comparing them 
with those obtained from a similar analysis applied to the various types of 
mixed events discussed above. An event displaying PMD hits (filled circles) 
and SPMD hits (open squares) within the overlap 7] — (f) zone (obtained by 
using the method described above) is shown in Fig. 39. The patch of size 
A(j) = 90° in azimuth having the highest neutral ratio / is also marked. The 
number of charged particles in this patch is only 12 as compared to 84 pho- 
tons, which corresponds to fmax — 0.77. Figure 40 shows histograms of the 
maximum value /max of the neutral fraction in each event, for data, mixed 
events, and V+G events. This corresponds to a patch size of A(f) = 60°. It 
is seen that the fmax distribution for the data extends to much larger values 
than for mixed and V+G events. 

To test the authenticity of these events, the neutral pion fraction / has also 
been calculated for the immediate preceding and succeeding events in the 
data in the same patch in r] — (f). The distribution of these / values is shown 
in Fig. 41. The solid line corresponds to the neutral pion fraction for patches 
having fmax > 0.55 with > 40 and the dashed histogram represents the / 
distribution for preceding and succeeding events in the data sample. The latter 
has its peak at about 0.35, which is characteristic of generic (normal) events. 
This shows that events with fmax > 0.55 are indeed special events, having large 
non- statistical charged-to-neutral fluctuations. In this sense , these events [155] 
may resemble the anti-Centauro events found in cosmic ray experiments [6]. 
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Fig. 41. Neutral pion fraction (/) distributions for events having fmax > 0.55 (solid 
line) and for preceding and succeeding events (dashed line) for 60° patch having 
> 40. Prom Ref. [155]. 

In summary, using this new event-by-cvcnt analysis the WA98 have found 
events with large charged-to-neutral fluctuations in their data, which exhibit 
fmax > 0.55 in certain A0 domains. These domains are uniformly distributed 
in azimuth. The fraction of such events are much larger in the data as com- 
pared to those seen in the mixed events and V+G events. The percentage of 
these events increases significantly with decreasing patch size. Also the frac- 
tion of such events increases significantly as one decreases the centrality. The 
special events selected by this procedure appear anti-Centauro like [155] simi- 
lar to those found in cosmic ray experiments [6] . A more careful analysis of the 
statistical significance of these events is needed to conclude about observation 
of DCC. 



3.4.3.2 The NA49 Experiment The NA49 Experiment, which detects 
only charged hadrons, has tried to put an upper limit on DCC production 
at SPS from the measurement of fluctuations in mean transverse momen- 
tum M{pt) [19]. They argue that fluctuations in M{pt) are also relevant for 
models of processes that lead to non-statistical fluctuations localized in trans- 
verse momentum, like the formation of DCCs. As we have seen previously 
in the theory review, pions emitted from DCC domains are expected to be 
preferentially produced at low transverse momenta. This provides a transla- 
tion of the number-fluctuations, expected from dynamical DCC simulations 
(see Sec. 2.4), into Pt fluctuations accessible to the NA49 experiment. For 
comparison purposes NA49 used the same type of simulated DCC events as 
the WA98 Collaboration [18], which we described in detail above. The DCC 
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Fig. 42. Limit on the amplitude of fluctuations in the px parent distribution as a 
function of the frequency of occurrence of the event showing the fluctuation. From 
Ref. [19]. 

production is characterized by the probability a = F to form a single DCC 
domain in an event and the fraction /? = ^ of pions coming from the DCC. 
They made the additional assumption that the DCC pions are produced with 
Pt < Pt'^^ = ^-K- The ratio of neutral to charged pions was chosen randomly 
according to the ideal DCC probability distribution of the neutral pion frac- 
tion, Eq. (10). The isospin fluctuations of pion production from DCCs then 
lead to multiplicity fluctuations of charged pions at low transverse momenta 
and, therefore, to non-statistical fluctuations in M{pt)- For DCCs occurring 
in every event {F = 1) the fluctuations observed in the data rule out DCC 
sizes of ^ > 0.35. The limit on the amplitude of fluctuations in the pt parent 
distribution as a function of the frequency of occurrence of the events showing 
the fluctuation is shown in Fig. 42. 



3.4.3.3 RHIC experiments Among the four RHIC experiments, both 
STAR [23] and PHENIX [158] have the capability to look for DCC-type fluc- 
tuations. So far the PHENIX Collaboration have reported [22] some results 
on search for Centauro and anti-Centauro type of events at RHIC energies. 
They use the photon information from the electromagnetic calorimeter and 
charged-track information using drift chamber and pad chamber detectors. 
One such typical possible Centauro event is shown in Fig. 43. The analysis 
is still in preliminary stage. In addition to looking for such events, they also 
carry out a multi-resolution analysis, with input being the asymmetry between 
the number of charged tracks and neutral clusters, on an event-by-event basis 
as a function of subdivided r] — (p phase-space. 
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Fig. 43. Centauro-type event seen in the PHENIX Experiment at RHIC. The soHd 
circles are photon hits and plus signs are charged-particle hits. The box shows the 
Centauro region in r] — (p phase-space. From Ref. [22]. 



3.5 Possibility and limitations of DCC search 



For detecting a DCC signal with very low probability of occurrence, it is nec- 
essary to understand, from an experimental point of view, all the factors that 
affect the detection of DCC domain during its transition from the time of for- 
mation to the time of detection. This has been investigated in Ref. [141] using 
the simple model for simulating DCC events described earlier. The analysis 
is carried out using the method based on the DWT technique and looking at 
the signal-to-background ratio. Possible factors that can affect the DCC signal 
are: effect of the 7r° decay; the presence of multiple DCC domains; detector 
related effects, such as efficiency and purity of particle detection; use (or lack) 
of pr-information of the detected particles; increase in particle multiplicity as 
one goes from SPS to RHIC and LHC energies. 



3. 5. 1 Effect of decay of vr" 

The formation of DCC leads to large event-by-event fluctuation in the neutral 
pion fraction /. To observe this in an ideal situation one has to count the 
number of neutral pions 7r° event-by-event. But the tt^'s decay into photons 
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Fig. 44. Typical / distributions inside a DCC domain before (solid line) and after 
decay (dotted line) of (/")• Also shown for comparison is the / distribution for 
normal events (dashed line). Prom Ref. [141]. 

by the time they reach the detector. As a result, photons coming from DCC 
7r°'s inside a given DCC domain may move out of the rj — (f) phase-space of 
the domain - depending on their momentum - before they are detected. Also, 
photons from outside the DCC domain may enter the 77 — phase-space of the 
domain. Both effects result in dilution of the strength of the signal. Considering 
the effect of 7r° decay, one modifies Eq. (181) as: 



where 8^ is the resultant number of photons that get removed from or added 
into the DCC domain and A^ch is the number of charged pions in the region of 
T] — phase-space under consideration. For multiplicity detectors (e.g. SPMD 
in WA98, or FMD in ALICE), which do not have the capability of particle 
identification, A^'ch denotes the multiplicity of all charged particles in the con- 
sidered phase-space region. Since 5^ is a non-zero number, the possibility of 
observing a pure Centauro-type event (/ <^ 1) is low. In fact the distribution 
of the ratio /" will be shifted away from zero, as illustrated in Fig. 44, which 
shows the probability distributions for / and /" inside a DCC domain in the 
model calculation. Here Ach = A^7r± • The effect of detecting other charged par- 
ticles along with pions will be discussed later. Clearly one can notice that, as 
a result of 7r° decay, the peak of the original DCC distribution is shifted to a 
non-zero value ~ 0.15. Also shown in the figure is the / distribution for normal 
events. In that case, the / and /" distributions are not very different because 
the relative population of 7r° is the same within and outside the domain and 
hence the loss due to decay is compensated by the gain due to decay from 
other phase-space regions. 



A'^/2 ± 5, 
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A"^/2 ± 5^ + A"ch ' 
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Fig. 45. Percentage decrease in the strength of DCC signal due to the effect of decay 
of tt" as a function of domain size Acp. From Ref. [141]. 



In order to quantify the decrease in the detectability of the DCC signal, the 
authors of [141] considered the following two cases: 

(i) Analysis using a 7r° detector: 

In the first case one puts two hypothetical detectors, one for detecting 7r*^'s 
and other for detecting charged pions, with 100% efficiency, within one unit 
of rj and full coverage. The detector effects will be discussed later taking 
reahstic efficiency and other parameters. The goal here is to study the 7r° 
decay effect only. Then one introduces DCC domains in the simulated events 
with domain size A0 varying from 30° to 180°. These events without 7r° 
decay and having DCC domain of a particular size are then analyzed using 
the DWT method to obtain the FFC distribution. The strength parameter 
C is calculated from Eq. (183). 

(ii) Analysis using an ideal photon detector: 

In the second case one replaces the 7r° detector with an ideal photon de- 
tector having 100% detection efficiency and no contamination from charged 
particles. The 7r°'s are allowed to decay and a similar analysis as mentioned 
above is carried out to obtain the strength values. 

From the two such obtained strength values for each DCC domain size, one 
calculates the percentage decrease in strength value of the second case (ii) 
compared to the first case (i). This is shown in Fig. 45 for various domain 
sizes in Ac/). The results show that the decrease in strength of the signal is 
more for smaller domains of DCC as compared to larger ones. We mention 
that for each DCC domain size, two cases have been studied: one in which 
all the events have DCC-typc fluctuations and another where only 20% of the 
events have DCC-type fluctuations introduced. It is found that the results are 
similar for both cases. 
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3.5.2 Multiple domains of DCC 

In heavy-ion collisions there is a possibility that more than one domain of 
DCC be formed, each with a chiral condensate pointing in independent di- 
rections in isospin space. Of course, the contribution to the neutral fraction 
from different independent domains is expected to reduce the DCC signal in 
isospin fluctuations [159]. In particular, in case of a large number of indepen- 
dent domains of equivalent size, one expects the distribution of the total 
neutral fraction / to approach a Gaussian centered at (/) = 1/3 and of width 
inversely proportional to the square root of the number of domains, as per the 
central limit theorem. This makes the dis-entanglement of DCC signal from 
that of normal events difficult. 

These effects have been studied in [141] using the multi-resolution DWT 
method. For simplicity, and in view of limitations of the model in terms of 
introducing DCC at the freeze-out stage, these authors have placed varying 
number of DCC domains in each event, depending on the domain size in A0. 
All the domains are assumed to extend over one unit in pseudo-rapidity. For 
example, wc can place up to four domains of A(j) = 90°, without any overlap 
among themselves. These authors have also carried out the analysis by placing 
up to six DCC domains with A0 = 30° in each event. Care has been taken in 
placing the domains randomly so that no two domains overlap. The / value of 
each domain in each event is randomly chosen following the 1/2 vT probability 
distribution. The aim is to sec how the signal changes in terms of the strength 
parameter as given in Eq. (183). 

Figure 46 shows the variation of strength of the DCC signal with the num- 
ber of domains in a given event for two domain sizes of Acf) — 90° and 30°. 
One observes that the strength of the DCC signal increases as the number of 
domains increases and saturates for larger number of domains. The increase 
in strength value with increase in number of domains is because the multi- 
resolution event-by-event analysis based on the DWT method has been able to 
pick up signals by looking at bin to bin (in (p) fluctuations in each event. It is 
worth emphasizing that this result is not in contradiction with the theoretical 
considerations described above [159,160]. Indeed, if the multiple DCC domains 
formed in the initial stages, in the course of their evolution, move towards the 
same part of phase space covered by a detector, then the strength of the signal 
will reduce. However, if the phase space separation is maintained in evolving 
DCC domains, then the strength will increase as found in the present analysis. 



Clearly, the situation is better if one of the domains is predominant [160]. 
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Fig. 46. Variation of the strength of DCC signal with increase in number of DCC 
domains in an event. All the events are DCC events. From Ref. [141]. 

3.5.3 Detector effects 



The detector limitations are basically those related to efficiency and purity of 
particle detection and to the acceptance in r/ — 0. The effect of efficiency of 
particle detection is trivial, higher the efficiency of particle detection, more 
reliable is the search for DCC. Similarly the role of higher acceptance in rj 
and full azimuthal coverage can be hardly over-emphasized (lower acceptance 
effectively reduces the total multiplicity of observable DCC pions and hence 
reduces the strength of the DCC signal). If the acceptance in t] is sizable, one 
can attempt DWT analysis using bins in 1] for limited azimuthal coverage, 
which would be complementary to the various analysis reviewed here. Such an 
analysis would certainly be of great interest. 

The purity of the photon sample, as measured in a detector like PMD, is 
around 60% — 70%, whereas that of charged particle samples measured using 
detectors hke FTPC or FMD is quite high, larger than 95%. Hence we discuss 
here only the effect of purity of photon sample. 



3.5.3.1 Effect of neutral pion fraction on the purity of photon sam- 
ple It is of interest to study the effect of the DCC neutral pion fraction on 
the purity of the photon sample. Consider the case where the neutral pion 

fraction for a DCC domain is / = 1. This corresponds to a situation where 
there are no charged pions. This will lead to greatly reduced number of charged 
particles falling on the photon detector. So the purity of the photon sample 
will be high. The variation of purity with neutral pion fraction can be investi- 
gated in the following manner [141] : The number of 7- like hits on the photon 
detector can be written as 

Ar^_iike = eA^^ + cA^ch , (196) 
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Fig. 47. Variation of purity of the photon sample as a function of neutral pion 
fraction. For / = 0.33 (generic event) purity is 0.74. From Rcf. [141]. 

where the first term denotes the contribution of actual photons and the sec- 
ond term denotes charged-particle contamination. The factor c is the fraction 
of charged-particles on the PMD acceptance treated as contamination to the 
detected photon sample. One assumes that the latter is given by a normal 
distribution with mean of 15% and width of 5%, the percentage being taken 
with respect to the total charged particles within the acceptance of the photon 
detector. This is a reasonable number, considering that the converter thick- 
ness of the preshower detector is ~ 10% of an interaction length and some 
interactions lead to multiple clusters (the effect of overlapping clusters is ig- 
nored). DCC-type fluctuations are introduced for fixed / values and domain 
size A(f) = 360° and Arj = 1 in a set of event. Several such sets of events 
were generated for / values varying from 0.025 to 1.0. The purity of the pho- 
ton sample is then calculated for each set of events with a fixed neutral pion 
fraction value. 

The results are shown in Fig. 47, where one clearly sees that the purity of 
photon sample increases with the increase in the neutral pion fraction value. 
However the purity does not reach a value of 1 for / = 1, because there 
are charged particle other than vr^ falling on the detector. Similarly it never 
reaches a value of zero for / = 0, as there are 7's from tt^'s decaying outside 
the detector acceptance and from other sources. 



3.5.3.2 Effect of charged-particle contamination in PMD Photon 
measurement in heavy-ion collisions has the problem of charged-particle con- 
tamination. The charged particles detected as photons have a correlation 
with those detected in the charged-particle multiplicity detector. This ad- 
ditional correlation suppresses the possible anti-correlation between photons 
and charged particles that may arise due to DCC formation. 
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Fig. 48. Variation of rms deviation of FFC distribution for pure VENUS events with 
increase in the percentage of contamination in the photon detector. The variation 
of rms deviation of FFC distributions of mixed events constructed form VENUS are 
also shown. From Ref. [141]. 

The effect of such a correlation has been studied in [141] with the help of 
VENUS (no DCC) events and a set of mixed events generated from them as 
described above. The percentage of charged-particle contamination was var- 
ied from 0% to 50%. The event-by-event A^-y-ukc was kept the same for all the 
cases. The result of the DWT analysis on these events is shown in Fig. 48. 
There is a decrease in the rms deviation of the FFC distribution for pure 
VENUS events with increase in the percentage contamination in the photon 
sample. Thus the presence of any additional anti-correlation due to DCC-like 
effect will have to overcome the opposing correlation effect due to contamina- 
tion in order to be observed. Figure 48 also shows the effect of contamination 
on Ml-type mixed events generated from the VENUS events. It is found that 
the rms deviations of the FFC distributions for the mixed events remain inde- 
pendent of the level of contamination. The mixed events appropriately break 
the additional correlation due to charged-particle contamination. 

The effect of charged-particle contamination can be corrected by knowing the 
level of contamination in photons for a given detector. Figure 49 shows the 
variation of rms deviation of FFC distribution for DCC-like events with in- 
crease in percentage of events being DCC-type in a given ensemble of nDCC 
events, described previously. It also shows the variation of rms deviation of 
FFC distribution of corresponding mixed events. The rms deviations of the 
FFC distribution of mixed events are found to be independent of the per- 
centage of events being DCC-type. This is along the expected lines. But for 
lower fraction of events being DCC-type it is above that of the parent sam- 
ple of nDCC events. This is because the A'^-uke — -^ch anti-correlation due 
to DCC-type effect is not sufficient to overcome the correlation between A'ch 



107 



0.25 



nDCC 



Mixed 

Mixed - corr 

I I 

1 10 10^ 

% of events being DCC type 

Fig. 49. Variation of rms deviation of FFC distribution of DCC-type events as a 
function of the percentage of events being DCC-type. Also shown are the results for 
corresponding mixed events and mixed events corrected for contamination effect. 
From Ref. [141]. 

and charged particles detected as contamination in the photon detector. How- 
ever as the fraction of events being DCC-type increases, the DCC-type effect 

dominates. The effect of the contamination can be corrected. This is done by 
taking into account the difference in rms deviation of normal events and mixed 
events. The rms deviation of FFC distribution of corrected mixed events are 
shown in Fig. 49. 

The results discussed in the above two sections indicate that it is better to 
look for anti-Centauro (photon excess) events in studies using photon and 
charged-particle multiplicity detectors. This is because the effect of decay is 
primarily to shrink the / distribution from the lower / side and high value of 
/ leads to a reduced effect of charged particle contamination on DCC search. 
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3.5.4 Effect of low-pT -information of particles 

It is expected that pr-information of particles would be very helpful in DCC 
search. This would also enable one to verify various expected features of DCC 
production, such as the fact that DCC-pions have low transverse momenta. 
In order to show the utility of p^^-information, the following study has been 
carried out in Ref. [141]. One assumes that all charged particles with px greater 
than 50 MeV/c are detected (the pr acceptance of preshower detector like 
PMD extends down to 30 MeV/c [144]). The pr-resolution for the charged- 
particle detector was taken as Apt/pt = 0.2 [161]. These realistic parameters 
correspond to a typical experiment for DCC search, as in STAR. DCC events 
are simulated as described above. Since DCC pions are believed to have low pt, 
one introduces DCC-type fluctuations in pions with pt < 150 MeV/c. Then 
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Fig. 50. Variation of C as a function of the fraction of events being DCC-type, 
for charged particle detector with pT-information and photon detector without 
pr-information (dotted Une) , for photon detector with pT-information and charged 
particle detector without j^y-information (dashed line) and for both charged parti- 
cle and photon detector with px-information (dot-dashed line). Also shown are the 
corresponding results with both charged particle detector and photon detector with- 
out pt information (solid line) . The horizontal straight line indicates the statistical 
error on C,. Prom Ref. [141]. 

a DWT analysis is carried out. The sample function in the DWT analysis 
is modified such that iVch taken corresponds to charged pions having pr < 
150 MeV/c. 

Figure 50 shows the strength C, of the DCC signal as a function of the percent- 
age of events being DCC-type. For comparison, also shown is the corresponding 
C, where the charged particle detector has no p^-information. We note that the 
absolute value of C, is lower compared to those presented in previous sections. 
This is because here we have taken only low-px pions as DCC pions in a given 
domain, whereas earlier all the pions in the domain were considered to be 
DCC-type. This reduces the strength as a result of reduction in multiplicity. 
The statistical error on C, for the present case is 0.07. 

Clearly one can see the increase in strength of the signal with the use of 
PT^-information. A similar analysis was carried out assuming that all photons 
have p2-information while the charged particles do not have any. The re- 
sults are also shown in Fig. 50. One sees that, although there is an increase 
in signal strength compared to the case with no pr-cut, it is much less com- 
pared to the case with p-r-information for charged particles. For the case where 
the analysis is carried out assuming that both photons and charged particles 
have pT-information, one finds that the results are close to those obtained for 
charged-particle detector with pT-information, except for events having higher 
percentage of DCC type events. 
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Fig. 51. Resultant pt distribution of charged particles with 20% of the events being 
DCC-type, obtained by subtracting the pT distribution of charged particles within 
2 o" of the FFC distribution from those beyond. Two cases are presented, additional 
low PT DCC pions being 60 and 100. The statistical errors are within the symbol 
size. From Ref. [141]. 

It is also interesting to look whether one can detect the enhancement of low- 
Pt pions due to the presence of a DCC. This can be done if the transverse 
momentum of charged particles is measured. In order to incorporate this effect 
in simulation, a number of low-pT pions {pr < 150 McV/c), with neutral 
fraction distributed according to the ideal DCC inverse square root law and 
having uniform px distribution were added within the chosen domain on top 
of the existing pions of a normal event [139]. The number of pions to be 
added depends on the size and energy density of the DCC domain. Two cases, 
with 100 and 60 additional low-p^ pions were considered, corresponding to 
an energy density within a DCC domain of about 50 MeV/fm^ and a domain 
radius of the order of 3 — 4 fm. As a test case a event sample with 20% of the 
events having a DCC domain of size A0 = 90° and Arj — 1 was generated. 

The event sample has been analyzed by the DWT method to obtain the FFC 
distribution, which was found to be near Gaussian for the present case [141]. 
One divides the FFC distribution into two parts: one within ±2a of the mean 
and the other beyond this. One then obtains the separate pt distributions of 
these two subsets of events. Finally, the histogram corresponding to events 
within 2cr of the FFC distribution is subtracted from those beyond, after 
proper normalization. The resultant spectrum is shown in Fig. 51. Clearly one 
sees that the above analysis technique should provide a useful tool to able to 
observe the low-pr pions enhancement due to possible DCC formation. 
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Fig. 52. Variation in the strength of DCC signal with multipHcity in a domain of 
A(f) = 90° and At] = 1. Here, is the average of the sum of mean photon and mean 
charged particles in a set of events. All the events are DCC events. Prom Ref. [141]. 

3. 5. 5 Effect of increase in center of mass energy 

The effect of higher multiplicity is to reduce the event-by-event statistical 
fluctuation associated with the multiplicity of photons and charged particles. 
A simple estimate of the strength of DCC-type fluctuations shows that for a 
DCC domain of azimuthal size A0 = 90° and A?7 = 1, the strength of DCC 
signal increases with increase in average multiplicity of photons and charged 
particles. Figure 52 shows the variation of the strength ^ of the DCC signal 
with increasing multiplicity, where N is the average of the sum of mean photon 
and mean charged particles in a set of events for one unit in r]. One observes 
that the increase is almost like a sfN effect. This is probably a consequence 
of going from a narrow Gaussian distribution to a wider 1/2^/f distribution. 
It may be mentioned that typical values of N in central collisions for Ary ^ 1, 
for the WA98 Experiment at SPS (combination of PMD and SPMD detectors 
for charged-particle detection) was about 320. It is expected to be about 500 
for the STAR Experiment at RHIC (PMD + FTPC) and about 2800 at the 
LHC, with the ALICE detector (PMD + FMD). 



4 Summciry and conclusions 



Since it has been proposed in the early 1990's, the idea that a DCC might 
be formed in high-energy collision experiments has triggered an intense activ- 
ity, both theoretical and experimental. Attempts to detect this phenomenon 
include various cosmic ray experiments, the T-864 MiniMAX experiment - 

This assumes a pseudo-rapidity density of ~ 8000 at mid-rapidity [141]. 
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a dedicated search in proton-antiproton collisions at the Fermilab Tevatron, 
or event-by-event analysis of Pb-Pb collisions data by the WA98 and NA49 
collaborations at the CERN SPS. No clear experimental evidence for DCC 
formation has been reported so far. This is a rather disappointing result in 
view of the fact that the phenomenon seems quite natural from a theoretical 
point of view. In particular, the growth of long- wavelength fluctuations due to 
a sudden cooling is a generic phenomenon, which, for instance, is commonly 
observed in ferromagnets in condensed matter physics. Numerical simulations 
of the non-equilibrium chiral phase transition in heavy-ion collisions suggest 
that this phenomenon might actually occur in existing or planned experiments 
at RHIC and LHC. 

The absence of experimental evidence may indicate that the appropriate con- 
ditions for DCC formation have not been met yet. For instance, it has been 
suggested that non-central collisions of large nuclei could be a more interest- 
ing place than central ones to look for DCC signatures. Another interesting 
suggestion is to look at very forward rapidities of high energy hadronic or 
nuclear collisions. But it should also be stressed that the present experimen- 
tal situation is actually consistent with present theoretical expectations: The 
estimated upper bound for the probability of DCC formation in heavy-ion 
collisions at SPS energies is at the edge of the current experimental limit. 
Moreover, present theoretical calculations suggest that the phenomenon we 
are seeking may actually be more difficult to observe than originally thought. 
In particular, numerical simulations of the chiral phase transition indicate that 
the DCC might not have the ideal isospin structure originally proposed. This 
should be taken into account in future theoretical investigations as well as 
experimental searches at RHIC and LHC. 

Besides providing a useful tool for the study of the chiral structure of QCD, 
the detection of a DCC is of fundamental interest on its own. Indeed, the 
possible formation of a coherent, classical pion field is to be expected on very 
general grounds as a direct consequence of the bosonic nature of pions. In 
analogy with the - nowadays commonly observed - phenomenon of Bose- 
Einstein condensation for non-relativistic systems, the DCC can be thought of 
as a macroscopic wave of QCD matter. The possibility that the latter might 
indeed be produced in existing experiments is very exciting. It is certainly 
worth have a look! 
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